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FOREWORD 


This  report,  which  presents  the  experimental  and  theoretical  results  of 
a  program  of  supersonic  flutter  testing,  was  prepared  by  the  Aeroelastic  and 
Structures  Research  Laboratory,  Massachusetts  Institute  of  Technology, 

Cambridge  39,  Massachusetts  for  the  Aircraft  Laboratory,  Wright  Air  Development 
Center,  Wright-Patterson  Air  Force  Base,  Ohio.  The  work  was  performed  at  the 
MIT  under  the  direction  of  Professor  R.  L.  Halfman,  and  the  project  was 

supervised  by  Mr.  J.  F.  McCarthy,  Jr.  The  research  and  development  work  was 

accomplished  under  Air  Force  Contract  No.  AF  33 (038) -22955,  Project  No.  1370, 
(Unclassified  Title)  "Aeroelastici ty,  Vibration  and  Noise,"  and  Task  No. 

13474,  (Unclassified  Title)  "Three-Dimensional  Supersonic  Flutter  Model  Tests 
at  Mach  Number  1.5".  'Mr.  Niles  R.  Hoffman  of  the  Dynamics  Branch,  Aircraft 
Laboratory,  is  task  engineer.  This  task  covers  a  continuing  effort  on  flutter 
research  at  supersonic  speeds.  Research  was  started  in  March  1951.  The  test 
data  presented  in  this  report  was  obtained  during  the  period  from  December 
1952  to  December  1954.  This  is  Part  1 1  of  a  report  to  be  issued  in  three 

separate  parts.  Part  I  of  this  report,  V/ADC  TR  54-113,  (Unclassified  Title) 

"Three-Dimensional  Supersonic  Flutter  Model  Tests  Near  Mach  Number  1.5,  Part  I. 
Model  Design  and  Testing  Techniques,"  was  issued  in  December  1955.  Technical 
Report  V/ADC  TR  54-114,  "(Unclassified  Title)"  A  Variable  Mach  Number  Supersonic 
Test  Section  for  Flutter  Research"  was  issued  in  December  1954. 

The  authors  are  indebted  to  Mr.  0.  V/atlin,  and  Mr.  C.  Fall  for  their 
help  in  the  model  construction  and  in  keeping  the  wind  tunnel  in  operation;  to 
Mr.  6.  M.  Falla  for  the  photograph;  to  Messrs.  A.  Heller  and  H.  Hagerup  for 
their  help  in  the  calculations;  to  Messrs.  J.  R.  Friery,  G.  Anitole,  and 
W.  Marchant  for  their  help  in  preparing  the  tables  and  figures;  and  to  Miss 
K.  Roberts  and  Mrs.  B.  Marks  for  their  help  in  typing  this  report. 

This  document  is  classified  CONFIDENTIAL  in  its  entirety  (excepting  the 
title)  because  results  of  supersonic  flutter  tests  generally  indicate  limiting 
performance  capabilities  of  present  and  future  military  aircraft  and  have 
application  in  the  form  of  design  criteria. 
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ABSTRACT 

Three-dimensional  supersonic  flutter  tests  were  made  on 
over  75  semi-span  models  in  the  MIT  blowdown  wind  tunnel  facil¬ 
ity.  The  testing  technique  involved  injecting  the  model  into  a 
stable  region  of  flow  and  decreasing  the  Mach  number  until  in¬ 
stability  occurred.  Experimental  flutter  stability  boundaries 
are  defined  for  bare  straight,  swept  and  delta-wing  planforms 
in  the  Mach  number  range,  1.3  -  2.0.  Exploratory  tests  were 
also  made  on  wings  with  ailerons  for  all  planforms,  and  on 
straight  and  swept  wings  with  tip  tanks  for  both  cantilever  and 
free-to-roll  root  conditions.  Except  for  absolute  stiffness, 
the  dimensionless  flutter  parameters  were  chosen  so  as  to  be 
typical  of  present-day  high-speed  aircraft. 

Extensive  theoretical  calculations  were  made  on  the 
straight -wing  planform  using  two-dimensional  supersonic  oscilla 
tory  aerodynamic  coefficients  and  three-dimensional  structural 
properties.  The  qualitative  prediction  by  the  theory  of  the 
effect  of  various  parameter  changes  generally  agrees  with  ex¬ 
periment,  but  the  quantitative  prediction  is  generally  poor. 

The  theoretical  calculations  are  unconservative  in  that  they 
predict  smaller  regions  of  instability  than  those  obtained  ex¬ 
perimentally  at  Mach  numbers  above  1.4.  No  theoretical  calcula 
tions  were  made  for  the  swept  and  delta  planforms. 

Comparison  of  the  experimental  data  with  the  results  of 
other  flutter  tests  shows  that  in  the  Mach  number  range  of  0.6 
to  2.0  for  wings  with  parameters  similar  to  those  tested  the 
following  conclusions  may  be  drawn: 
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1) .  At  constant  altitude, the  critical  flutter  region 

for  straight  wings  lies  at  about  M  =  1.7 

2) .  At  constant  altitude,  the  critical  flutter  region 

for  swept  wings  lies  in  the  transonic  regime  at 
about  M  s'  1 . 1 

3) .  At  constant  altitude,  the  critical  flutter  region 

for  6o°  delta  wings  lies  at  the  highest  Mach  number 
tested,  M  =  2.0 

4) .  At  constant  dynamic  pressure,  the  critical  flutter 

region  lies  in  the  transonic  regime  close  to 
M  =  1.0,  for  all  the  straight,  swept,  and  delta 
planforms  tested 

A  complete  tabulation  of  the  design  properties  for  all 
the  models  tested  is  presented  along  with  the  results  of  static, 
vibration,  and  flutter  tests. 


PUBLICATION  REVIEW 


This  report  has  been  reviewed  and  is  approved. 

FOR  THE  COMMANDER: 

DANIEL  D.  McKEE 
Colonel,  USAF 

Chief,  Aircraft  Laboratory 
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Static  unbalance  of  tip  tank  about  elastic  axis 
of  wing,  positive  nose  up 


Dimensionless  static  unbalance  of  tip  tank, 

rr\  b-£ 

Static  unbalance  per  unit  span  of  wing  section 
about  the  elastic  axis,  positive  nose  up 


Static  unbalance  per  unit  span  of  aileron  about 
the  aileron  hinge  line 


Absolute  temperature 
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LIST  OF  SYMBOLS  (Contd.) 


Time 

Volume  of  tip  tank 
Free-stream  velocity 
Work 

Coordinate  measured  perpendicular  to  elastic 
axis  from  mid-chord  of  airfoil,  positive  aft 

Fraction  of  chord  elastic  axis  of  wing  is 
behind  leading  edge 

Dimensionless  distance  the  center  of  gravity 
of  the  wing  section  lies  aft  of  the  elastic 
axis , 

*  W\  i-i 


Dimensionless  distance  the  center  of  gravity 
of  the  aileron  lies  aft  of  the  aileron  hinge 


Coordinate  measured  parallel  to  elastic  axis 
from  root  chord  of  right  wing 

Frequency  ratio,  (Toy 

Vertical  displacement  of  wing,  positive  down 

Twist  about  elastic  axis,  positve  nose  up 

Amplitude  of  assumed  torsion  mode  at  wing  tip 
station 

Rotation  of  aileron  about  its  hinge  line, 
positive  for  increasing  angle  of  attack 

Ratio  of  specific  heats 

Dimensionless  spanwise  variable,  -|— 
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Subscripts 

A 

BA 
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LIST  OF  SYMBOLS  (Contd.) 

Rigid-body  angular  rotation  about  wing  root, 
positive  right  wing  down 

Amplitude  of  assumed  0  -mode  at  wing  tip  station 

Wing  taper  ratio,  ratio  of  wing  tip  chord  to 
wing  root  chord 

Wing  mass-density  ratio, 

Wing  mass-density  ratio, 

Kinematic  viscosity 
Free-stream  air  density 
Density  of  balsa  wood 
Thickness  ratio  of  airfoil 
Frequency  ratio,  (^7) 2  (1  +  ig) 

Flutter  frequency 

?  W.  K1 

Reduced  frequency,  parameter, 

First  bending  frequency  of  bare  wing 
Second  bending  frequency  of  bare  wing 
First  torsional  frequency  of  bare  wing 
Second  torsional  frequency  of  bare  wing 
First  aileron  frequency  of  wing 


w* 

TT/ofe*- 


xyo 


VY1 


Airplane 
Balsa  Wood 
Wing  chord 
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Center  of  gravity 
Flutter 

First  bending  mode 

Model 

Wing  root 

Refers  to  stagnation  conditions  in  wind  tunnel 
Tip  tank 

First  torsional  mode 
3/4  span  station 
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SECTION  I 
INTRODUCTION 

The  phenomenon  of  flutter,  the  self-excited  oscillation  of 
an  elastic  structure  in  an  airstream,  appeared  early  in  the  de¬ 
velopment  of  aircraft.  At  first,  designers  could  apply  only 
crude  corrective  measures,  and  it  was  not  until  the  early  1930's 
that  aeronautical  scientists  were  able  to  fashion  promising  theo¬ 
retical  approaches  to  the  problem  of  flutter.  Although  experi¬ 
mental  work  lagged  behind  the  theoretical  attack,  the  flutter 
problem  in  incompressible  flow  had  become  quite  tractable  by  the 
end  of  the  Second  World  War.  Both  the  theoretical  and  experi¬ 
mental  approaches  were  well  developed  and  understood.  Until  the 
advent  of  transonic  arid  supersonic  aircraft,  these  techniques 
were  adequate  for  the  airplane  designer. 

Unfortunately,  the  trend  towards  higher  speeds,  increased 
structural  flexibility  and  lower  aspect-ratio  lifting  surfaces  as 
well  as  the  growing  use  of  large  external  stores  so  aggravate  the 
flutter  problem  that  it  is  now  often  a  primary  design  considera¬ 
tion  rather  than  an  occurrence  that  can  be  remedied  fairly  easily. 
Because  of  the  questionable  reliability  of  existing  methods  of 
flutter  analysis  in  the  high-speed  range,  the  airplane  designer 
calls  upon  the  experimentalist  to  provide  data  that  are  immedi¬ 
ately  useful  and  that  can  be  used  to  confirm  theory.  A  reason¬ 
able  amount  of  experimental  flutter  data  exists,  but  attempts  to 
correlate  this  data  with  calculations  made  with  incompressible 
unsteady  aerodynamic  coefficients  have  not  been  too  successful 
(see  Ref .  4) . 

In  the  low  supersonic  speed  range,  M  =  1.2  to  M  =  2.0, 
basic  flutter  theory  may  be  used  successfully  with  linearized 
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aerodynamic  forces  provided  the  perturbation  velocities  on  the 
system  to  be  studied  are  small  compared  to  the  free  stream 
velocity.  Flutter  analyses  using  two-dimensional  supersonic 
aerodynamic  coefficients  (Ref.  7)  have  been  done  for  some  time, 
but  there  has  been  little  confirmation  of  the  theory  with  ex¬ 
periment  over  a  range  of  Mach  numbers  and  for  dimensional  para¬ 
meters  which  are  typical  of  present-day  aircraft,  especially  in 
the  case  of  the  mass  ratio,  which  has  generally  been  much  higher 
for  the  models  tested  than  that  encountered  in  practice 
(Refs.  8-11). 

Experimental  investigations  of  flutter  in  the  supersonic 
speed  range  at  first  glance  appear  to  be  as  difficult  as  those 
in  the  transonic  range.  Models  mounted  on  rockets,  bombs,  or 
sleds  must  go  through  the  transonic  range  before  encountering 
supersonic  speeds  so  that  the  same  problems  of  complexity  and 
expense  (mounting,  data  recording,  expendable  models,  etc.)  are 
still  present.  For  tests  conducted  in  the  wind  tunnel,  there 
are  the  large-power  requirements  for  continuous -flow  tunnels  of 
reasonable  size,  the  aerodynamic  problems  of  obtaining  uniform 
flow  in  the  test  section,  and  the  potential  damage  that  could 
be  inflicted  on  the  testing  facility  by  models  which  are  lost. 
For  supersonic  flow,  the  Mach  number  in  the  test  section  is  a 
function  only  of  the  geometry  of  the  nozzle,  so  it  might  appear 
that  many  nozzles  of  fixed  geometry  would  have  to  be  used  to  get 
useful  data.  The  model  designer  is  again  confronted  with  the 
problem  of  building  efficient  structure  into  thin  wings  in 
order  to  obtain  the  very  high  natural  frequencies  required  when 
Mach  number  must  be  simulated  (see  Reference  12). 

Some  of  the  complexity  and  expense  of  testing  flutter 
models  in  the  low  supersonic  speed  range  was  eliminated  by  de¬ 
signing  the  facility  described  in  Reference  13.  Briefly,  a 
supersonic  nozzle  of  variable  geometry  was  built  for  installa- 
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tion  in  a  blowdown  wind  tunnel.  It  is  of  the  asymmetric  sliding- 
block  type,  and  the  Mach  number  can  be  varied  through  the  com¬ 
plete  range  of  the  nozzle  during  a  run  (M=1.2  -  2.1)  without  any 
change  in  the  dimensions  of  the  test  section.  The  testing  tech¬ 
nique  involves  injection  of  the  model  into  the  air stream  in 
order  to  avoid  destruction  of  the  model  by  the  violent  starting 
shock.  Thus,  the  problems  of  large  power,  damage  from  broken 
models,  and  testing  at  fixed  Mach  number  were  immediately  solved. 
Also,  the  difficulty  of  obtaining  low  mass  ratios  is  somewhat 
alleviated  because  of  the  high  air  density  in  the  test  section, 
which  is  characteristic  of  a  blowdown  wind  tunnel.  The  approach¬ 
es  for  designing,  building,  and  testing  inexpensive  supersonic 
flutter  models  with  desired  parameters  are  discussed  in  Refer¬ 
ences  14-1?. 

The  planforms  shown  in  Figure  1.1  were  chosen  for  investi¬ 
gation  (Ref.  18).  These  are  typical  of  present-day  high-speed 
fighters  and  proposed  supersonic  bombers.  Although  most  of  the 
work  was  done  on  the  bare  wings,  some  exploratory  tests  were 
made  on  models  with  ailerons  and,  in  the  case  of  the  straight 
and  swept  wings,  on  models  with  tip  tanks  for  both  cantilever 
and  free-to-roll  root  conditions.  The  10-inch  root  chord  rep¬ 
resents  the  model  of  maximum  size  that  can  be  tested  in  the 
facility  without  shock  interference  (see  Reference  14).  The 
symmetrical  double  wedge  airfoil  section  was  used  because  of  its 
simplicity.  The  range  of  flutter  parameters  built  into  the 
models  are  typical  of  current  high-speed  aircraft. 

Theoretical  flutter  calculations  were  made  only  for  the 
straight-wing,  since  it  was  felt  that  two  dimensional  aerody¬ 
namic  coefficients  could  be  used  successfully  for  this  planform. 
No  theoretical  work  was  done  on  the  swept  or  delta  planforms 
because  the  primary  emphasis  of  the  research  program  was  on  ex¬ 
perimental  results.  However,  enough  theoretical  work  was  done 
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on  the  straight  wing  so  that  the  complete  range  of  experimental 
parameters  was  covered  and  a  reasonable  comparison  between 
theory  and  experiment  could  be  made  for  this  planform.  The 
theoretical  trends  exhibited  by  the  straight  wing  were  compared 
with  those  obtained  experimentally  for  the  swept  and  delta  wings. 

In  subsonic  and  low  transonic  flutter  testing  wings  are 
brought  from  a  stable  to  an  unstable  region  by  increasing  the 
speed  or  Mach  number.  Figure  1.2,  which  shows  a  hypothetical 
flutter  boundary,  demonstrates  a  peculiarity  of  supersonic  flut¬ 
ter  testing.  A  wing  is  generally  brought  from  a  stable  region 
to  an  unstable  region  (see  Reference  14)  by  decreasing  the  speed 
or  Mach  number  along  a  tunnel  operating  curve.  From  Fig.  1.2 


M 

FIGURE  1.2  HYPOTHETICAL  FLUTTER  BOUNDARY 

it  can  also  be  seen  that  decreasing  the  region  of  instability 
over  the  whole  Mach  number  range  will  lower  the  Mach  number  for 
flutter  in  the  supersonic  speed  range  and  raise  the  Mach  number 
for  flutter  in  the  subsonic  and  low  transonic  speed  range. 
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SECTION  II 
FORMULATION 

The  dimensionless  answers  which  one  hopes  to  obtain  from 
a  flutter  model  test  are  the  values  of  the  Mach  number,  M,  the 
reduced  frequency,  k,  and  the  frequency  ratio,  Z,  at  flutter  and 
the  flutter  mode  shape.  These  dimensionless  quantities  can  be 
derived  from  basic  flutter  theory  (see  Reference  19  and  Appendix 
A.l)  along  with  other  parameters  which  define  the  physical  prop¬ 
erties  of  the  lifting  surface.  For  bare  wings  of  reasonably 
large  aspect  ratio,  the  physical  properties  evolve  as  dimension¬ 
less  parameters  which  must  be  defined  at  every  spanwise  station, 
viz. , 


planform,  ^ 

o 

mass  distribution,  — 

mo 

location  of  the  chordwise  center  of  gravity,  x^ 

dimensionless  moment  of  inertia  in  pitch,  r^ 

El 

bending  stiffness  distribution, 

GJ 

torsional  rigidity  distribution, 

location  of  the  elastic  axis,  a 

For  wings  of  low  aspect  ratio  or  for  surfaces  where  chord- 
wise  deformations  are  appreciable,  the  concepts  of  bending  stiff¬ 
ness,  torsional  rigidity  and  elastic  axis  are  not  valid,  and  the 
chordwise  as  well  as  the  spanwise  distributions  of  mass  and  stiff¬ 
ness  must  be  defined. 
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For  given  distributions  of  mass  and  stiffness,  the  model 
has  a  discrete  set  of  eigenfrequencies  with  associated  eigen¬ 
functions  and  a  distinct  set  of  influence  coefficients.  For  both 
models  and  full-scale  airplanes,  these  latter  parameters  are  rel¬ 
atively  easy  to  obtain  experimentally  compared  to  obtaining  mass 
and  stiffness  distributions,  especially  when  the  chordwise  as 
well  as  the  spanwise  variations  must  be  considered.  Therefore, 
in  the  usual  formulation  of  the  flutter  problem,  natural  fre¬ 
quencies  and  mode  shapes  or  matrices  of  influence  coefficients 
replace  stiffness  distributions.  Through  the  years,  flutter 
engineers  have  become  used  to  thinking  of  the  flutter  problem 
in  terms  of  natural  frequencies  and  mode  shapes  rather  than  in 
terms  of  stiffnesses.  The  concept  of  attacking  the  problem  with 
matrices  of  influence  coefficients  is  relatively  new,  since 
lifting  surfaces  of  low  aspect  ratio  have  only  recently  become 
popular.  The  tendency  is  to  treat  this  latter  type  of  planform 
in  the  same  way,  i.e.,  in  terms  of  natural  frequencies  and  mode 
shapes  when  attacking  the  problem  physically  or  when  interpre¬ 
ting  results.  This  tendency  is  still  valid  since  coupled  fre¬ 
quencies  and  mode  shapes  evolve  as  theoretical  solutions  to  the 
flutter  equations  at  zero  airspeed,  even  though  the  problem  has 
been  formulated  in  terms  of  influence  coefficients. 

It  is  the  intention  of  this  report  to  adhere  to  the 
classical  concepts  of  natural  frequencies  and  mode  shapes  and 
elastic -axis  location,  rather  than  to  consider  matrices  of  in¬ 
fluence  coefficients  or  stiffness  distributions.  Valid  criticism 
may  accompany  the  concept  of  an  elastic  axis  for  any  of  the  plan- 
forms  considered  because  of  sweep  and  low  aspect  ratio.  Also, 
in  the  interest  of  simplicity,  all.  models  were  designed  to  have 
identical  spanwise  distributions  of  mass  and  stiffness.  The 
former  varies  as  the  square  of  the  chord,  and  the  latter  varies 
as  the  fourth  power  of  the  chord.  These  distributions,  which 
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were  accomplished  by  tapering  all  dimensions  linearly  from  root 
to  tip,  are  typical  of  present-day  aircraft.  The  absolute  value 
of  bending  and  torsional  stiffness  are  measured  by  the  magnitude 
of  the  first  bending  and  the  first  torsional  frequency,  respec¬ 
tively,  and  that  of  the  mass  by  the  value  of  the  flow  parameter, 
jx.  .  The  models  were  designed  to  have  constant  values  of 
chordwise  center-of -gravity  location,  dimensionless  moment  of 
inertia  in  pitch  and  elastic-axis  location  at  every  spanwise 
station.  No  attempt  was  made  to  design  chordwise  distributions 
of  mass  or  stiffness  into  the  models  although  matrices  of  in¬ 
fluence  coefficients  were  measured  in  some  cases. 

With  these  simplifications,  the  parameters  which  were  con¬ 
sidered  for  each  bare  wing  planform  are: 

mass  ratio,  ^ 

location  of  the  chordwise  center  of  gravity,  x^ 

dimensionless  moment  of  inertia  in  pitch,  r^ 

first  torsional  frequency, 

LOh 

frequency  ratio, 

location  of  the  elastic  axis,  a. 

The  ranges  of  the  values  of  these  parameters  were  chosen  to  be 
typical  of  present-day  aircraft.  By  virtue  of  the  design,  the 
quantities ,  /a.  ,  x^  ,  r^,  and  a  are  constant  at  all  spanwise 
stations;  the  frequencies  and  ,  are  three-dimensional 
structural  properties.  Natural  frequencies  higher  than  first 
bending  and  first  torsion  were  not  considered  separately  in  the 
model  design  since  their  values  are  dictated  by  the  choice  of 
the  mass  and  the  stiffness  distributions  already  mentioned.  We 
also  note  that  the  mass  and  stiffness  distributions  determine 
the  values  of  the  mode  shapes  associated  with  each  natural  fre- 
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quency. 


are: 


Other  parameters  which  evolve  from  basic  flutter  theory 


airfoil  shape 

structural-damping  coefficients 
Reynolds  number 
Prandtl  number 
ratio  of  specific  heats. 

For  simplicity,  the  airfoil  section  was  taken  as  a  symmetrical 
double  wedge  for  all  the  models,  since  it  is  probably  of  second¬ 
ary  importance  in  the  flutter  problem.  For  the  structural-damp¬ 
ing  coefficients,  the  correct  order  of  magnitude  was  obtained 
through  discriminate  choice  of  structural  material.  As  high  a 
Reynolds  number  as  possible  was  obtained  by  using  as  large  a 
scale  model  and  as  high  a  fluid  density  as  practical  with  the 
available  facility  (Reference  13).  Control  of  the  values  of 
Prandtl  number  and  the  ratio  of  specific  heats  was  not  consider¬ 
ed  since  air  was  used  as  the  testing  medium. 

For  wings  with  control  surfaces,  additional  parameters 
evolve  out  of  model  theory  (see  Appendix  A.l),  viz., 

location  of  the  aileron  hinge  line,  c 

chordwise  location  of  the  aileron  center  of  grav> 

ity, 

dimensionless  radius  of  gyration  of  the  aileron 
about  its  hinge  line,  r^ 

CA>3 

frequency  ratio,  — — 

OJ^ 
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Again,  the  problem  has  been  formulated  in  terms  ,o£  mode  shapes 
and  frequencies,  and  for  the  control  surface,  the  stiffness 
chordwise  is  considered  large  with  respect  to  the  stiffness 
spanwise.  The  values  of  c,  x^  and  r^  were  made  constant  at 
each  spanwise  station  in  the  model-design  procedure. 


For  wings  with  tip  tanks  (see  Appendix  A. 2)  :Ln  the  can¬ 
tilever  or  free-to-roll  condition,  we  must  define  the  geometry 
and  mass  of  the  tip  tank  in  addition  to  the  bare  wing  parameters 
already  considered: 


volume , 
location 


JT 

4b^I 

n,j, 

’  b~ 


M 

mass,  _T 

mX 


static  unbalance,  ST=  ST 

mb  jt 

X 

moment  of  inertia  in  pitch,  T 

mb^X 


moment  of  inertia  in  roll  of  the  root  support, 


These  parameters  have  been  non-dimens ionalized  in  a  somewhat 
arbitrary  fashion  because  of  the  simple  theoretical  model  con¬ 
sidered. 


Having  formulated  the  problem  in  terms  of  explicit  para¬ 
meters  for  all  configurations,  it  now  remains  for  us  to  devise  a 
scheme  of  model  design  whereby  desired  values  of  the  parameters 
can  be  obtained.  Also,  techniques  of  testing  the  models  at  zero 
airspeed  and  in  the  wind  tunnel  must  be  developed  to  verify  our 
design  values  and  to  obtain  the  desired  answers.  Detailed  con- 
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sideration  of  model  design  and  testing  techniques  are  presented 
in  Reference  14.  A  cursory  glance  at  the  highlights,  sufficient 
for  this  report,  is  given  in  Appendix  D. 

Since  the  emphasis  of  the  research  program  was  on  exper¬ 
imental  results,  the  theory  was  used  primarily  as  a  guide  for 
model  design.  Only  the  straight-wing  planform  was  treated  theo¬ 
retically  since  it  is  most  amenable  to  analysis.  The  calcula¬ 
tions  were  based  on  three-dimensional  structural  properties  and 
two-dimensional  supersonic  aerodynamic  forces  since  experience 
has  shown  that  a  similar  procedure  gives  reasonable  results  for 
subsonic  flutter  below  the  transonic  range.  -Also,  other  more 
complicated  methods  of  analysis  which  are  presently  being  de¬ 
veloped,  (e.g.  Reference  20)  were  considered  far  too  premature 
to  be  used  on  the  present  program.  Some  calculations  based  on 
Piston  Theory  were  made  because  of  its  simplicity.  In  formulating 
the  theory,  only  those  degrees  of  freedom  which  experience  had 
shown  to  be  essential  in  the  analysis  were  included.  Thus,  for 
the  bare-wing  and  cantilever  tip-tank  calculations,  only  first 
bending  and  first  torsion  were  included.  For  wings  with  control 
surfaces,  the  aileron  degree  of  freedom  was  added,  and  for  wings 
with  freedom  to  roll,  the  rigid-body  roll  motion  was  added. 

The  choice  of  eigenvalues  was  such  that  most  useful  in¬ 
formation  for  a  given  effort  could  be  obtained  from  a  model- 
design  viewpoint.  For  example,  in  the  theoretical  analysis  of 
the  bare  wing,  the  problem  was  set  up  so  that  the  bending  and 
torsional  frequencies,  which  are  measures  of  stiffness  level  were 
obtained  as  results  for  each  set  of  assumed  conditions  (fre¬ 
quency  and  Mach  number  at  flutter).  In  the  analysis  of  wings 
with  ailerons  and  wings  with  tip  tanks,  the  torsional  frequency 
along  with  the  aileron  frequency  for  the  former  and  the  tip-tank 
static  unbalance  for  the  latter  were  chosen  as  unknowns.  This 
choice  of  eigenvalues  still  allows  for  a  comparison  between 
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theory  and  experiment  so  long  as  the  complete  range  of  experi¬ 
mental  parameters  is  covered. 

The  theoretical  results  are  presented  in  terms  of  the 
flutter  coefficient,  >  rather  than  in  terms  of  a  velocity 

ratio,  as  has  been  done  in  the  past.  This  procedure  eliminates 
the  necessity  of  treating  velocity  (v) ,  size  (b)  and  stiffness 
level  (  uV)  as  separate  parameters.  Also,  this  dimensionless 
quantity  evolves  out  of  the  theory  (see  Appendix  C.l). 
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SECTION  III 

PRESENTATION  OF  RESULTS 


3 . 1  Bare  Wings 

3.1/1  Straight 


In  designing  the  bare  straight-wing  models,  an  attempt 
was  made  to  vary  the  first  torsional  frequency  of  the  wing  keep¬ 
ing  all  other  pertinent  dimensionless  parameters  constant,  viz., 


mass  ratio,  yu. 

location  of  chordwise  center  of  gravity, 

dimensionless  mass  moment  of  inertia  in  pitch, 

frequency  ratio,  -~ 

location  of  elastic  axis,  a 

spanwise  mass  and  stiffness  distributions. 


In  this  way,  experimental  curves  of  the  flutter  coefficient, 

V  f. 

,  and  the  reduced  frequency,  k,  could  be  determined  as  a 
function  of  Mach  number.  These  experimental  curves  could  then 
be  compared  with  those  obtained  by  theory  in  order  to  test  the 
validity  of  the  theory  and,  if  necessary,  aid  in  the  development 
of  a  criterion  for  torsional  rigidity. 


This  choice  of  parametric  variation  was  particularly 
difficult  to  accomplish  insofar  as  model  design  was  concerned 
because  of  the  interdependence  of  the  parameters  and  the  re¬ 
quirement  that  the  frequency  ratio,  ,  be  held  constant,  so 

that  some  variations  occurred  in  the  parameters  which  were  to 
be  held  constant.  Control  of  the  mass  ratio,  yu.  ,  was  difficult 
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because  of  the  change  in  air  density  with  Mach  number  and  the 
fact  that  the  Mach  number  at  flutter  could  not  be  accurately 
predicted  beforehand.  The  values  of  other  mass  parameters, 
center-of-gravity  location,  c.g.,  and  moment  of  inertia  in  pitch, 
r^  ,  as  well  as  the  spanwise  distributions  of  mass  and  stiff¬ 
ness  could  be  controlled  accurately.  Insofar  as  elastic  axis  is 
concerned,  there  was  some  question  as  to  how  this  concept  should 
be  handled  for' real  wings.  After  some  research  late  in  the  pro¬ 
gram,  it  was  decided  that,  based  on  classical  flutter  theory, 
the  elastic  axis  should  be  treated  as  the  locus  of  shear  centers 
rather  than  as  that  point  on  the  wing  where  bending  and  torsion 
is  statically  uncoupled,  "apparent"  elastic  axis  (see  Reference 
21).  Experiments  showed  that  the  locus  of  this  latter  point 
varied  slightly  along  the  wing  because  of  sweep  and  root  effects. 
Furthermore,  accurate  control  of  the  locus  of  shear  centers  could 
be  maintained  by  careful  construction  of  the  models,  since  the 
location  of  the  shear  center  is  a  function  of  the  position  of 
the  structural  elements  of  the  model.  Both  measured  "apparent 
elastic  axis"  and  calculated  locus  of  shear  centers  are  tabu¬ 
lated  in  Appendix  D  where  available.  Good  control  of  the  fre¬ 
quency  ratio,  ,  also  evolved  with  experience.  Careful 

selection  of  balsa  wood  and  realistic  estimates  of  the  effect  of 
glue  were  made  late  in  the  program.  A  detailed  tabulation  of 
the  parameters  for  all  the  models  tested  is  given  in  Appendix  D. 

In  order  to  determine  the  effects  on  flutter  of  those 
parameters  which  were  difficult  to  control  accurately  in  the 
models,  theoretical  straight  wing  studies  were  made  varying 
these  parameters.  Figure  3-1  presents  a  systematic  variation  of, 

Cx>. 

frequency  ratio,  - 
mass  ratio , 
elastic  axis,  a 
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FIGURE  3.1(a)  THEORETICAL  RESULTS  FOR  BARE  STRAIGHT-WING  PLAN- 

FORM,  e.a.  AT  39.3#  CHORD,  (^)2  =  0.09 
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FIGURE  3.1(b)  THEORETICAL  RESULTS  FOR  BARE  STRAIGHT-WING  PLAN- 

FORM,  e. a.  AT  39.3#  CHORD,  (~)2  =  0.16 
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FIGURE  3.1(c)  theoretical  results  for  bare  straight -wing  plan- 

form,  e.a.  AT  39.3^  CHORD,  (^)  2  =0.25 
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FIGURE  3.1(d) 


THEORETICAL  RESULTS  FOR  BARE  STRAIGHT -WING  PLAN 
FORM,  e.a.  AT  44.3#  CHORD,  (-^) 2  =  0.09 
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MACH  NUMBER 

FIGURE  3.1(e)  THEORETICAL  RESULTS  FOR  BARE  STRAIGHT -WING  PLAN- 

FORM,  e.a.  AT  44.3#  CHORD,  (—j')  2  =  0.16 
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FIGURE  3.1(f)  THEORETICAL  RESULTS  FOR  BARE  STRAIGHT-WING  PLAN- 

FORM,  e.a.  AT  44.3#  CHORD 2  =0.25 
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for  fixed  spanwise  mass  and  stiffness  distributions,  fixed 
moment  of  inertia  in  pitch  and  fixed  chordwise  center-of-gravity 
location.  Taper  was  taken  into  account  in  determining  the  aero¬ 
dynamic  forces.  All  the  curves  of  Figure  3.1  were  based  on  zero 
structural  damping.  However,  spot  checks  showed  that  inclusion 
of  a  small  amount  of  structural  damping  (g  =  0.01)  in  the  analy¬ 
sis  did  not  alter  the  answers  appreciably,  probably  because  all 
of  the  calculated  cases  had  reasonable  amounts  of  positive  static 
unbalance.  Two-dimensional  aerodynamic  coefficients  and  three- 
dimensional  structural  properties  were  used  in  the  theoretical 
calculations.  The  detailed  formulation  of  the  theory  is  given 
in  Appendix  A.l,  and  an  example  analysis  is  presented  in  Appen¬ 
dix  C.2.  Theoretical  results  for  values  of  frequency  ratio, 

,  other  than  those  presented  in  Figure  3.1  can  be  obtain¬ 
ed* 

ed  by  cross-plotting  the  curves  of  Figure  C.l.  Because  avail¬ 
able  tabulated  values  of  the  supersonic  oscillatory  aerodynamic 
coefficients  were  limited  (Ref.  7),  portions  of  the  curves  of 
Figure  3.1  could  not  be  defined  accurately.  These  doubtful 
portions  are  presented  as  dashed  lines. 

Figure  3.1  shows  that,  in  the  range  of  practical  inter¬ 
est  (  ~~ —  ^5  ),  the  theoretical  stability  boundaries  are  not 

Oft  <“'<*. 

very  sensitive  to  changes  in  the  parameters,  mass  ratio  and 
elastic-axis  location  (for  the  ranges  considered)  but  do  vary 
somewhat  with  frequency  ratio,  J— .  We  also  notice  that  the 
trends  are  comparable  to  those  to  be  expected  from  experience  in 
subsonic  flow  (Ref.  22),  i.e.,  the  region  of  instability  in¬ 
creases  with, 

decreasing  mass  ratio,  yu 

increasing  distance  between  the  elastic  axis  and 
the  center  of  gravity, 

increasing  frequency  ratio, 
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The  theoretical  curves  of  Fig.  3-1  may  be  used  to  predict 
the  onset  of  flutter  which  will  be  given  by  the  intersection  of 
an  operating  line  and  the  flutter  boundary.  The  operating  line 
for  a  model  with  a  given  set  of  mass  and  stiffness  parameters  de¬ 
pends  on  the  velocity-  Mach  number  relationship  of  the  environment 
in  which  the  model  is  tested.  Velocity  at  a  given  Mach  number  is 
a  function  of  the  ambient  temperature  only.  Velocity  versus  Mach 
number  for  extreme  ranges  of  the  atmosphere  and  of  the  facility 
used  for  flutter  testing  on  this  program  is  given  in  Figure  3-2. 

In  the  tunnel  the  stagnation  temperature  is  roughly  that  of  the 
atmosphere,  and  at  a  given  Mach  number  the  static  temperature,  and 
hence  the  speed  of  sound  and  velocity,  is  less  than  atmospheric  in 
accordance  with  isentropic  flow  relations.  Figure  3*2  shows 


graphically  the  difference  between  the  velocity-  Mach  number  re¬ 
lationship  in  the  atmosphere  and  the  wind  tunnel.  A  given  model 
rfill  then  have  a  different  operating  line  in  the  tunnel  than  it 


ATMOSPHERIC  AND  WIND  TUNNEL  CONDITIONS 
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has  in  the  atmosphere.  Figure  3-3 s  illustrates  this  point,  it 
shows  the  operating  lines  for  a  model  in  the  wind  tunnel  and  in  the 
standard  atmosphere  at  sea  level  and  at  the  tropopause.  Since 
these  operating  lines  are  not  the  same  they  intersect  the  flutter 
boundary  at  different  values  of  Mach  number.  It  is  also  interest¬ 
ing  to  note  that  if  the  flutter  boundary  is  to  be  approached  from 
a  stable  region  at  either  constant  altitude  in  the  atmosphere  or  in 
the  tunnel,  it  must  be  approached  by  decreasing  the  Mach  number. 

For  all  of  the  straight  wing  models  tested  on  the  present 
program  the  theoretical  Mach  number  of  flutter  has  been  determined 
from  curves  similar  to  Figure  3.3.  The  flutter  Mach  number  has 
been  determined  for  conditions  in  the  tunnel  and  in  the  standard 
atmosphere.  The  theoretical  Mach  numbers  of  flutter  so  determined 
are  compared  with  the  flutter  Mach  number  of  the  actual  tests. 

These  results  are  given  in  Table  3.1. 

An  even  more  graphic  comparison  between  theory  and  exper¬ 
iment  is  given  in  Figure  3-4.  All  the  legitimate  experimental 
flutter  points  were  obtained  by  injecting  the  model  into  a  stable 
region  and  decreasing  the  Mach  number  until  flutter  occurred.  All 
the  straight  wings  built  on  this  program  fluttered,  and  no  single¬ 
degree -of -freedom  torsion  flutter  was  encountered  experimentally 
(see  Ref.  23,  p.  6).  The  scatter  of  the  data  is  small  and  most  of 
the  deviations  from  the  mean  can  be  explained.  The  models  which 
fluttered  during  injection  would  have  their  marginally  stable  con¬ 
dition  at  higher  Mach  number.  These  injection  flutters  generally 
occurred  at  frequencies  closer  to  the  first  torsional  frequency 
than  would  otherwise  be  expected.  Models  ST- Id  and  ST-ld-1  were 
designed  with  elastic  axes  forward  so  that  they  had  slightly  high¬ 
er  values  of  the  coefficient,  ,  than  the  other  models,  veri¬ 

fying  the  trend  predicted  by  the  theory.  Model  ST-12  was  a  low- 
density  =30,  wing  designed  to  have  a  margin  of  safety  against 
flutter  based  on  the  test  results  of  models  with  mass  ratio  on  the 
order  of  65  and  on. the  theoretical  effect  of  lowering  the  mass 
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FIGURE  3-3  EXAMPLE  USE  OF  THEORETICAL  CURVES 
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TABLE  3.1 


MACH  NUMBER  AND  FREQUENCY  AT  FLUTTER  FOR  STRAIGHT-WING  PLANFORM,  THEORY  AND  EXPERIMENT 


Experiment 


Wind  Tunnel 


ST-ld-1  1.52 
ST-2*  1.71 
ST-4  1.52 
ST-4-1  1.30 
ST-5  1 .44 
ST-5-1  1.47 
ST-6  1.72 
ST-7-1*  1.88 
ST-7-2*  1.92 
ST-7-3  1.91* 
CT-12  1.45 


ST- la*  1.80 


ST-lc*  1.33 


ST -4a  1.43 


ST-4b*  1.83 


Wind  Tunnel 


Theor 


Tropopause 


Bare  Wrings 


Sea  Level 


.80  200.0 

.72  |  87.9 

8".  5 


Wings  with  Tip  Tanks  (Cantilever) 


unstable 

unstable 

1.56 

1.43 

1.40 

stable 

unstable 

unstable 


unstable 

unstable 


1.48 

stable 

unstable 

unstable 


Wings  with  Tip  Tanks  ( free-to-roll ) 


c-1 

c-2 

Kg® 

33.9 

30.0 

1.57 

1.57 

50 

49 

c-3 

1.32 

32.7 

1.57 

50 

Wings  with  Ailerons 


Structural- 
Damping 
Coefficient 
in  Theory 


52 

52 

1.62 

1.62 

59 

59 

!  “ — 

0 

52 

1.62 

59 

0 

MIIIOI 

i  i  m  my&J  [ssisdi 


‘Injection  flutter 

Note1  tip-tank  solutions  attempted  only  inside  interval  1.3  <  M  <1.8 


WADC  TR  54-113,  Part  II 


CONFIDENTIAL 


ONFIDE 


f-2b„H 

JILL 

23  — 1  b.  f— 

THEORY 

EXPERIMENT 

mmm&) 

)S.  ■  65 

/i?65 

e  g.  <a  50% 'CHORD 

eg  @  50%  CHORD 

rJ  -0.230 

.42C 

ri  -0.250 

°.42C 

«.a.@  42.3  %  CHORD 

e.a. »  42%  CHORD 

Sf1011 

O 

g  -0 

'  a  - 

g  s  0.015 

1  _  ..... 

*  INJECTION 

FLUTTER 

MS— 


MACH  NUMBER 

FIGURE  3.4  THEORETICAL  AND  EXPERIMENTAL  FLUTTER  STABILITY 
BOUNDARIES  FOR  BARE  STRAIGHT-WING  PLANFORM 
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ratio.  When  the  model  fluttered,  the  experimental  results  were 
re-interpreted  on  the  basis  of  the  parameter,  •  Theo¬ 

retically,  the  flutter  coefficient  and  the  mass  ratio  are  inde¬ 
pendent  parameters  and  consequently  cannot  be  combined  into  a 
single  quantity.  Piston  theory  shows,  however,  that  for  the  case 
of  high  supersonic  bending-torsion  flutter  of  a  representative 
two-dimensional  section,  this  combination  can  be  demonstrated 
analytically  for  some  special  cases  (Ref/  24,  p.  81  ff ) .  The 
semi-empirical  flutter  formula  proposed  in  Reference  22  (p.  17) 
for  subsonic  flow  also  suggests  this  combination. 


Figure  3-5  shows  the  results  of  plotting  the  experimen¬ 
tal  data  on  the  basis  of  this  parameter.  We  see  that  the  scatter 
of  Model  ST-12  is  considerably  reduced  with  no  appreciable  loss 
of  consistency  in  the  other  data  points.  The  range  of  mass  ratio 
is  not  large  enough  to  determine  if  this  trend  is  general  however. 


Figure  3-4  emphatically  shews  the  value  of  being  able 
to  vary  the  Mach  number  during  a  flutter  model  test.  If  a 
fixed  nozzle  of  low  supersonic  Mach  number  were  used,  we  might 
conclude  that  the  theory  gave  good  agreement  with  experiment, 
and  there  would  be  some  doubt  as  to  the  location  of  the  neutrally 
stable  region.  The  gross  deviation  between  theory  and  experiment 
at  higher  Mach  numbers  led  to  a  theoretical  investigation  using 
Piston  theory.  The  theoretical  formulation  is  presented  in 
Appendix  A. 3,  and  a  numerical  example  is  given  in  Appendix  C.6. 
The  results  of  this  investigation  are  shown  in  Figure  3.6.  It 
should  be  emphasized  that  this  theory  is  valid  only  at  high  Mach 
numbers,  well  out  of  the  range  of  the  experimentation  of  this 
program.  We  see  that  if  thickness  is  taken  into  account,  the 
results  of  Piston  Theory  give  better  agreement  with  experiment 
than  the  results  of  the  classical  theory  near  Mach  2.  Piston 
Theory  certainly  does  not  apply  at  Mach  numbers  below  2.0., 
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FIGURE  3.6  THEORETICAL  RESULTS  FOR  BARE  STRAIGHT-WING  PLANFORM 

BASED  ON  PISTON  THEORY 
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FIGURE  3.6  (Continued)  THEORETICAL  RESULTS  FOR  BARE  STRAIGHT-WING 

P'LANFORM  BASED  ON  PISTON  THEORY 
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Reference  25,  based  on  potential  flow  solutions,  shows  that  the 
effect  of  thickness  on  the  damping  of  an  oscillating  airfoil  in 
the  range  of  Mach  number  from  2.0  to  about  1.2  may  be  either 
,'tabilizing  or  destabilizing  depending  on  the  Mach  number,  re¬ 
duced  frequency,  and  other  parameters.  Reference  33,  based  on 
Piston  Theory,  shows  that  the  effect  of  thickness.,  for  the  cases 
considered  there,  is  destabilizing  and  that  the  effect  increases 
with  increasing  Mach  number.  It  should  also  be  mentioned  that 
even  though  the  Piston  Theory  shows  a  tremendous  effect  due  to 
thickness,  this  effect  may  not  be  as  large  as  the  results  in¬ 
dicate,  because  terms  have  been  neglected  in  the  analysis  which 
could  be  of  the  same  order  of  magnitude  as  the  thickness  terms. 

The  results  are  presented  in  order  to  determine  trends  for  a 
better  insight  into  what  happens  at  high  Mach  number. 

For  purposes  of  comparison,  Figure  3.7  gives  the  re¬ 
sults  of  a  theoretical  flutter  analysis  based  on  incompressible 
aerodynamic  coefficients  (Ref.  26).  A  numerical,  example  is 
given  in  Appendix  C.l.  It  is  interesting  to  note  that  the  in° 
compressible  theory  gives  almost  the  same  results  as  the  super¬ 
sonic  theory  near  Mach  number  1.35. 

It  is  interesting  to  compare  the  flutter  parameters  of 
the  models  tests  on  this  program  with  those  of  some  actual  high¬ 
speed  airplanes.  Table  3-2  gives  flutter  parameters  for  three 
representative  airplanes.  The  data  in  Table  3-2  has  been  derived 
from  actual  airplane  data,  averaged  to  give  representative 
straight,  swept  and  delta  wing  airplanes.  The  torsional  frequency 
of  the  scaled  model,  is  included  in  Table  3-2  so  that  a 

direct  comparison  could  be  made  of  the  model  and  actual  airplane 
stiffness.  From  dimensional  analysis  it  can  be  shown  that  if  Mach 
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FIGURE  3.7  THEORETICAL  RESULTS  FOR  BARE  STRAIGHT-WING  PLANFORM 

BASED  ON  INCOMPRESSIBLE  THEORY 
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number  is  to  be  kept  constant  from  model  to  airplane,  the  follow¬ 
ing  relationship  holds, 

=  k&j  Eq.  (3.1) 

where 

Id#  is  the  first  torsional  frequency  of  the  wing 
L  is  length 
T  is  temperature 

The  subscripts,  M  and  A,  refer  to  model  and  airplane,  respective¬ 
ly.  The  scaled  model  frequencies,  ( of  Table  3-2  were  cal¬ 
culated  assuming  a  scaled  model  of  the  airplane  with  a  root  chord 
equal  to  that  of  the  wings  tested  on  this  program  (10  inches). 

Air  was  assumed  as  the  model  testing  medium.  The  air  was  assumed 
to  expand  adiabatically  from  a  stagnation  temperature  equal  to 
room  temperature  (70°F)  (See  Reference  12). 


3.1.2  Swept 

The  experimental  data  for  the  bare  swept  wings  are  pre¬ 
sented  graphically  in  Figure  3-8;  no  theoretical  work  was  done 
for  this  planform.  Less  extensive  tests  were  conducted  on  the 
swept  than  on  the  straight  wings,  and  some  of  the  models  tested 
did  not  flutter  (see  Table  D.8).  Again,  all  legitimate  flutter 
points  were  obtained  by  injecting  the  model  into  a  stable  region 
and  decreasing  the  Mach  number  until  flutter  occurred.  The  re¬ 
marks  made  above  pertaining  to  control  of  the  parameters  for  the 
straight  wings  are  also  applicable  to  this  planform. 
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FIGURE  3.8  EXPERIMENTAL  FLUTTER  STABILITY  BOUNDARY  FOR  BARE 

SWEPT-WING  PLANFORM 
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In  Figure  3*8  those  models  which  fluttered  during  in¬ 
jection  would  probably  be  neutrally  stable  at  higher  values  of 
Mach  number,  and  the  model  which  fluttered  during  retraction 
would  probably  be  neutrally  stable  at  a  lower  Mach  number. 

Models  SW-3c,  SW-3  and  SW-3d  were  purposely  designed  to  have 
three  different  elastic-axis  locations  (loci  of  shear  centers). 
All  other  parameters  were  held  very  nearly  constant  (see  Tables 
D.5  &  D.8) .  The  theoretical  and  experimental  trend  evidenced  by 


Model 

e.a. 

%  chord 

Mf 

to^ 

cps 

SW-3C 

37-0 

1.48 

135 

SW-3 

42.0 

1.41 

132 

SW-3d 

47.0 

1.25 

137 

the  straight-wing  planform  holds,  i.e.,  moving  the  elastic  axis 
forward  from  a  given  center-of-gravity  location  is  equivalent 
to  a  decrease  in  stiffness.  Therefore,  on  Figure  3-8,  Model 
SW-3c  should  have  a  slightly  higher  and  Model  SW-3d  a  slightly 
lower  value  of  the  flutter  coefficient,  r  — —  ,  than  the  mean. 


bD 


The  destabilizing  effect  of  lowering  the  mass  ratio, 
is  evidenced  by  Mod^l  SW-8  (See  Table  D.8)  which  was  tested  once 
without  flutter.  In  a  second  test  during  which  the  mass  ratio 
was  lowered  by  raising  the  air  density  in  the  test  section,  the 
model  fluttered  during  injection.  If  the  experimental  results 
are  plotted  on  the  basis  of  the  parameter,  —  ,  (Figure 

3.9),  there  is  no  appreciable  increase  in  scatter,  but  again 
there  is  no  sound  theoretical  basis  for  this  choice  of  parameter, 
and  it  should  be  used  cautiously  because  of  the  scantiness  of 
the  data. 
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The  experimental  stability  boundary  for  the  swept  plan- 
form  occurs  at  about  the  same  values  of  the  flutter  coefficients, 
and  Jj=  ,  as  for  the  straight  planform  at  low 

supersonic  Mach  number  while  the  slope  of  the  boundary  with  in- 
creasing  Mach  number  does  not  appear  to  be  as  steep  as  that  for 
the  straight  wing.  The  same  effect  of  sweepback  has  also  been 
obtained  from  flutter  tests  at  transonic  speeds  (see  Reference  4, 
p.  27  and  Reference  6).  The  numerical  equivalence  of  flutter 
coefficients  probably  means  that  the  swept  wing  is  slightly  less 
desirable  from  a  flutter  viewpoint  than  the  straight  wing. 

While  it  is  difficult  to  compare  wings  of  different  geometry, 
some  remarks  can  be  made.  From  elementary  rod  considerations, 
the  following  relation  can  be  deduced  for  a  wing  of  given  root 
stiffness  and  of  given  root  chord, 


where 


•£(/*)  I  cos  4- 


Eq.  (3.1) 


f/'  is  the  sweep  of  the  elastic  axis 

/  is  the  semi- span  of  the  wing,  measured  from  the  root 
chord. 


The  effect  of  taper,  A  >  on  the  torsional  frequency  has  been 
estimated  from  the  relation, 

tapered  ~  untapered  ^ *  ..P-2) 

where  the  function,  fj  (A),  is  obtained  from  Reference  27  as 

1.6 

fx  (A)  =  1  +  1.87  (1  -  A  )  Eq.  (3.3) 
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It  should  be  mentioned  that  the  effect  of  taper  on  the  bending 
frequency  is  more  difficult  to  estimate,  and  the  methods  of 
Reference  28  are  suggested.  On  the  basis  of  equation  3.2,  for 
the  experimental  parameters  of  Figures  3.4  and  3*8,  the  straight 
wing  should  have  a  higher  torsional  frequency  by  about  10  percent 
than  the  swept  wing  for  the  same  root  stiffness  and  the  same  root 
chord. 


3.1.3  Delta 

Considerable  difficulty  was  encountered  precipitating 
delta-wing  flutter.  Some  twenty  models  were  flown  before  suc¬ 
cessful  flutter  was  obtained,  utilizing  the  extreme  values  of 
center-of-gravity  and  elastic-axis  locations  available  with  the 
model-design  procedure.  Even  models  made  only  of  balsa  and  lead 
were  flutter-free.  Lead  was  added  to  some  of  the  wings  which 
had  been  tested  without  event  in  an  attempt  to  lower  the  second 
bending  frequency  without  significantly  changing  first  bending 
or  first  torsion,  since  there  was  some  evidence  that  a  critical 
condition  might  be  one  where  the  second  bending  frequency  is  be¬ 
low  first  torsion  (Ref.  29). 

Successful  flutter  was  finally  obtained  by  reducing  the 
thickness  ratio  of  the  wings  to  4  percent,  thus  reducing  the 
stiffness  well  below  that  which  would  be  encountered  in  practice. 
To  illustrate  this  point  one  need  only  compare  the  scaled  tor¬ 
sional  frequency  of  the  delta-wing  airplane  presented  in  Table 
3.2  with  those  of  the  delta-wing  models  tested  on  this  program 
(Table  D.9).  The  frequency  of  the  model  of  the  actual  aircraft, 
where  Mach  number  is  kept  constant  from  model  to  airplane,  is 
an  order  of  magnitude  higher  than  those  of  any  model  built  on 
this  program.  A  similar  conclusion  can  be  drawn  by  comparing 
the  flexibility  influence  coefficients  of  the  models  tested. 
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(Table  D,18)  with  those  of  the  XF-92A  airplane  (Table  D.19) 
using  the  relation, 

(C1J)M  /?«  T„  Eq> 


where 


is  the  force-deflection  influence  coefficient 
is  the  air  density. 


Equation  3-4,  which  is  based  on  dimensional-analysis  considera¬ 
tions,  assumes  that  the  mass  ratio  as  well  as  the  Mach  number  is 
held  constant  from  model  to  airplane.  It  is  interesting  to  note 
that  the  ratio  of  influence  coefficients  depends  only  on  flow 
parameters  and  is  independent  of  the  length  scale  factor.  This 
ratio  is  equal  to  1/3  if  the  XF-92A  airplane  is  simulated  at  an 
altitude  of  35,000  in  a  wind  tunnel  where  air  is  expanded 
isentropically  from  a  stagnation  temperature  equal  to  room  tem¬ 
perature  (70°)  and  a  stagnation  pressure  equal  to  10  psig. 


Figure  3.10  shows  the  flutter  stability  boundaries  ob¬ 
tained  experimentally.  Again,  each  successful  flutter  point, 
except  those  otherwise  indicated,  represents  a  model  which  was 
injected  into  the  airstream,  flutter-free.  Flutter  was  then 
approached  by  decreasing  the  Mach  number  until  instability 
occurred.  Figure  3.11  presents  the  results  of  plotting  the  ex- 
perimental  data  on  the  basis  of  the  coefficient,  -r —  J  — 

r  b.  LOrt.  v*  f- 

Experimental  evidence  of  the  destabilizing  effect  of  lowering 
the  mass  ratio,  yu.  ,  was  obtained  in  the  case  of  Model  De-2d-2 
which  was  tested  once,  flutter-free,  at  higher  mass  ratio  than 
that  obtained  at  flutter  during  a  second  test  (see  Table  D.9). 
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The  curves  of  Figures  3-10  and  3*H  are  not  as  well  defined  as 
those  of  the  straight  and  swept-wing  planforms  because  of  the 
inherent  difficulty,  with  a  simple  design  procedure,  of  accurate¬ 
ly  controlling  the  natural  frequencies  of  delta-wing  models  where 
the  effects  of  plate- type  vibratory  modes  become  significant. 
While  there  is  little  justification  for  quantitatively  comparing 
the  stability  boundary  of  the  delta  wing  with  those  of  the 
straight  and  swept  wings,  even  on  the  basis  of  equation  3.1>  the 
delta  wing  should  have  its  stability  boundary  at  much  lower 
values  of  the  coefficient,  >  t^an  obtained  experimen¬ 

tally  (about  1/4  that  of  the  straight  wing)  for  the  same  root 
chord  and  the  same  root  torsional  rigidity.  It  i-s  interesting 
to  note  that  the  same  general  changes  in  stability  boundary  with 
sweepback  as  those  obtained  on  this  program  have  also  been  ob¬ 
tained  from  transonic  tests  (see  Reference  4,  p.  27  and  Refer¬ 
ence  6) . 

3 . 2  Wings  with  Ailerons 
3.2.1  Straight 

In  the  theoretical  formulation  of  the  straight  wing 

00/2 

with  aileron,  the  frequency  ratio,  —  ,  was  considered  as  an  un¬ 
known  in  the  analysis  with  the  fixed  parameters  shown  on  Figure 
3.12.  The  detailed  formulation  of  the  theory  is  given  in 
Appendix  A. 2,  and  a  numerical  example  is  given  in  Appendix  C.3. 
Only  three  models  were  tested  so  that  a  direct  graphical  com¬ 
parison  between  theory  and  experiment  could  not  be  made  as  in 
the  case  of  the  bare  straight  wing.  A  tabulated  comparison  of 
the  Mach  number  and  the  frequency  at  flutter  is  given  in  Table 
3.1,  and  a  detailed  tabulation  of  the  experimental  parameters 
is  given  in  Appendix  D. 
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MACH  NUMBER 

FIGURE  3.12  THEORETICAL  RESULTS  FOR  STRAIGHT  WING  WITH  AILERON 
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Again,  t he  theory  shows  that  for  flutter  at  super¬ 
sonic  Mach  numbers,  instability  must  be  approached  from  high 
Mach  number.  decreasing  the  frequency  ratio,  ~  ,  is  an  effec- 
tive  decrease  in  torsional  stiffness,  and  no  significant  change 
occurs  near  the  condition,  =  1.0  ,  as  would  be  expected  from 
subsonic  experience. 

These  trends  have  been  verified  by  experiment  for  the 
three  models  tested. 


Model 

COtf 

aJ>c< 

Mf 

ST-lb 

0.691 

1.80 

ST-le 

1.49 

1.72 

ST-lf 

00 

1.65 

*  Injection  Flutter 

Model  ST-lf  was  the  result  of  the  repairing  Model  ST-le  and 
locking  the  aileron.  Increasing  the  frequency  ratio,  , 
lowers  the  Mach  number  at  flutter  or  has  the  same  effect  as  in¬ 
creasing  the  absolute  stiffness,  .  Although  the  trends  are 
correct,  the  theory  is  again  unconservative  in  that  experiment 
shows  a  larger  region  of  instability  than  the  theory  (see  Table 
3.1).  The  experimental  flutter  frequencies  are  also  lower  than 
those  predicted  by  the  theory  except  in  the  case  of  Model  ST-lb 
which  fluttered  during  injection  near  the  first  torsional  fre¬ 
quency. 


In  the  design  of  the  aileron  models,  the  aileron  fre¬ 
quency  was  controlled  by  flexures  even  though  such  a  design  was 
difficult,  especially  since  all  the  models  had  sealed  gaps, 
typical  of  high-speed  aircraft.  If  the  aileron  frequency  were 
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controlled  from  outside  the  wind  tunnel,  a  real  wing  would  not 
be  simulated  because  the  aileron  would  rotate  relative  to  a  line 
fixed  in  space  rather  than  relative  to  the  wing  itself. 

3.2.2  Swept  and  Delta 

The  trends  obtained  for  the  straight  wing  with  aileron 
were  also  evidenced  for  the  swept  and  delta  planforms,  although 
no  theoretical  work  was  done  and  only  one  model  with  aileron  was 
tested  for  each  of  these  planforms.  The  data  show  that  Model 
SW-3b,  which  had  roughly  the  same  parameters  as  Models  SW-3  and 
SW-8-1  except  for  finite  aileron  frequency  (see  Tables  D.5  and 
D.8) ,  fluttered  at  a  higher  Mach  number  than  either  of  these 
latter  two.  The  delta-wing  elevon  model,  DE-2e,  fluttered  during 
injection  in  a  mode  in  which  aileron  motion  predominated  at  much 
higher  Mach  number  than  the  experimental  results  of  the  bare 
delta  wings  would  indicate. 

3.3  Wings  with  Tip  Tanks 

In  addition  to  tests  on  .bare  wings  and  on  wings  with 
ailerons,  some  exploratory  tests  were  made  on  straight  and  swept 
wings  with  tip  tanks  for  both  cantilever  and  free-to-roll  root 
conditions.  Some  theoretical  work  was  done  on  the  straight  wing 
for  a  simplified  model  (see  Figure  3.13)  so  that  trends  useful 
for  model  design  could  be  obtained,  In  the  theoretical  formula¬ 
tion  of  the  problem,  the  dimenionless  static  unbalance  of  the 
tip  tank,  "3^ ,  and  the  first  torsional  frequency  of  the  bare  wing, 
co^  ,  were  treated  as  unknowns.  The  frequencies  of  the  bare 
wing  were  chosen  as  independent  parameters  because  those  of  the 
wing  with  tip  tank  depend  on  the  tip-tank  static  unbalance, 
which  was  to  be  varied  both  theoretically  and  experimentally. 
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3.3.1  Cantilever  Root  Condition 

The  results  of  the  theoretical  calculations  on  the 
simplified  model  of  the  straight -wing  planform  for  zero  struc¬ 
tural  damping  are  presented  in  Figure  3-13.  The  detailed  formu¬ 
lation  of  the  theory  is  given  in  Appendix  A. 2,  and  a  numerical 
example  is  given  in  Appendix  C.4.  It  should  be  mentioned  that 
the  aerodynamic,  forces  on  the  tip  tank  were  taken  into  account 
in  the  calculations  (see  Appendix  B)  and  were  significant, 
especially  insofar  as  the  pitching  moment  was  concerned.  Theo¬ 
retical  results  for  values  of  the  dimensionless  tip-tank  static 
unbalance,  1?T,  other  than  those  presented  in  Figure  3.13  can  be 
obtained  by  cross-plotting  the  curves  of  Figure  C.4. 

The  curves  of  Figure  3.13  are  at  best  confusing.  It 
can  be  shown  mathematically  that  the  curves  must  be  continuous, 
so  that  the  theory  seems  to  indicate  that  the  wing  with  tip 
tank  is  unstable  even  for  infinite  stiffness  (  - , =0).  Con- 
siderable  difficulties  were  encountered  in  the  calculations  be¬ 
cause  of  the  limited  tabulated  values  of  the  supersonic  oscilla¬ 
tory  aerodynamic  coefficients  (Ref.  7)  and  the  choice  of  eigen¬ 
values. 


The  extreme  sensitivity  of  the  theoretical  curves  and 
the  disappointing  results  led  to  the  decision  to  do  some  of  the 
calculations  including  a  small  amount  of  structural  damping. 

The  results  of  these  calculations  are  presented  in  Figure  3.14. 
Theoretical  curves  for  values  of  dimensionless  tip- tank  static 
unbalance,  other  than  those  presented  can  be  obtained  by 
cross-plotting  the  curves  of  Figure  C.5.  The  curves  of  Figure 
3.14  are  more  encouraging,  following  the  trends  typified  by 
bare  wings  and  by  wings  with  ailerons. 
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FIGURE  3.13  (Continued)  THEORETICAL  RESULTS  FOR  STRAIGHT  WING 

WITH  TIP  TANK,  NO  STRUCTURAL  DAMPING 

WADC  TR  5^-113,  Part  II  49 


CONFIDENTIAL 


CONFIDENTIAL 


MACH  NUMBER 

FIGURE  3.14  THEORETICAL  RESULTS  FOR  CANTILEVER  STRAIGHT  WING  WITH 

TIP  TANK,  STRUCTURAL  DAMPING  INCLUDED 
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The  theory  shows  that  instability  increases  with  in¬ 
creasing  tip-tank  static  unbalance ,  Srp,  corresponding  to  rear¬ 
ward  travel  of  the  tip-tank  center  of  gravity.  Few  real  roots 
were  obtained  for  negative  tip-tank  static  unbalance,  i.e.  for 
tip-tank  center-of-gravity  locations  ahead  of  the  wing  elastic 
axis  (see  Figure  C.5),  indicating  that  the  wing  is  always  stable 
if  the  center  of  gravity  of  the  tip  tank  is  far  enough  forward. 
The  trends  predicted  by  the  theory  were  verified  by  experiment, 
but  again  the  theory  is  unconservative  if  quantitative  results 
are  compared.  Four  cantilever  tip- tank  models  were  tested,  from 
the  most  aft  to  the  most  forward*  tip-tank  center-of-gravity 
location  obtainable  with  the  model  design.  All  the  configura¬ 
tions  with  tip-tank  center  of  gravity  aft  of  the  wing  elastic 


Model 

Tip  Tank 
c.g. 

%  Wing 

Tip  Chord 

Mf 

ST-4b* 

I.05 

77.5 

1.83 

ST-la* 

0.401 

56.0 

1.80 

ST-lc* 

0.0506 

39-1 

1.83 

ST -4a 

0.5H 

25.0 

1.43 

*  Injection 

.  Flutter 

axis  fluttered  during  injection,  even  the  model  with  a  tip-tank 
static  unbalance  of  very  nearly  zero  (Model  ST-lc).  The  only 
successful  flutter  was  obtained  with  a  model  whose  tip-tank 
center  of  gravity  was  well  forward  (Model  ST-4a) .  According  to 
the  theory,  this  latter  model  should  have  been  stable  for  all 
Mach  numbers  above  1.3.  A  tabulated  comparison  of  theory  and 
experiment  is  presented  in  Table  3.1,  and  detailed  parameters 
for  the  models  tested  are  given  in  Appendix  D. 
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The  stabilizing  effect  of  moving  the  tip-tank  center 
of  gravity  forward  has  also  been  demonstrated  experimentally  at 
high  subsonic  speeds  (Ref.  30). 

Three  swept-wing  models  with  tip  tanks  were  tested;  no 
theoretical  studies  were  made  for  this  planform.  Although  the 
tip-tank  center-of-gravity  location  was  not  varied  as  above  for 
the  straight  wings,  the  models  did  have  different  values  of  the 
torsional  frequency  of  the  bare  wing,  with  essentially  constant 
values  of  other  parameters  (see  Table  D.ll).  The  trend  predict¬ 
ed  by  the  theory  was  verified,  i.e.  increasing  the  torsional 
frequency  of  the  wing  lowered  the  Mach  number  at  flutter. 


Model 

sT 

Tip -Tank 
c.g. 

%  Wing 

Tip  Chord 

<4* 

cps 

Mf 

SW-3a* 

-0.51 

25.0 

145 

1.92 

SW-7 

-0.51 

25.0 

161 

1.44 

SW-6 

-0.53 

25.0 

I89 

1.30 

* 


*  Injection  Flutter 

3.3.2  Free-to-Roll  Root  Condition 

The  results  of  the  theoretical  calculations  for  the 
simplified  straight-wing  model  with  free-to-roll  root  condition 
are  presented  in  Figure  3- 15-  The  theoretical  formulation  is 
given  in  Appendix  A. 2,  and  a  numerical  example  is  given  in 
Appendix  C.5.  The  calculations  were  made  for  zero  structural 
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damping  only. 

Again  the  theoretical  results  are  confusing,  but  it  did 
not  seem  worthwhile  to  expend  additional  effort  redoing  the  cal¬ 
culations  with  structural  damping  included  since  the  major  effort 
of  the  research  program  was  on  experimental  results.  The  solu¬ 
tion  of  the  three-degree-of-freedom  flutter  determinant  for  the 
free-to-roll  root  condition  involves  tedious  graphical  solutions 
of  simultaneous  equations  for  each  set  of  assumed  conditions,  i.e. 
the  Mach  number  and  the  frequency  at  flutter  (see  Appendix  C.5). 

Some  preliminary  calculations  indicated  that  with  the 
inclusion  of  structural  damping  in  the  analysis,  one  of  the 
branches  of  the  curves  would  disappear  (see  Figure  C.6)  and  a 
general  shift  similar  to  that  for  the  cantilever  root  condition 
would  occur.  The  condition  is  one  which  is  fairly  common  in 
flutter  analyses,  where  small  changes  in  input  give  erratic 
changes  in  results. 

In  the  theoretical  formulation,  no  account  was  taken  of 
the  canted  hinge,  which  was  used  to  give  the  models  aerodynamic 
stability  in  roll,  although  this  would  be  possible  in  a  more  re¬ 
fined  analysis.  The  amount  of  cant  was  kept  as  low  as  possible, 

2  1/2°,  consistent  with  the  change  in  flow  direction  with  Mach 
number,  about  +0.1°  (see  Reference  13).  For  small  angles  it 
can  be  shown  that  there  is  a  simple  relation  between  cant  angle, 
roll  angle,  and  angle  of  attack, 

0  ^  Eq.  (3.6) 
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FIGURE  3.15  THEORETICAL  RESULTS  FOR  STRAIGHT  WING  WITH  TIP  TANK, 


FREE  TO  ROLL 
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FIGURE  3.15  (Continued)  THEORETICAL  RESULTS  FOR  STRAIGHT  WING 

WITH  TIP  TANK,  FREE  TO  ROLL 
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where 

c*  is  the  angle  of  attack 

0  is  the  angle  of  roll 

£  is  the  angle  of  cant  in  radians. 

In  any  wind  tunnel  test  where  a  semi-span  model  is 
given  freedom  to  roll,  the  classical  question  of  whether  or  not 
the  anti-symmetrical  degrees-of-freedom  are  adequately  simulated 
always  arises.  This  uncertainty  comes  from  the  fact  that  the 
model  plus  its  mirror  image  are  actually  represented  (see  Figure 
3.16).  Because  of  this  doubt  along  with  the  canted  hinge,  the 
results  for  the  models  which  were  tested  probably  lie  somewhere 


MIRROR 

IMAGE 


FIGURE  3.16  SEMI-SPAN  MODEL  WITH  FREEDOM  TO  ROLL  IN  WIND 

TUNNEL 

between  the  cantilever  and  the  ideal  free-to-roll  configurations. 

The  experimental  results  for  the  free-to-roll  tip-tank 
models  are  as  inconclusive  as  the  theoretical  results.  The  four 
straight-wing  models  tabulated  in  Appendix  D  actually  represent 
one  model  which  was  repaired  after  each  flutter  so  that  small 
changes  in  the  experimental  parameters  occurred  (see  Appendix  D). 
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The  straight-wing  free-to-roll  model  fluttered  at  a  lower  Mach 
number  than  its  cantilever  counterpart,  while  the  opposite  is 
true  for  the  swept  wing. 


Model 

Root  Condition 

Mf 

ST-4a 

Cantilever 

1.43 

ST-4c-2 

Free-to-Roll 

1.32 

SW-7 

Cantilever 

1.44 

SW-7a 

Free-to-Roll 

1.92 

3 . 4  Comparison  with  Other  Experimental  Results 

In  assessing  the  value  of  any  experimental  program  of 
high  speed  testing  the  question  always  arises  as  to  whether  the 
results  obtained  are  representative  of  what  would  be  obtained  for 
a  full  scale  airplane,  or  merely  represent  the  peculiarities  of 
the  testing  facility  or  the  models  tested.  Little  or  no  full 
scale  airplane  flutter  data  is  available,  but  a  considerable  body 
of  flutter  experience  for  models  of  various  geometry  and  construc¬ 
tion  methods  is  available.  In  particular,  Reference  4,  5,  47,  48 
49  and  50  give  flutter  data  for  straight,  swept,  and  delta  wings 
whose  characteristics  are  similar  to  those  tested  under  the  pres¬ 
ent  program.  A  variety  of  model  construction  methods  are  repre¬ 
sented;  from  the  solid  plate  models  of  Reference  47  to  the 
wrapped  aluminum  foil,  foam  plastic  core  models  of  Reference  5-  A 
number  of  different  testing  facilities  have  been  used  to  obtain 
the  test  data.  Figures  3-17,  3-18,  and  3-19  are  plots  of  the  par¬ 
ameter,  Was*  ^  if versus  Mach  number  for  the  straight  swept 
and  delta  wing  models  of  this  report  and  for  similar  models  from 
References  4,  5,  47,  48,  49,  and  50.  This  parameter,  suggested  in 
Reference  47,  is  an  extremely  convenient  one  for  the  plotting  of 
flutter  results,  but  it  must  be  remembered  that  the  use  of  the 
factor  has  no  firm  theoretical  basis.  Figures  3-17,  3-18 

and  3-19  represent  a  collection  of  data  for  wings  that  are  similar 
but  not  identical.  Therefore  some  differences  in  the  data  are  bound 
to  occur.  Even  so  the  correlation  of  the  various  test  results  is 
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FIGURE  3-17  COMPARISON  OF  STRAIGHT  WING  RESULTS  WITH  SIMILAR  TESTS 
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quite  good. 

Figure  3-17  shows  the  results  for  the  straight  wing 
models.  The  data  of  Reference  47  is  too  sparse  to  draw  a  curve 
since  only  three  points  are  available  from  M  =  O.67  to  M  =  2.0. 
The  level  of  this  data  agrees  quite  well  with  that  of  the 
present  report.  The  data  of  References  4  and  48  agree  very 
well  with  each  other  in  the  transonic  range.  In  the  low  super¬ 
sonic  range,  the  data  of  Reference  4  shows  that  a  smaller  00 ^ 
is  needed  at  a  given  Mach  number  to  reach  the  stable  region  than 
for  the  models  of  this  report.  The  trend  of  the  data  obtained 
in  this  report  with  increasing  Mach  numbers  up  to  about  1.7  is 
somewhat  disturbing  since  it  indicates  that  higher  and  higher 
cu^'s  are  needed  to  prevent  flutter.  Since  no  data  was  obtain¬ 
ed  at  Mach  numbers  below  1.3,  it  is  not  possible  to  say  whether 
there  is  a  peak  in  the  boundary  at  about  M  =  1.0  which  is  higher 
than  the  peak  in  the  boundary  at  about  1.7.  The  data  of  Refer¬ 
ence  47  is  too  widely  separated  in  Mach  number  to  determine 
whether  or  not  a  peak  exists.  The  extended  curve  for  the 
straight-wing  models  of  this  report,  given  as  the  heavy  dash 
line  on  Fig.  3.17,  shows  a  slight  peak.  The  extension  is  based 
on  the  results  of  Ref.  48  and  some  unpublished  data.  Such  an 
extension  shows  that  an  airplane,  flying  at  constant  altitude 
or  constant  density  may  encounter  bending-torsion  flutter 
in  the  low  supersonic  range  at  about  M  =  1.6  after  having  passed 
through  the  transonic  range  without  flutter.  On  the  whole  the 
agreement  of  the  straight  wing  data  obtained  on  this  program 
with  that  of  References  4,  ^7,  and  48  is  quite  good.  There 
seem  to  be  no  major  discrepancies  that  cannot  be  explained  by 
the  differences  in  the  model  parameters. 

Figure  3. 18  shows  similar  curves  of  the  parameter, 
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vs.  Mach  number  of  the  swept  wing  models.  The 
agreement  between  the  present  data  and  that  of  References  li  and 
is  very  good,  The  three  sets  of  data  could  easily  be  repre¬ 
sented  by  a  single  curve.  In  spite  of  a  slight  tendency  towards 
increasing  ' s  necessary  to  prevent  flutter  at  the  higher 
Mach  numbers  covered  by  the  present  data,  the  critical  region 
for  flutter  appears  to  be  the  transonic  range. 


Figure  3-1°  shows  the  comparison  of  the  delta  wing  data 
of  the  present  report  with  that  of  Refs.  Lq f  and  50.  Less 
delta  wing  data  is  available  for  comparison  because  of  the  ap¬ 
parently  heavy  dependence  of  delta  wing  flutter  speed  on  second 
bending  to  first  torsion  and  first  bending  to  first  torsion  fre¬ 
quency  ratios  noted  in  Figs.  3-10  and  3 -If-  This  heavy  depend¬ 
ence  on  these  frequency  ratios  means  that  if  delta  wing  flutter 
data  is  to  be  compared  on  a  rational  basis  these  ratios  must  be 
similar.  Data  for  only  those  models  from  Refs.  ^7,  ^9,  and  50 
which  have  frequency  ratios  and  mode  shapes  similar  to  those  of 
the  delta  wings  of  the  present  program  have  been  used  in  Fig. 
3.19-  In  general  the  agreement  is  not  too  bad.  The  data  of 
Ref.  I17  shows  a  level  of  the  flutter  parameter  about  equal  to 
that  of  the  present  report  at  Mach  numbers  of  1.3  and  2.0.  The 
data  of  Refs.  ^9  and  50  indicates  that  there  is  a  peak  in 
boundary  around  M  =  1.0. 

The  curves  of  Figs.  3-17,  3-18,  and  3-19,  are  also 
extremely  useful  for  practical  airplane  design  in  that  they 
show  clearly  the  critical  conditions  for  flutter.  A  straight 
line  from  the  origin  M  =  0,  jfe/75 — )  -  0  of  the  curves 

can  easily  be  shown  to  be  a  line  of  constant  dynamic  pressure. 

A  straight  line  parallel  to  the  abscissa  represents  a  line  of 
constant  altitude. 


It  is  then  apparent  that  for  the  straight  wing  data 
show  in  Fig.  3.17  that  the  critical  flutter  region  for  a  con- 
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FIGURE  3-19  COMPARISON  OF  DELTA  WING  RESULTS  WITH  SIMILAR  TESTS 
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stant  altitude  airplane  condition  occurs  at  about  M  =  1.7.  At 
constant  dynamic  pressure  the  critical  flutter  region  occurs  at 
Mach  numbers  slightly  less  than  1.0. 

For  the  swept  wing  data  shown  in  Fig.  3-18  the  critical 
flutter  region  for  both  constant  altitude  and  constant  dynamic 
pressure  conditions  occurs  at  about  Mach  number  1.1. 

For  the  delta  wing  data  shown  in  Fig.  3-19,  the  peak  in 
the  boundary  at  a  Mach  number  of  about  1.1  is  very  close  to 
being  at  the  same  level  as  the  boundary  at  a  Mach  number  of  about 
2.0.  Therefore,  it  is  difficult  to  say  with  any  assurance 
whether  the  most  critical  flutter  region  at  constant  altitude  is 
in  the  transonic  region  or  at  higher  Mach  numbers.  The  data  of 
Ref.  47  indicates  that  the  boundary  continues  to  rise  beyond 
M  =  2.0  since  the  flutter  point  obtained  in  Ref.  47  for  M  =  3.0 

occurs  at  a  higher  value  of  7S  than  the  point  at 

M  =  2,0  shown.  It  would  appear,  then,  that  the  critical  region 
for  flutter  at  constant  altitude  lies  at  the  higher  Mach  numbers 
around  M  -  2.0  in  the  range  covered  by  Fig.  3-19-  This  con¬ 
clusion  should  be  used  with  care  because  of  the  scatter  in  the 

data  at  the  higher  Mach  numbers  and  the  general  scarcity  of  data 
on  equivalent  delta  wings  in  the  transonic  range. 

As  a  further  illustration  of  the  usefulness  of  the 
boundaries  of  Figs.  3-17,  3-18,  and  3-19  for  airplane  design,  a 
simple  design  example  is  discussed  below.  Figure  3-20  shows  the 
consolidated  flutter  boundaries  for  straight  wings,  drawn  from 
the  data  of  Fig.  3-17-  The  critical  constant  altitude  and  con¬ 
stant  dynamic  pressure  lines  are  shown.  The  following  airplane 
characteristics  and  desired  performance  are  given: 

-  20.0 

£>0.7i'  -  5-0  (ft.  ) 

M™.x.  =  2.0  at  30,000  (ft.) 
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FIGURE  3-20  CONSOLIDATED  STRAIGHT  WING  FLUTTER  BOUNDARY 
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The  level  of  wing  torsional  frequency  required  to  prevent  flutter 
for  the  high-speed  condition  will  be  determined,  and  then  the 
speed  and  dynamic  pressure  restrictions  that  occur  will  be  noted. 
The  standard  atmosphere  is  assumed. 

For  the  constant  altitude  condition  the  critical  value  of 
'jo  ns  is  0.28  at  M  =  1.7b  as  can  be  seen  from 
Fig.  3  .20.  The  resulting  value  of  <4*  needed  to  just  prevent 
flutter  is  9-8  cps.  For  this  value  of  cx>^  ,  the  airplane  must 
be  limited  in  both  speed  and  dynamic  pressure  to  prevent  flutter 
at  altitudes  below  30,000  ft.  For  example,  with  this  value  of 
,  the  Mach  number  for  flutter  at  sea  level  is  0.81  and  the 
critical  dynamic  pressure  is  670  pounds  per  square  foot. 

If  the  example  airplane  is  to  be  flutter  free  at  a  sea 
level-  Mach  number  of  1.1,  then  the  critical  value  of  the  para¬ 
meter  jT3—  77^0  is  0.23.  The  value  of  tO#  necessary 
to  prevent  flutter  for  this  condition  is  1^.8  cps,  about  50$ 
greater  than  the  previous  case. 
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SECTION  IV 
CONCLUSIONS 


It  is  possible  to  draw  a  number  of  conclusions  from  the 
data  presented  in  this  report.  The  experimentally  obtained 
flutter  boundaries  for  the  straight  wing  planform  can  be  compared 
with  the  boundaries  obtained  from  theoretical  studies.  This  com¬ 
parison  can  serve  to  indicate  in  what  Mach  number  ranges  theo¬ 
retical  calculations  can  be  expected  to  give  good  quantitative 
correlation  with  experimental  results.  Figure  3.-4  shows  graph¬ 
ically  that  good  correlation  between  theoretical  studies  and  ex¬ 
perimental  results  is  found  at  Mach  numbers  between  1.3  and  1.4. 
Outside  of  this  Mach  number  regime  the  correlation  appears  to  be 
progressively  poorer  with  the  theoretical  results  being  uncon¬ 
servative,  showing  too  small  a  region  of  instability,  above 
M  =  1.4,  and  conservative  below  M  =  1.3.  The  results  of  Table 
3.1  show  that  because  of  the  steeper  slope  of  the  theoretical 
curve  on  Figure  3-4  nearly  all  the  bare  wings  have  theoretical 
flutter  Mach  numbers  between  1.3  and  1.4. 

Some  other  conclusions  can  be  drawn  from  the  trends  of  the 
theoretical  and  experimental  data.  Both  the  theory  for  the 
straight  wing  and  the  experimentation  for  all  three  planforms 
show  that,  for  the  ranges  of  parameters  considered,  which  are 
typical  of  present-day  high-speed  aircraft,  the  region  of  in¬ 
stability  increases  with  the  following  independent  parametric 
variations : 

1.  decreasing  torsional  frequency 

2.  decreasing  mass  ratio 

3.  rearward  movement  of  the  tip-tank  center  of  gravity 

4.  decreasing  aileron  frequency 

h,  forward  movement  of  the  elastic  axis 
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These  items  are  listed  in  order  of  importance  based  on  the  ex-, 
perimental  results. 

In  general,  the  straight  wings  fluttered  in  bending-tor¬ 
sion;  the  swept  wings  sometimes  had  a  small  amount  of  second 
bending  in  the  flutter  mode.  It  is  difficult  to  say  what  modes 
were  involved  in  the  flutter  of  the  delta  wings  since  all  the 
motion  was  near  the  tip,  and  very  few  cycles  were  required  to 
damage  the  models.  One  delta-wing  with  aileron  fluttered  on 
injection  in  a  mode  in  which  aileron  motion  predominated,  but 
even  in  this  case  the  mode  shape  of  flutter  is  not  too  clear 
since  wing  torsion  and  bending  motions  are  involved  in  the 
flutter.  The  experimental  values  of  the  reduced  frequency  at 
flutter  were  much  lower  than  those  characteristic  of  subsonic 
flow.  Some  of  the  trends  noted  above  and  the  same  general 
changes  in  the  stability  boundaries  with  sweepback  as  those 
obtained  on  this  program  have  also  been  obtained  from  transonic 
tests . 


In  Section  3-^  the  data  obtained  on  the  present  program 
is  compared  with  that  obtained  in  other  tests  on  similar  models 
with  similar  parameters.  In  that  section  the  flutter  boundaries 
are  drawn  from  M  =  0.60  to  M  =  2.0  using  all  this  data.  Some 
scatter  exists  but  the  general  trends  may  be  noted.  Conclusions 
can  only  be  made  for  models  whose  parameters  compare  well  with 
chose  presented.  The  following  general  conclusions  may  be 
reached : 

(1)  For  straight  wings  with  the  parameters, 

Panel  =1.50 

center  of  gravity  at  about  50$  chord 
elastic  axis  at  about  42$  chord 


WADC  TR  5 h. - 1 1 3  ,  part  IT  68 


CONFIDENTIAL 


CONFIDENTIAL 


The  critical  flutter  condition  at  constant 
altitude  lies  at  Mach  numbers  of  about  1.7. 

(2)  For  swept  wings  with  parameters  similar  to 
the  straight  wings  given  above,  the  critical 
flutter  condition  at  constant  altitude 
probably  lies  in  the  transonic  region. 

(3)  For  60°  deltas  with  mode  shapes  and  natural 
frequency  ratios  similar  to  those  of  the 
subject  program,  the  critical  flutter  con¬ 
dition  probably  lies  near  M  =  2.0.  As 
noted  in  Section  3.4  this  conclusion  must 
be  used  with  care  because  of  the  scarcity 
of  delta  wing  flutter  data  for  models  with 
characteristics  similar  to  those  tested  on 
the  subject  program. 

(4)  For  all  those  planforms,  straight,  swept, 
and  delta,  the  critical  flutter  condition 
at  constant  dynamic  pressure  lies  in  the 
transonic  regime  at  Mach  numbers  close  to 
1.0. 
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APPENDIX  A 

THEORETICAL  FLUTTER  ANALYSES  FOR  STRAIGHT-WING  PLANFORM 

IN  SUPERSONIC  FLOW 

A. 1  Aileron  and  Bare  Wing 

Using  the  standard  nomenclature  of  Figure  A.l,  the  elements 
of  the  three-dimensional  flutter  determinant  of  Reference  26 
can  be  written  directly  for  our  case  by  a  straightforward  sub¬ 
stitution  of  the  two-dimensional  supersonic  aerodynamic  co¬ 
efficients  of  Reference  7-  A  comparison  of  the  expressions 
fo_  the  lift  and  the  moments  in  incompressible  flow  (Ref.  26, 
p.  25  ff)  and  supersonic  flow  (Ref.  7>  p.  7)  shows  that  the 
following  aerodynamic  terms  are  equivalent, 


Incompressible  Supersonic 
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FIGURE  A. 1  NOMENCLATURE  AND  CONVENTIONS  FOR  THEORY 
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Equations  (A.l)  through  (A. 9)  assume  that  the  leading  edge  of 
the  aileron  lies  on  the  hinge  line.  The  coefficient,  M,  ,  is 
taken  equal  to  1/2  (Ref.  '7,  p.  29).  With  this  equivalence,  the 
flutter  determinant  can  be  written  as  (Ref.  7,  pp.  62,  fa]), 

SEE 
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G  H  I 


-A  (N,  +  lNz) 

Eq.  (A 

~~A  (  f  l  N<t ) 

Ev  (A 

-  4  (  ^  N4^ 

El.  (A 

where 


A=  [}' (far)  (fau  )(•  +  ^)] / K fy)  A% 

o 

J  2  2 

-  y  <L<  l  Lz )  far  i  w 

B-  (  faulty 

~Y  f(LJ +t  Li)  fa )  fa]  ifa ,Jy 


Eq.  (A. 11) 
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Eq.  (a. 13) 


Eq.  (A. 14) 


Eq.  (A. 15) 


Eq.  (A. 16) 
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The  assumed  bending  and  torsion  mode  shapes  are  taken  as 
approximations  to  the  uncoupled  first  bending  and  first  torsion 

•mode  shapes  of  an  ideal  beam;  the  aileron  mode  is  assumed  to  be 
a  constant  rigid-body  motion,  i.e., 
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Eq.  (A. 20) 
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Eq.  (A. 21) 


Eq.  (A. 22) 


The  mass  parameters  of  the  models  were  built  to  be  typical 
of  actual  wings,  hence, 
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where 


A  is  the  ratio  of  the  tip  chord  to  the  root  chord.  The 
subscript,  o,  refers  to  conditions  at  the  root  chord  of  the 
wing.  In  the  case  of  the  aileron,  these  are  fictitious  values 
since  the  aileron  extends  over  the  outboard  half  of  the  wing 
only.  Taking  A  =  1/2,  introducing  the  dimensionless  variable, 


7  =  1 
L  JL 


Eq.  (A. 29) 


and  dividing  the  columns  and  rows  of  the  flutter  determinant  by 
appropriate  constants,  the  elements  of  the  flutter  determinant 
can  be  rewritten  in  the  following  dimensionless  form, 

A=  [l- (2£ f (%-fo+^Sl/L  jo-  %)Y*1  Eq.  (A>30) 
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Eq.  (A. 31) 


C  -  yj- 
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Eq.  (A. 33) 


Eq.  (A. 34) 
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where 

yU  is  the  wing  mass-density  ratio  in  supersonic  flow,  ^ oyi. 
It  should  be  noted  that  by  virtue  of  the  linear  taper  of  all 
dimensions  (see  Figure  1.1)  and  the  spanwise  mass  variations  of 
equations  (A. 23)  through  (A. 28),  the  parameters,^  ,  ,  iq  , 

f  ,  c  and  a,  are  constant  along  the  span. 

r 

The  integrals  multiplying  the  mass  parameters  are  definite 
integrals  and  can  be  evaluated  in  closed  form.  Their  values 
are : 


Eq.  (a. 35) 
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Eq.  (A. 41) 


Eq.  (A. 42) 


Eq.  (A. 43) 


Eq.  (A. 44) 


The  integrals  containing  the  aerodynamic  coefficients  cannot  be 
evaluated  directly  because  the  aerodynamic  coefficients  are 
functions  of  the  reduced-frequency  parameter  ,  ,  which  varies 

along  the  span. 


The  flutter  determinant  for  the  case  of  the  bare  cantilever 
wing  can  be  obtained  by  setting  the  aileron  terms  equal  to  zero. 
Thus,  the  determinant  is, 


A 

D 


Eq.  (A. 45) 


where  the  elements,  A,  B,  D,  E,  are  as  defined  for  the  wing- 
aileron  flutter  determinant. 
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A. 2  Tip  Tank;  Free-to-Roll  and  Cantilever  Conditions 

Assuming  that  the  motion  of  the  wing  at  flutter  can  be 
represented  by  the  superposition  of  the  fundamental  uncoupled 
bending  mode  (h) ,  the  uncoupled  torsion  mode  (<*)  and  the  rigid 
body  rotational  mode  about  the  wing  root  (0  ),  the  vertical  de¬ 
flection  of  the  wing  illustrated  in  Figure  A. 2  may  be  written 
as 


?  ( rx}y)t)  -  yd(i)  4-  +  (tc  -  ba.)c<(yti)  Ecl-  (A‘^6) 

Equation  (A. 46)  assumes  that  the  elastic  axis  of  the  wing  is 
straight  and  is  perpendicular  to  the  wing  root  chord. 


FIGURE  A. 2  AXIS  SYSTEM  FOR  WING  WITH  TIP  TANK 


With  the  introduction  of  assumed  modes,  the  kinetic  energy 
of  the  mechanical  system  becomes 
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KE-  t  m 


WU  ^  +  £  $  'h  lt)  +  <()], 


o  -t 

+  4-  Is  L§  «)J  a 


Jxtt 


Eq.  (A. 4?) 


+  £  J  +  Ui)  +(*-U)iUi)]  d-K 


where 


nij,(x  )  is  the  mass  of  the  tip  tank  per  unit  chordwise 
distance 

I  is  the  rolling  moment  of  inertia  of  the  wing 
support 

Each  motion  has  been  separated  into  a  time -dependent  variation 
(superscript—)  and  a  space-dependent  variation  (mode  shape). 
For  the  case  of  simple  harmonic  motion,  equation  (A. 47)  may  be 
differentiated  with  respect  to. time  and  integrated  along  the 
chord  to  give  the  following  results 


T§T~)  -  -  a)z«-  J Qjmtf  n  %  +  Ir^Of)  {(^  Jy 

f * 

+  *  -+  Ol,  +0/ Mr 

+  b  £  My  +  oL  £  Sr  J 


Eq.  (A. 48) 


t  J]  U Jrvtj)  {(,) 


-h  <*  j +  OJMr  -t-hMr 


Eq.  (A. 49) 


+  u  5t 
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where 

Mj,  is  the  total  mass  of  the  tip  tank 

is  the  static  unbalance  of  the  tip  tank  about  the 
elastic  axis  of  the  wing,  positive  nose  up 

IT  is  the  mass  moment  of  inertia  of  the  tip  tank  about 
the  elastic  axis  of  the  wing. 


The  assuming  of  fundamental  mode  shapes  to  describe  the 
flutter  motion  is  a  Rayleigh-type  approximation,  and  it  can  be 
shown  (Ref.  31 ,  P*  19)  that  the  variation  of  the  potential 
energy  with  respect  to  the  assumed  modes  may  be  written  as 


}(Pe) 


o 


Eq.  (A. 51) 


since  the  rigid-body  rolling  mode  cannot  affect  the  internal 
potential-energy  level  of  the  system,  and 


Eq.  (A. 52) 


3  (Pe) 

l P  ex) 


Eq.  (A. 53) 
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Because  the  potential  energy  is  a  function  only  of  the  elastic 
deflection  of  the  wing,  it  .is  convenient  to  choose  the  mode 
shapes  and  thus  the  frequencies  in  equations  (A. 52)  and  (A. 53) 
as  those  of  the  bare  wing  without  the  tip  tank.  This  assump¬ 
tion  implies  that  the  uncoupled  mode  shapes  are  the  same  for 
the  wing  with  and  without  the  tip  tank.  By  eliminating  the 
necessity  of  estimating  new  mode  shapes  and  frequencies  for 
each  change  in  the  parameters  of  the  tip  tank,  theoretical 
trends  for  variations  in  tip-tank  parameters  can  be  obtained 
with  a  minimum  of  effort. 


The  equations  of  motion  are  derived  by  applying  Lagrange's 
equation, 


J  /^)(ke))  __  2( KE )  ~^(PE)  _  ~ 

dt  /  Pg,  Pft  =  v; 


Eq.  (A. 54) 


th 


where  is  the  generalized  force  of  the  i  mode  and  represents 
all  the  forces  not  included  in  the  potential-energy  or  kinetic- 
energy  functions,  and  q^  is  the  magnitude  of  the  1  generalized 
coordinate  at  the  wing  tip  station.  Expressions  for  the  gen¬ 
eralized  force,  Q^,  are  obtained  from  virtual-work  considera¬ 
tions. (Ref.  31,  p.  56).  The  virtual  work  done  on  the  wing  as 
it  moves  through  the  virtual  displacements,  S0  ,  5'U  and  W  is 


Eq.  (A. 55) 


—  / 


^W~ Ja  -h  J  J7)  +  j SoC  Eq.  (A. 56) 


Sw*  Qg  s 6  +  QrSk  +  a  S5 


Eq.  (A. 57) 
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where 

L  is  the  lift  on  the  wing  per  unit  span,  positive  down 

I 

M  is  the  moment  about  the  elastic  axis  of  the  wing  per 
unit  span,  positive  nose  up. 

These  generalized  forces,  which,  for  a  system  with  zero 
damping,  consist  only  of  aerodynamic  forces  may  be  written  as 


<3? 

=  ®  +  9eh  ^ 

-h 

9e*  * 

Eq.  (A. 58) 

H 

"  9hQ®  +  Qhhh 

+ 

Qhu* 

Eq.  (A. 59) 

-  +  SU* 

QoioL* 

Eq.  (A. 60) 

th 

where  Q..  is  the  generalized  force  in  the  1  mode  per  unit 
1J  th 

displacement  of  the  j  mode.  For  the  aerodynamic  forces  act¬ 


ing  on  the  wing  and  on  the  tip  tank, 

W  /  +  Eq.  (A. 61) 

+  Lhr£  Eq'  <A-62) 

=  l  Eq.  (A.63) 

r  J’  __  _ 

^ he  4 £0  i +  Eq'  ^A‘6^ 

*Lkr  Eq-  (A-65) 

+  Eq-  (A'66) 

+  fy-T  Eq-  (A-67) 

Q«4-  l  +  'V  Eq.  (A. 68) 
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jf  *V/)  6'j)  Eq-  (A-69) 

where  explicit  expressions  for  and  M*  are  obtained 

from  Reference  7  (p.  7)  for  the  wing  and  those  for  L,  ,  , 

L Irn  rp 

R  and  M  from  equations  (B.2o)  and  (B.27)  for  the  xtip  tank, 

Llfjri  rj-1 


— -  yY)  Ltxj~t  /  \ 

I— h  ^  d  (t-t+i-Lz)  Eq.  (a .70) 

—  LCot 

L*  - -J?  *>'!>*  fa 

W  --  +if%  ) 

h 

Ah  =-^ 

I  Tr  Z  LOjt 

*-hr  ~  f°  '  C  Eq.  (A. 7*0 

j  —  — -  2  ^  ^  /  ■  1/  \ 

LoiT  -yo  VT  co  C  (>1T  -  l  Eq.  (A. 75' 


Eq.  (A. 71) 


Eq.  (A. 72) 


Eq.  (A. 73) 
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K  V 


VTCU^ 


LOlt  ,  V  ) 

(  «r  +  l  ~E^J 


Eq.  (A.v 


- -  L  y  p  7- 

=^Vr^e  K  +  -' 


f  +  .l'  +  £r/  iw-  (A.77) 


where 

V^,  is  the  /olume  of  the  tip  tank 

T  is  the  distance  that  the  geometrical  center  of  the 
tip  tank  lies  aft  of  the  elastic  axis  of  the  wing 

v  is  the  free-stream  velocity 

I,  is  the  volume  moment  of  inertia  in  pitch  of  the  tip 
tank  about  its  geometrical  center. 

To  simplify  the  theory,  the  wing  is  assumed  to  be  uniform 
with  constant  chord  (see  Figure  C.3).  If  the  non-rigid  bending 
and  torsion  modes  are  taken  as  the  fundamental  bending  and  tor¬ 
sion  modes  of  a  uniform  cantilever  beam  (see  Reference  32, 
pp .  of ,  66 ) ,  then 


I A  V  J3  *  4 


Eq.  (A. 78) 


=  0.338,931 


Eq.  (A. 79) 
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Using  these  results,  the  generalized  aerodynamic  forces 


Eq 
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.  (A. 80). 


.  (A.81). 


(A. 82)  . 

become 


.  (A. 83) 


.  (A.  84) 


.  (A. 85) 


(A. 86) 
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uzmSz  [-  y  ^  (am) ] 

■+■  —  /—  )4^ — i-  -A.)  7 

^  U^/l  *  Jz  JL  j  J 

9*0  =  c*z  e  /"  (Sj ^  +  tMz) 


Eq.  (A. 87) 


Eq.  (A. 88) 


7- (*&)(? +tz)]  ^  <A'89> 


+  ^  (W)(t  "T  t)] 


Eq.  (A. 90) 


Eq.  (A. 91) 


The  equations  of  motion  for  simple  harmonic  motion  of  the 
system  can  be  derived  by  the  application  of  Lagrange's  equation 
(Eq.  (A. 54)  with  the  following  results, 


f  2  e3 

L 


4-  CO 2  I3  -t-  COZ 


+  In  [o.2Q4j4I^7-^^Z+  cozMr-#  +- 

+  5  |?r 7=o 


Eq-  (A. 92) 
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Q  f  0-284)41  o>Z  m£Z -i-  (v*  Mr -l  4- 

Qhs>  j 

! 

+  h 

^  +  co-Mr  + 

<?44  1 
ei<J£  J 

Eq. 

(A. 93) 

~b  oi.  i 

f 0. 338,93  +  0)2Sr  + 

0^  7 

J 

i 

=  O 

Q  [^Tz 

^  5<  iZ  + 

7 

-h  V\  \ 

,0.338,03  Wir  -f  - 

1 

Q*kl 

tLu,tJ 

Eq . 

(A. 94) 

+  «[ 

(loz-u£^  4  ^  ^  ^ 

1 

eW  ' 

=  0 

Introducing  equations  (A. 83)  through  (A.91)  into  equations 
(A. 92)  through  (A. 94),  the  equations  of  motion  become,  after 
proper  non-dimensionalization  and  simplification, 

Kee^  +  Kh(r)  +  =0  Eq‘  (A'95) 

+  (t)  =  0  Eci-  (A*96) 

K^q  (0) +1(^(2)  +  J-)  =  0  Eq.  (A. 97) 
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where 


w  _  i  Is  hr  i  /  \  i 
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Eq.  (A. 99) 


Eq.  (A. 100) 


Eq.  (A. 101) 


Eq.  (A. 102) 


Eq.  (A. 103) 
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t  ^ 


+  -  4^  fa'*1  w  J 
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D<C><. 


1 

(-Hr  +.L)  Eq- 

,  A  k  / 
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0. 338,93  /  j  \ 
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yH  7I2?  (  h  JT)  q 

,  Ir  - 
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m 

Eq.  (A. 104) 


Eq.  (A. 105) 


/  ~Vy- 


f/Jk)z*.  1 1+  _L  -JL  Eq.  (a.  106) 

[{ tj  ej  yq 


where 


m  btf 


Eq.  (A. 107) 


Eq.  (A. 108) 


It  should  be  remembered  that  the  problem  has  been  formulated  so 
that  the  frequencies ,  and  ,  are  those  of  the  bare  wing 
without  the  tip  tank.  Structural  damping  can  be  introduced  by 
merely  replacing  the  frequency  ratio,  2  ,  by  (Ref.  31,  p.  196) 
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-(£-)('•♦>)  for 


K- 


hh 


Eq.  (A. 109) 


_n_ 


for  K 


«* 


Eq.  (A. 110) 


For  the  cantilever  condition,  0=0,  and  the  ecn’^Hons  of 
motion  degenerate  to 


(a  )  ~h  K*  \  )  ° 

—  I 

u  (Jl  x  [/  (-  i-)  *  o 

l  j?  /  ■  "<**  V~  JL  / 


Eq.  (A. Ill) 


Eq.  (A. 112) 


A. 3  Piston  Theory 

It  has  been  shown  in  Reference  33  that  the  application  of 
Piston  Theory  (Ref.  3*0  to  the  flutter  analysis  of  a  typical 
section  results  in  great  simplifications.  The  technique  is 
applied,  here  to  the  straight  tapered  wing  of  this  report. 

Using  the  first  bending  and  the  first  torsion  modes,  the 
flutter  determinant  from  Appendix  A.l  is 


A 

D 


Eq.  (A. 113) 
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where 


h-%krl'-(^tT  ?]  - 1. 


Eq.  (A. 114) 


2  _ 


J5  ~  3s  /*  ** 


Eq.  (A. 115) 


35*  /* 


X3 


Eq.  (A. 116) 


■33  -2 

6~60  /U  04 


(l  -  z) 


Eq.  (A. 117) 


and  the  aerodynamic  terms  are 


i 

(m, k)  =  f  (i,  +  i  Lz)0-  -r)  %A clZ 


Eq.  (A. 118) 


i 

j ~  f (l3 ^ 


Eq.  (A. 119) 


Is(M,k)-  f  (W  +i M*)#-  §)3fJZ 


Eq.  (A. 120) 
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V*,*}-  IT f  h 

If  any  one  of  the  conditions, 


Eq.  (A. 121) 


/v|2  »  |  MZk  »  \  P^k  »  I 

Eq.  (A. 122) 

holds,  Piston  Theory  can  be  applied,  and  the  aerodynamic  terms 
are  greatly  simplified.  From  Reference  33  (p.  6),  the  aerody¬ 
namic  coefficients  for  a  symmetrical  double-wedged  airfoil  at 
zero  angle  of  attack  are 


L  — 


M/e 


r 


Eq.  (A. 123) 
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-+  L  — 


Mkc 


F  ■+■ 


Mk 


[-  &  +  F(i  -  ?'X0jl 


Eq.  (A. 124) 


[-6+rfi-z-x.j] 

+  Fk  l~ZG  ^ 


E|.  (A. 125) 


Eq.  (A.  126). 


where 


F  =  /  + 


2T+  ( 
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(nr)2' 


Eq.  (A. 127) 
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7<p  is  the  fraction  of  the  chord  that  the 
elastic  axis  of  the  wing  is  behind 
the  leading  edge 


(A. 128) 


/T  is  the  thickness  ratio  of  the  airfoil 


For  a  flat  plate,  equations  (A. 126)  and  (A. 127)  reduce  to 


F  =  1 
G  =  0 


For  the  straight-wing  planform  of  this  report  (Fig.  1.1), 

y  «  K  Eq.  (A. 129) 

so  that 

k  =  kt  0-  i)  Eq.  (A. 130) 

and  the  integrals  of  equations  (A. 118)  through  (A. 121)  can  be 
evaluated  exactly, 


i 


Eq.  (A. 131).. 


I2(Wtk)  = 


Mk £  20 


F  + 


L  II 
Mko  120 


[-G  +F(i-:*S]  Eq-  (A. 132) 


6  // 

MK  120 


-hr 


Eq.  (A. 133) 
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For  a  given  set  of  wing  parameters  and  a  given  Mach  number, 
all  the  elements  of  the  flutter  determinant  (Eq.  (A.113)  are 
known  as  simple  functions  of  the  parameters,  £  and  kQ.  There¬ 
fore,  the  flutter  determinant:  can  be  solved  directly  for  its 
eigenvalues  £  and  kQ. 
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APPENDIX  B 

AERODYNAMIC  FORCES  CN  THE  TIP  TANK 

Expressions  for  the  lift  and  the  moment  on  the  tip  tank 
for  simple  harmonic  motion  in  supersonic  flow  can  be  obtained 
by  applying  the  results  of  Reference  35-  In  order  to  be  con¬ 
sistent  with  the  notation  of  Reference  35,  a  separate  set  of 
symbols  is  defined  for  this  appendix. 

The  basic  assumptions  of  the  analysis  are: 

(a)  The  medium  is  continuous. 

(b)  The  flow  is  frictionless. 

(c)  There  is  no  heat  transfer. 

(d)  There  are  no  shock  waves  of  finite  strength. 

(e)  There  are  no  body  forces  acting  on  the  fluid, 

(f)  The  undistrubed  medium  is  homogeneous. 

(g)  The  effect  of  wing  interference  on  the  tip  tank  is 
neglected. 

(h)  5  <<  l 


where 


where 


where 


S  is  the  body  fineness  ratio. 

i  $  «  I 

o  c*J  ( 2  c) 

k  is  the  reduced  frequency,  — jj — 

co  is  the  frequency  of  the  body  motion 

Zc  is  the  body  length 

U  is  the  free-stream  velocity 

(  M$)Z In  $  <<  \ 

M  is  the  free-stream  Mach  number. 
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(k) 

With  these  assumptions,  it  is  shown  in  Reference  35  (Eq .35) 
that  linearization  of  the  equations  of  motion,  which  results  in 
the  wave  equation,  in  the  neighborhood  of  the  body  for  the  gen- 
eral  transient  case  further  reduces  the  equations  of  motion  to 

=  °  Eq.  (B.l) 

where 

0  is  the  velocity  potential 

z  is  the  complex  coordinate,  x  +  iy  (see  Figure  B.l) 
z  is  the  complex  conjugate  to  z. 

•Equation  (B.l)  is  Laplace's  equation  in  the  x-y  plane,  and 
solutions  are  known  since  it  governs  the  case  of  irrotational 
incompressible  steady  flow.  Solutions  to  the  wave  equation, 
which  is  applicable  at  large  distances  from  the  body,  are  need¬ 
ed  for  the  determination  of  the  drag  only  and  need  not  concern 
us  here. 

Since  the  unsteady  supersonic  flow  at  any  instant  of  time 
is  approximated  by  an  incompressible  steady-state  flow,  ex¬ 
pressions  for  the  complex  potential  function  can  be  obtained 
by  resorting  to  classical  incompressible  aerodynamic  theory. 

For  the  case  at  hand,  i.e.,  a  body  of  revolution  in  unacceler¬ 
ated  forward  flight  but  with  transverse  motion,  the  flow  can 
be  considered  as  made  up  of  two  parts,  viz., 

(1)  the  radial  flow  due  to  the  free  stream,  U,  and  the 
change  in  radius  along  the  s-axis, 

(2)  the  transverse  flow  due  to  the  motion  of  the  body. 
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The  complex  potential  function  for  the  radial  flow  can  be 
obtained  by  considering  a  point  source  of  such  strength  that 
the  boundary  conditions  are  satisfied  on  the  surface  of  the 
cylinder  (see  Figure  B.2).  The  complex  potential  function  at 


y 


FIGURE  B.2  RADIAL  FLOW  DUE  TO  A  SOURCE 

the  point,  z p  for  a  point  source  at  the  origin  of  the  z^~ 
vector  is  (Ref.  36,  p.  212  ff  and  Ref,  37,  P-  196  ff ) , 

y  =  In  Eq.  (B.2) 

where 

W  is  the  complex  potential  function,  0  i  fs 
if/  is  the  stream  function 

H  is  the  volume  of  fluid  spilled  by  the  source  in  unit 
time . 
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The  strength  of  the  source,  H,  is  determined  by  applying  the 
boundary  condition, 

=  J£_  =  TJ 

Vr  "  ZVR  J*  Eq.  (B.3) 

where 

vr  is  the  radial  velocity 

R  is  the  radius  of  the  circular  cross-section 
at  station,  s. 

Therefore , 

H  =  S  '  17  Eq.  (B. 4) 

where 

S  is  the  cross-sectional  area  of  the  body. 


The  prime  denotes  differentiation  with  respect  to  s.  Hence, 

4-  in*. 

Eq.  (B.5) 

since  the  velocity,  U,  is  taken  equal  to  one  in  Reference  35. 


w  - 


The  complex  potential  function  for  the  transverse  flow  due 
to  the  motion  of  the  body  can  be  found  explicitly  in  Reference 
37  (p.  246). 


w=  - 


R Z  mtA 


Eq.  (B.6) 


where 

— — -Yfjis  the  substantial  derivative  ot  z 
Dt'31  g 


Therefore,  the  solution  to  equation  (B.l)  for  the  case  at 
hand  is 

T  -  5<s)  .  g'(s)  '  V_/r  ) 

V  ”  Zfr  (ni'-  qr  z,  Dt{  a>  Eq.  (B.7) 

which  compares  with  equation  106  of  Reference  35  except  that  it 
is  for  transverse  motion  in  unaccelerated  flight. 
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Reference  35  derives  an  explicit  expression  for  the  lift 
distribution  associated  with  a  more  general  complex  potential 
(Eq.  40  of  Ref,  35) , 


VY-S ac/5’/'0^'12  +  k  ^Sj  ^  +  ah(s>^  Z.  £q_  (B.8) 


Here  the  term,  b  (s,t),  represents  the  solution  which  allows 
the  inclusion  of  the  radiation  condition  (zero  disturbance  at 
infinity),  and  for  this  analysis  it  need  not  be  determined 
since  it  does  not  contribute  to  the  lift. 


In  terms  of  the  above  series,  the  lift  distribution  may  be 
written  as  (Eq,  92  of  Ref.  35), 

Eq.  (B.9) 

where 

F  is  the  lateral  force  per  unit  length  of  the  body 

b 

q  is  the  dynamic  pressure. 

The  circulation,  F  ,  is  given  by  equation  99  of  Reference  35, 

$(*) ]  Eq.  (B.lO) 

It  should  be  noted  that  the  lift  distribution  depends  only  on 
the  first  term  of  the  semi-infinite  series  of  equation  (B.8). 

It  now  remains  for  us  to  express  the  complex  potential  for 
the  case  at  hand,  equation  (B.7),  in  the  form  of  equation  (B.8), 
From  Figure  B.l, 


so  that 


2,  =  Z-  ~  Eq.  (B.  11) 

in  a,  Eq-  (B'12) 
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_  / 

2'  ’  2  -  H 

Expansion  of  equations  (B.12)  and  (B.13)  yields, 


Eq.  (B.13) 


2,  =  -fn  a 


Eq.  (B. 14) 


2, 


— -f 
Z.2 


Eq.  (B.15) 


Comparison  of  equation  (B.8)  with  equations  (B.7),  ( B . 14 ) ,  and 
(B.15)  shows  that  the  coefficient  of  interest,  a-p  must  be, 


(s) 

Z'iT 


kI(s) 

~¥ 


Eq.  (B. 16) 


Inserting  equations  (B.16)  arid  (B.lO)  into  equation  (B.9)  and 
performing  the  indicated  operations  give, 

Fs( s)=  Z%[ $(s)k  2^6)  -  560  60 

-i  £  &(s)  2j(s)  -  5>  Eq.  (B.17) 

- 

Simple  harmonic  motion  has  been  introduced  by  assuming 


_  .4 

*1  c 


Eq.  (B. 18) 
Eq.  (B.19) 


In  the  above  presentation,  the  dimensionless  parameters  defined 
in  Reference  35  have  been  used,  i.e.,  all  lengths  are  referred 
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to  the  body  length,  and  all  velocities  are  referred  to  the  free 
stream  velocity.  In  dimensional  form,  equation  (B.17)  becomes 

~  772  ~  K>(s)  -Zy  (3)-'  JJ-  S(s)gj(s) 

a  *  d  Eq.  (B.20) 

-2;^.  *)|e)  -&)  ’(,) J 


The  motion  of  the  body,  z  ,  can  be  expressed  in  terms  of 

D 

the  assumed  modes  of  vibration.  For  bending  of,  and  torsion 
about  the  elastic  axis  of  the  wing, 


Eq.  (B.21) 


Lcot 
e 

Eq.  (B. 22) 

where 

h  is  the  vertical  displacement  of  the  wing  elastic 
axis,  positive  down 

cC  is  the  twist  about  the  wing  elastic  axis,  positive 
nose  up. 


The  mode  shapes,  f^  and  f ^  ,  are  normalized  to  be  one  at  the 
wing  tip.  The  frequency  of  oscillation  in  equations  (B.21)  and 
(B.22)  is  <0,  not  k  as  above,  because  the  equations  have  been 
written  in  dimensional  form.  It  should  be  noted  that  time,  t, 
in  equations  (B.21)  and  (B.22)  is  dimensional,  whereas  time,  t, 
in  equations  (B.18)  and  (B.19)  is  dimensionless.  The  motion  of 
the  centroid  of  the  body  in  the  complex  plane  is 

=  -L  [h  + 
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The  subscript,  e.a.  ,  refers  to  the  elastic  axis  of  the  wing. 


The  amplitude  of  the  lift  is 


L  =  i  -  f: 


Eq.  (B. 24) 


and  the  amplitude  of  the  pitching  moment  about  the  center  of  the 
body  is 

(1 

/vj  =  lJFs  (j~c)  cfs 

°  Eq.  (B.25) 

When  equations  (B.20)  and  (B.23)  are  substituted  into  equations 
(B.24)  and  (B.25),  certain  integrals  appear  which  are  easily 
evaluated  if  one  remembers  that  the  body  under  consideration  is 
symmetrical  about  s  =  c  and  that  it  is  closed  at  both  ends.  The 
final  expressions  for  the  amplitudes  of  the  total  oscillatory 
lift  and  moment  are 


Eq.  (B.26) 


Mc  4.  =  nr  L  ■+  L/°  V  t'r  ^  h 

+■(/>«>*: t,  -bc^oV^ojnr  +/oVzvr)*  Ec*-  (B-27) 


where 


jO  is  the  air  density 

is  the  volume  of  the  body 

n^,  is  the  distance  that  the  geometrical  center  of  the 
body  lies  aft  of  the  elastic  axis  of  the  wing 

I  is  the  volume  moment  of  inertia  in  pitch  of  the  body 
about  its  geometrical  center, 
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The  above  results  agree  with  those  obtained  by  quasi-steady 
momentum  theory  for  incompressible  flow  (Ref.  38,  pp.  60,  61), 
as  expected. 
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SYMBOLS  FOR  APPENDIX  B 

am  Coefficient  in  expansion  of  W,  Eq .  (B.8) 

bQ  Coefficient  in  expansion  of  W,  Eq.  (B.8) 

c  Half  length  of  body 

D/Dt  Substantial  derivative, 

f^  Uncoupled  bending  mode  of  wing 

f^  Uncoupled  torsion  mode  of  wing 

F  Lateral  force  per  unit  length  of  body 

F  Amplitude  of  F  in  simple  harmonic  motion 

o  o 

h  Vertical  displacement  of  wing  elastic  axis,  positive  down 

h  Amplitude  of  wing  bending  mode  at  wing  tip  station 

H  Strength  of  source,  i.e.,  the  volume  of  fluid  spilled 

by  the  source  in  unit  time 

1-^  Volume  moment  of  inertia  in  pitch  of  the  body  about  its 

geometrical  center 

k  Reduced  frequency, 

L  Amplitude  of  oscillatory  lift,  positive  down 

M  Free -stream  Mach  number 
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M 


n 


T 


q 

R 


s 

t 

u 


V 


r 


V 


T 


W 


z 


z 


z 


g 


z 


cr 

o 


SYMBOLS  FOR  APPENDIX  B  (Contd.) 

Amplitude  of  oscillatory  pitching  moment  about  center 
of  body,  positive  nose  up 

Distance  the  geometrical  center  of  the  tip  tank  lies  aft 
of  the  elastic  axis  of  the  wing 

Dynamic  pressure, 

Radius  of  body  cross-section 

Coordinates  defined  in  Figure  B.l 

Cross-sectional  area  of  body 

Time 

Free-stream  velocity 

Radial  velocity 

Volume  of  the  body 

Complex  potential  function,  0  +  i  ^ 

Complex  coordinate, 

Complex  conjugate  to  z 

Coordinate  describing  motion  of  centroid  of  area  S 
Amplitude  of  z  in  simple  harmonic  motion 
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SYMBCLS  FOR  APPENDIX  B  (Contd.) 

Complex  coordinate  referred  to  center  of  area  (See 
Figure  B.l) 

ot  Twist  about  wing  elastic  axis,  positive  nose  up 

Amplitude  of  wing  torsion  mode  at  wing  tip  station 
P  Circulation 

£  Body  fineness  ratio,  i.e.,  ratio  of  maximum  cross- 

sectional  length  to  body  length 

p  Air  density 

0  Velocity  potential 

( Y  Stream  function 

aj  Frequency  of  the  body  motion 

subscript  e.a.  Elastic  axis  of  the  wing 
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APPENDIX  C 

ARITHMETICAL  EXAMPLES  ILLUSTRATING  USE  OF  THEORY 


C . 1  Bare  Wing,  Incompressible 


Expressions  for  the  elements  of  the  three-dimensional 
flutter  determinant  in  incompressible  flow  are  given  on  page  65 
of  Reference  26.  Thus,  for  the  combined  bending-torsion  case, 
the  flutter  determinant  is  given  by 


A 

D 


Eq.  (C.l) 


The  coefficients  of  the  flutter  determinant  are  the  same  as  for 
the  bare  wing  in  supersonic  flow,  as  given  by  equations  (A.30), 
(A. 31),  (A. 33)  and  (A. 34) ,  except  that  the  aerodynamic  coeffi¬ 
cients  are  the  incompressible  coefficients  tabulated  in  Refer¬ 
ence  26,  and  jji,  is  the  incompressible  wing  mass  ratio, 

/ 

\m 

'h' /o  Eq.  (C.2) 

The  mass  parameters  are  assumed  to  vary  in  the  manner  expressed 
by  equations  (A. 23)  through  (A. 28). 


The  eigenvalues  of  the  two-degree-of-freedom  system  were 
chosen  as  H  and  •  Since  the  wing  is  tapered,  the  aero¬ 

dynamic  coefficients  are  functions  of  spanwise  location,  so  that 
the  aerodynamic  integrals  must  be  evaluated  numerically  for 
every  value  of  k  .  Reference  points  for  the  numerical  integra¬ 
tions  were  chosen  as  0,  20,  40,  60  and  100$  span.  Simpson's 
Rule  was  used  between  20$  and  100$  span  since  it  requires  an 
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odd  number  of  stations,  while  the  Trapezoidal  Rule  was  used  be¬ 
tween  0$  and  20$  span.  In  evaluating  the  aerodynamic  coeffi¬ 
cients,  L^,  L  ^  and  M  ^  ,  for  the  specified  spanwi.se  stations, 
it  was  found  that  interpolation  between  tabulated  values  of  re¬ 
duced  frequency,  k,  could  be  avoided  by  using  the  basic  ex¬ 
pressions  given  on  page  29  of  Reference  26  and  by  discriminately 
choosing  the  values  of  kQ.  Convenient  tables  of  Theordorsen 1 s 
function  can  be  found  on  page  3^2  of  Reference  39. 

For  kQ  =  0.20,  the  aerodynamic  integrals  are 


Jth  (i-  -  -0JS%  46  -0.343,940 

Jo[L<  -Lj,(4 +«)]/>-%)  ¥ 

—  (o. /39,5k  3-  0.433,36  t 


Eq.  (C .  3 ) 


Eq.  (C.4) 


O.OIOSH,  -h  0.3/3,  14  i 


Eq.  (C.5) 


,t° 


I  fa  -  Mh  (i  +«)  -  u  (i  ^)-hLk  a  4a)2](/-  *)  fjn  s 


Z.'o5,31  -  Q.577,/7  i  E9-  (c-6) 


Using  the  parameters, 


/•  =  65 


x*  =  0.114 
a  =  -0.114 
V  =  0.237,40 

§h  “  ^  =  °. 


Eq.  (C.7) 
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the  flutter  determinant  can  be  expanded,  and  the  real  and  imagi¬ 
nary  parts  can  be  set  equal  to  zero,  giving 


4.051,38  (3.473,00  (zk)*  Z 


-3.330,33  2.  -h  iS-.j-OQjGf  ~  Q 


Eq.  (C.8) 


2.6'&Z)Q5'(^fl-t-0.855J842-3.77lJ2O-O  Ecl-  (c-9) 


From  the  definition  of  Z, 


4  fe„/IT  Eq.  (C.10) 

Therefore,  the  following  table  can  be  calculated  for  the  real 

solutions  of  equations  (C.8)  and  (C.9). 

/  \  1  v+ 

Z  Ui7  *0  4 

4.026,63  C.031s25  5  2.492 


0.778,04 


C. 031,25 
1.545,28 


5.669 


Repeating  the  steps  for  different  values  of  kQ  results  in 
the  curves  shown  in  Figure  3-7. 


C , 2  Bare  Wing,  Supersonic 

The  flutter  determinant  for  the  bare  wing  in  supersonic 
flow  is  given  by  equation  (A. 45).  The  parameters,  Z  and  (~ 
were  again  chosen  as  eigenvalues.  The  method  of  solution  is  the 
same  as  for  the  bare  wing  in  incompressible  flow  except  that  a 
value  of  Mach  number  as  well  as  frequency  must  be  chosen.  For 
high  values  of  frequency,  the  aerodynamic  coefficients,  Lp  Lg , 
L^,  Lp  Mp  Mg,  ,  and  M^  were  obtained  from  Reference  7  and 
linear  interpolation  was  used.  For  low  values  of  frequency,  it 
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was  found  that  linear  or  three-point  interpolation  between  the 
tabulated  values  of  the  reduced-frequency  parameter,  00  could 
not  be  used  to  evaluate  the  aerodynamic  coefficients  with  suffi¬ 
cient  accuracy  because  they  are  highly  non-linear  in  the  low- 
co  range.  Also,  it  has  been  pointed  out  in  Reference  40  that 
a  few  errors  exist  in  the  tables  of  Reference  7  associated  with 
the  smallest  values  of  co 

Therefore,  for  low  values  of  co>  (  caj  £-1/2,  approximately) 
and  for  Mach  numbers  not  tabulated  in  Reference  7  (M  =  ^  for  the 

O  J 

bare  wing,  M  =  75-  for  the  wing  with  tip  tank),  resort  was  made  to 
the  basic  expressions  for  the  aerodynamic  coefficients  given  in 
Reference  7  and  to  the  method  of  Reference  41  for  evaluating 
the  basic  function,  f  ^  (M,  ) .  It  should  be  noted  that  the 

method  of  Reference  41  allows  the  evaluation  of  the  function, 

,  directly  without  the  need  of  recursion  formulae  and  that 
the  labor  required  for  a  given  accuracy  decreases  with  decreas¬ 
ing  to  .  More  extensive  tables  of  the  aerodynamic  coefficients 
than  those  of  Reference  7,  particularly  for  the  argument  of  Mach 
number,  are  given  in  Reference  23.  This  reference  became  avail¬ 
able  after  the  computations  for  the  prescribed  Mach  numbers  were 
well  underway.  The  same  method  of  integrating  the  aerodynamic 
coefficients  across  the  span  as  outlined  above  in  Appendix  C.l 
was  used. 

For  M  =  U)  =  0.72835  and  xq  =  0.443,  the  aerodynamic 
integrals  are 

j A, =  0.001,68  +  0.601, so i  Ei-  <c-n 
J  (L,  *  cL, S3  -  0.380, 07 i  E<1-  (c'12> 
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[m,  +  L  Mj,-  =  Q.OZZ,<oZ  +0.047,  asi 


Eq.  (C. 13) 


,  i  ■ 

(M3  +  i  fffi?  0.372, 06  -  0.02 1 , 96  i 


Eq.  (C. 14) 


Noting  that  the  wing  mass  ratio  is  now  that  for  the  supersonic 


case , 


r  r. 

/•'/4 


Eq.  (C.  15) 


for  the  same  parameters  given  in  equation  (C.7),  expanding  the 
flutter  determinant  and  setting  the  real  and  imaginary  parts 
equal  to  zero  give 

(<£)  0.477Sl(-^)i  -0.381,111 

Eq.  (C. 16) 

+  0.939,43  =  (9 

2-  -  O.  l&030&(%+f 1-1.007,85  Eq-  (c-17) 

The  real  solutions  of  equations  (C.16)  and  (C.17)  are 


1.031,19 


0.129,64 


From  the  definitions  of  Z  and  <^J0 

£  =  9.384,  =  5-300. 

O  b0  CO U 

Repeating  the  steps  for  different  values  of  M  and  oo  re¬ 
sults  in  the  curves  shown  in  Figure  C.l  (h) .  The  parameters 
used  for  the  remaining  curves  of  Figure  C.l  were  chosen  so  that 
the  complete  experimental  range  of  the  models  tested  was  covered. 
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FIGURE  C. 1(a)  FLUTTER  COEFFICIENT  AND  REDUCED  FREQUENCY  VERSUS 
FREQUENCY  RATIO  FOR  BARE  STRAIGHT-WING  PLANFORM, 
e.a.  AT  39.3^  CHORD ,  =  30 
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FIGURE  G . 1(b)  FLUTTER  COEFFICIENT  AND  REDUCED  FREQUENCY  VERSUS 
FREQUENCY  RATIO  FOR  BARE  STRAIGHT-WING  PLANFORM , 
e .  a.  AT  39-3$  CHORD,  yU.  =  35 
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FIGURE  C. 1(c)  FLUTTER  COEFFICIENT  AND  REDUCED  FREQUENCY  VERSUS 
FREQUENCY  RATIO  FOR  BARE  STRAIGHT-WING  PLANFORM, 
e . a.  AT  39-3$  CHORD, =  65 
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FIGURE  C.l(d)  FLUTTER  COEFFICIENT  AND  REDUCED  FREQUENCY  VERSUS 
FREQUENCY  RATIO  FOR  BARE  STRAIGHT-WING  PLANFORM, 
e .  a.  AT  39. 3^  CHORD ,yu,  =  75 
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FIGURE  C.l(e) 
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FLUTTER  COEFFICIENT  AND  REDUCED  FREQUENCY  VERSUS 
FREQUENCY  RATIO  FOR  BARE  STRAIGHT-WING  PLANFORM, 
e.a.  AT  44.3$  CHORD ,  jx  =  30 
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FIGURE  C.l(g)  FLUTTER  COEFFICIENT  AND  REDUCED  FREQUENCY  VERSUS 
FREQUENCY  RATIO  FOR  BARE  STRAIGHT-WING  PLANFORM, 
e.a.  AT  44.3$  CHORD,  u.  =  55 
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In  all  cases,  not  enough  real  solutions  were  obtained  for  Mach 
number  ^  to  define  the  curves  in  the  range  of  interest,  so  that 
M  =  j  was  used  when  information  below  M  =  was  considered 
necessary.  The  curves  of  Figure  C.l  were  then  cross-plotted  to 
obtain  those  of  Figure  3  •  1  • 


It  should  be  noted  that  the  problem  has  been  formulated 
completely  on  a  dimensionless  basis  so  that,  in  the  computations, 
there  is  no  need  for  using  a  curve  of  velocity  versus  Mach  num¬ 
ber  for  the  wind  tunnel  or  for  the  atmosphere  (see  Figure  3*2). 


C . 3  Wing  with  Aileron 


The  coefficients  of  the  determinant  for  combined  bending- 
torsion-aileron  flutter-  are  given  by  equations  (A. 30)  through 
(A. 38).  In  this  case  it  was  decided  that  the  maximum  amount 
of  information  for  a  given  amount  of  effort  could  be  obtaine.' 
by  choosing  the  frequency  ratios,  Z  and  ,  as  eigenvalues. 

The  values  of  the  integrals  multiplying  the  mass  parameters  are 
given  by  equations  (A. 39)  through  (A. HH) .  When  the  integration 
extended  over  the  entire  wing,  the  integrals  containing  aero¬ 
dynamic  coefficients  were  evaluated  numerically  for  each  value 
of  the  Mach  number,  M,  and  the  reduced-frequency  parameter,  to  , 
as  above  for  the  bare  wing.  The  integrals  containing  aerody¬ 
namic  coefficients  but  with  the  integration  extending  only  over 
the  aileron  were  evaluated  by  assuming  that  to  could  be  taken 
as  a  constant  for  a  representative  section  of  the  wing.  This 
representative  section  varied  for  each  integral  and  was  found 

by  considering  a  weighted  value  of  the  integral.  To  illustrate, 

IQ  /  to 


/ 0-g)  Eq-  (c.i8) 


The  complex  coefficient, 


+  1 


l6)~ 
b  z 


was  evaluated  for  the 
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value  of  <o  taken  at  station  £  ,  where  £  was  located  at  the 

centroid  of  the  function, 


r?  —  — — 
L  r'° 


~~To  ”  ‘ 

/  (,-?) 


Eq.  (C. 19) 


This  method  would  give  exact  results  if  the  aerodynamic  coeffi¬ 
cients  varied  linearly  over  the  inter-  .1  of  integration. 

For  M  =  cc>o  =  0.849,76  and  xQ  =  0.443,.  the  aerodynamic 
integrals  are, 


f(l<  ^  <-  Lz)(\-  Jr)  = 


0.06/,  8S  -/-  0.6-//,  OS  i 


Eq.  (C . 20) 


f  t  ^ 

J  (l3  +lQ(/-  £)  =  3.03/,34  -  0.32S ;  92  l 


Eq.  (C.21) 


2  3  2 

J  (L5'hi  L^)(l-  i)  =  0.100, 76  +  O. OOOJ6  0 


Eq.  (C.22) 


(M,  +iMt)(-z)  =  0.0Z2i  01 +  0.03Q,8ZL 


Eq.  (C.23) 


j  (rt3  =  O.ZSl]  90  -  0. 013 ,  QOi 


Eq.  (C.24) 
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j  (Ms.*  §)  ?c,1  =  0.522,20  -h  0.000,841  Eq.  (C.25) 


i 

j (M  +  l  =  0.004,16  +  O.Ol7,Z7C 


Eq.  (C.26) 


f  (N3  =  0.104,33  -  0.009,37c 


Eq.  (C.27) 


j  (s/s  +lMc)(i-%)  cf{  =  0.153,3(0  +  0.  000,35 c 


Eq.  (C.28) 


Using  the  parameters, 

wing  plus  aileron 

aileron 

/C  =  65  (^r) 

V  =  0.114 

x^  =  0.012 

P 

a  =  -  0.114 

c  =  o .  6o 

r  ot^=  0.237,40 

(S)Vo-10 

r  2=  0.003,24  Eq.  (C.29) 
f3 

8h  “  “  0 

g  =0 

v 

expanding  the  determinant  and  setting  both  the  real  and  the 
imaginary  parts  equal  to  zero  give 


- 0.00  3 


~  3  Afc  f  Z 

Z  +  0.035,89^/1 


Eq.  (C.30) 


0.032,83  h^j  t  -0.0%£6Z  -  0.316,93  i.  -  0.398,4}  =  O 
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FIGURE  C . 2  FLUTTER  COEFFICIENT  AND  REDUCED  FREQUENCY  VERSUS 
FREQUENCY  RATIO  FOR  STRAIGHT  WING  WITH  AILERON 
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0. 049,  (8 (■%£■')  f  +  O'OSi , 

—  0.000 j  63 ■£ 2  —  0.454,292.  +  0.534,83  =  O 


Eq.  (c .31) 


Combining  these  two  equations  results  in  a  quartic  which  can  be 
solved  by  Graeffe's  root-squaring  process  (Ref.  42,  p.  484). 
Only  the  real,  positive  solutions  of  equations  (C.30)  and  (C.31) 
have  any  physical  significance;  they  are 


1.159  1.772 


From  the  definitions  of  Z  and  , 

£  =  3.332.  =  3.095 

•  0  4  Ok 

By  repeating  this  procedure  for  other  values  of  M  and  (x>0  , 
it  is  possible  to  obtain  the  curves  shown  in  Figure  C.2.  The' 
curves  of  Figure  3-12  were  obtained  by  cross-plotting  those  of 
Figure  C.2. 

C . 4  Wing  with  Tip  Tank,  Cantilever 

The  equations  of  motion  for  the  cantilever  wing  with  tip 
tank  are  given  by  equations  (A. Ill)  and  (A. 112).  For  simplifi¬ 
cation,  the  wing  used  for  the  theoretical  calculations  was 
assumed  to  have  constant  chord  and  constant  mass  properties  in 
the  formulation  of  the  theory.  The  values  for  the  chord  and 
the  mass  properties  were  taken  as  representative  of  the  seventy- 
percent-spanwise  station  of  the  actual  planform  (See  Figure  1.1). 
The  seventy-percent-spanwise  station  was  chosen  because  experience 
has  shown  this  location  to  be  best  for  aerodynamic  purposes.  The 
tip-tank  parameters  were  unchanged  from  the  actual  values  so  that 
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the  ratio  of  the  tip-tank  mass  to  the  total-wing  mass  was  great¬ 
er  for  the  theoretical  model  than  for  the  actual  models.  The 
fact  that  the  theoretical  bare  wing  is  lighter  and,  as  a  con¬ 
sequence,  would  have  higher  vibratory  frequencies,  is  partially 
rectified  by  having  a  greater  proportion  of  the  mass  near  the 
wing  tip.  The  configuration  of  the  theoretical  model  is  given 
in  Figure  C . 3 • 

It  was  decided  that  the  most  useful  information  could  be 
obtained  by  leaving  the  static  unbalance  of  the  tip  tank  and  the 
first  torsional  frequency  of  the  bare  wing  as  unknowns  in  the 
flutter  equations.  Thus,  the  parameters,  !3L  and  Z,  were  taken 
as  the  two  eigenvalues  of  the  flutter  determinant. 


FIGURE  C.3  TIP-TANK  MODEL  FOR  CALCULATIONS 
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The  parameters  chosen  for  the  analysis  are 


m  =  0.020,31  slugs/ft 

yU  =  65  (if) 


( 


CO. 


x^=  0.114 


r  =  0.237,40 
-j2'=  0.10 


a) 


vT  -  0.008,19  ft3 


T 

Mr 

h 


0.030,47  slugs 
0.000,46  ft.5 


IT  =  0.001,77  slug-ft‘ 


n 


T  =  0.114 


Eq.  (C.32) 


With  these  parameters,  the  coefficients  of  the  flutter  deter¬ 
minant  become 


=  2./70,72  -0.025",  002  -  0.004,9 1  (l,  +lLz ) 

Eq.  (C.33) 

=  0.036,7 Z  -h  $T  -  0.006,66 (L3+cL4)-  0.000,10 

Eq.  (C.34) 

• 

~  0*058,12.  +  5r  -  0.006,66  (d,  ■+  0.000,10 

Eq.  (G.35) 

K^  =  1.6^1^3-0.118,702  -Q.009,Q2(M3i;M4)+ 

k 

Eq.  (0.36) 

For  M  =  to  =  0.34  and  xQ  =  0.443-,  the  aerodynamic 

are 

coefficients 

L,  +  tlL  =  adzt,s 7f  n. i47, 60c 

Eq.  (C.37) 

L.3  4- 1/.^  =  126.373,49-3.713,93  i 

Eq.  (0.38) 

M,  +  iMz  =  0.4  06, 87  3-  I.I9Z,  13  i 

Eq.  (C.39) 
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FIGURE  C.4(a)  FLUTTER  COEFFICIENT  VERSUS  TIP-TANK  STATIC-UNBALANCE 
PARAMETER  FOR  STRAIGHT-WING  PLANFORM,  NO  STRUCTURAL 

DAMPING 


WADC  TR  54-113,  Part  II  ■ 


FIGURE  C.5  FLUTTER  COEFFICIENT  AND  REDUCED  FREQUENCY  VERSUS  TIP 


TANK  STATIC -UNBALANCE  PARAMETER  FOR  STRAIGHT -WING 


PLANFQRM,  STRUCTURAL  DAMPING  INCLUDED 
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M3  + «-  Ml  =  13. 1G3>, 8G  -  O.lAQ^zi 


Eq.  (C. 40) 


Expanding  the  determinant  and  setting  the  real  and  imaginary 
parts  equal  to  zero  give 


0.001,97  -  0.297,  iZ  £  +  0.160 ,  f2  5r 


—  z. 


-  5T  +  3 .4  3^73  =  O 


Eq.  (C.41) 


0.06 3,/S-Z-  G.53Q}SSSr  -  0,134 j8Z  -  O 


Eq.  (C.42) 


The  real  roots  of  equations  (C.41)  and  40-42)  are 

Z  Srn 


14.298,94 


0.312,28 


Hence , 

i  =  11.534,  =  3.050 

K  b  cJ* 


Choosing  other  values  for  M  and  oj  and  repeating  the  steps 
above  result  in  the  curves  shown  in  Figure  C.4.  A  small  amount 
of  structural  damping  was  included  in  the  analysis  to  obtain 
the  curves  of  Figure  C.5.  The  curves  of  Figures  C.4  and  C.5 
were  then  cross-plotted  to  obtain  those  of  Figures  3-13  and  3-14, 
respectively. 
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C . 5  Wing  with  Tip  Tank,  Free-to-Roll 

The  equations  of  motion  for  the  straight  wing  with  tip  tank 
are  given  by  equations  (A. 95)  through  (A. 97).  Again,  the  wing 
is  assumed  to  have  constant  chord  and  constant  mass  properties 
(see  Figure  C.3),  and  the  quantities,  S^,  and  Z,  were  chosen  as 
eigenvalues.  With  the  parameters  of  equation  (C.32)  and  the 
value  of  the  mass  moment  of  inertia  of  the  roll  support  about 
the  roll  axis, 


I  =  0.001,50  slug-ft2 
s 

the  coefficients  of  the  flutter  determinant  become 


=  2.400,93- C>.OQej5'S'(L,  -f  lL?)  Eq.  (C.43) 

“  Z  ZQS)  /3  -  0.005,59  (L,  -i  l  Lz)  Eq<  (C*44) 


-  4D46,2©  +-5y  -  o. 007, 9&(l3  +  c  i4 J  - 0.000,10 


Eq.  (C.45) 


=  Z-ZOS,  13- O'OOSjSS  Cl, +iLz)  £q.  (c-46) 


2MOJ£-O.OZSJOOZ  -  0.O04fi\ (l,  +cl2)  Eq.  (C.47) 


WADC  TR  54-113,  Part  II 


138 


CONFIDENTIAL 


1(^  =  0.038,72  +  5r -0. 006^6 (l3+lL4)  ~  O.OOOJO  -~- 

Eq.  (C.48) 

0.04(6,28  4^  "0.007,36(^+6^)+ 0.000,70  ~ 

Eq.  (C. 49) 

K<k  =  °‘OSb>12-  +^>T  -  0.006,66  (m,  +0  M2)4  0.000,70 

Eq.  (C.50) 

K^=  1.641  99-  0./I8, 702-0.009,82 CfHj+cM*)  4  --0Q^iZ^ 

Eq.  (C.51) 

For  M  =  2, 03=  0.10  and  xq  =  0.443,  the  aerodynamic 
are  found  from  Reference  7  (pp.  8  and  24)  to  be 

coefficients' 

L,  +  ‘Li  =  <^92,  II  +  !S~,  386,40  J 

Eq.  (C.52) 

U+cLj  =4/0.262- 3. 363,85*4 

Eq.  (C.53) 

=  O.OS£T  66  4-  /.749,  452 

Eq.  (0.54) 

H3  4-  =  46.653  4  3. 039,95c 

Eq.  (C. 55) 

Expanding  the  flutter  determinant  and  separating  the  real  and 
imaginary  parts  give 


0.63*3,64  -  Oj4e>,3l2L+  0.001,  I4?z  1-  0.4  3lj7Z  5r 

_  _  Eq.(C.56) 

-  0-169,33  5r  -  0.079,372  5r  +  0.025",  00  ?Sr  =0 
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FIGURE  C.6  GRAPHICAL  SOLUTION  OF  SIMULTANEOUS  EQUATIONS 
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-0/1 69,97  4-  0.082,24  £  -O.OOZP92.2  -0.07^88  Sr 

-2  _  Eq.  (C.  57) 

+  0,643,67  5-7-  4-0.003,21  £  Sr  =  0 

This  set  of  simultaneous  equations  is  most  easily  solved  graph¬ 
ically  (see  Figure  C.6).  Additional  accuracy  can  be  obtained 
by  iteration.  The  solutions  are 


z 

1 

IT 

V4 

b  cO^ 

5.300 

-1.074 

26 . 67 

11.583 

7.270 

-0.405 

26.67 

9.890 

19.46 

0.479 

26.67 

6.045 

23.75 

2.191 

26.67 

5.472 

v<  1 

Therefore,  four  points  on  the  curves  of  and  r-  versus 

K  _ 

have  been  determined.  By  assuming  different  values  of  co  , 
the  process  can  be  repeated  until  the  curve  for  M  =  2  has  been 
defined  (see  Figure  C.7).  As  explained  in  Appendix  C.2  above, 
interpolation  of  the  tables  in  Reference  7  is  inaccurate,  so  the 
aerodynamic  coefficients  should  be  tabulated  by  the  method  of 
Reference  41  (or  Reference  40  for  low  co  )  for  values  of  M  and 
oo  not  found  in  Reference  7-  The  curves  of  Figure  3*15  were 
obtained  by  cross-plotting  those  of  Figure  C.7* 

C.6  Piston  Theory 

The  coefficients  of  the  two-degree-of -freedom  flutter  de¬ 
terminant  for  the  straight-wing. planform  of  this  report  (Fig, 
1.1)  are  given  by  equations  (A. 114)  through  (A, 117).  The  aero¬ 
dynamic  terms  have  been  derived  for  the  tapered  planform  with  a 
symmetrical  double-wedged  air  foil  section  at  zero  angle  of 
attack  in  equations  (A.131)  through  (A. 134). 
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C.7(a)  FLUTTER  COEFFICIENT  VERSUS  TIP-TANK  STATIC -UNBALANCE 
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FIGURE  C.7(b)  REDUCED  FREQUENCY  VERSUS  TIP-TANK  STATIC -UNBALANCE 
PARAMETER  FOR  STRAIGHT-WING  PLANFORM,  FREE-TO-ROLL 
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With  the  following  parameters, 


> 

II 

CTn 

/ =  1.4  (air) 

x.  =  0. 114 

/V  =  0.06 

rj  =  0.237,^0 

M  =  2.0 

f—f=  0.25 

x  =  0.443 

0 

8h  =  S*  -  0 

the  coefficients  become 

A  =  69  -  0.878,672.  -  O.OSSjQA  ~H~ 

f  S  Q 

B  --  0.331,  S3  -  O.Q7S,6S -  0.OO/.97  -~- 

k 

D  *  0.33^5-9  -  0,00^97  -j- 

_  It 

E  -  0.112.,  10  -  0.1(2^0  2.  -  0.002,87  ■j£a  -  0.014,5*6 

Expanding  the  determinant  and  setting  the  real  and 
equal  to  zero  give 


-o.sn.ei  -v  o.5-46j95'Z  +- 


7  1 

I03  ^ 
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Eq.  (C.59) 

Eq.  (C.60) 

Eq. (C .61) 

Eq.  (C . 62 ) 
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imaginary  parts 

Eq.  (C.63) 
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Q.&Z5  70  TLZ  -  3./2B,  49Z.  +  £l 

J  to3  kZ 


E4. 


+ 


14.  /SS.43 


(03  k 


2  -t~  2.392,84  *  0 


The  real,  positive  solutions  to  this  set  of  equations  are 

Z 


1 

¥ 


1.657 


0 

7.651 


whence 


* 


=  5.9^ 


Repeating  the  process  for  other  parameters  (Eq .  C.58) 
suits  in  the  curves  of  Figure  3.6. 
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APPENDIX  D 

DETAILED  TABULATION- FOR  DATA  FOR  MODELS  TESTED 

D. 1  General 

In  this  appendix,  a  detailed  tabulation  of  data  is  present¬ 
ed  for  all  the  models  tested.  The  number  of  significant  figures 
quoted  is  consistent  with  the  experimental  accuracy.  See  Refer¬ 
ence  14  for  a  discussion  of  the  testing  techniques  used  to  ob¬ 
tain  these  data. 

The  basic  construction  of  balsa,  aluminum  and  lead  shown 
in  Figure  D.l  is  the  same  for  all  the  models  tested.  The  balsa 
gives  the  desired  aerodynamic  shape.  The  airloads  are  trans¬ 
mitted  by  the  balsa  to  the  aluminum  spar,  whose  dimensions,  B 
and  H,  can  be  varied  to  control  the  torsional  rigidity  and  the 
bending  stiffness  of  the  model,  and  whose  chordwise  location 
determines  the  elastic -axis  position.  75ST  aluminum  alloy  was 
selected  for  the  spar  because  it  is  light  and  strong.  Also,  its 
yield  strength  is  close  to  its  ultimate  strength,  and  thus  a 
maximum  linear  range  is  available.  The  rectangular  cross- 
section  allows  for  ease  of  computation  during  the  design  and 
ease  of  fabrication  in  the  shop.  Two  pairs  of  strain  gages  are 
mounted  on  the  spar  for  use  in  the  vibration  and  flutter  tests; 
one  pair  along  the  elastic  axis  to  pick  up  bending,  and  the 
other  at  45°  to  the  elastic  to  pick  up  torsion.  Lead  weights, 
fore  and  aft  of  the  spar,  are  designed  so  as  to  give  the  de¬ 
sired  mass  and  inertial  properties  at  each  spanwise  station. 

The  lead  is  slit  before  gluing  so  that  its  contribution  to  the 
model's  stiffness  is  negligible. 

Because  the  thickness  ratio  of  the  wing  is  constant,  all 
dimensions  are  tapered  linearly  in  thickness  and  in  width. 


\l±f 

L  .  w 
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eALSA 


LEAD 


75  ST  ALUMINUM 
SPAR 


ELASTIC  AXIS 


B  ond  H  are  measured 
parallel  to  the  roof  chord 


FIGURE  D.l  A  TYPICAL  SUPERSONIC  FLUTTER  MODEL 


Therefore,  the  mass  varies  as  the  square  of  the  chord,  and  the 
stiffnesses  vary  as  the  fourth  power  of  the  chord, 


m  =  m 

o 

EI=  <EI>o  (s  )4 

GJ=  (CJ!0  )4 

'  O  / 


Eq. 

(D.l) 

Eq . 

(D.2) 

Eq. 

(0.3) 

WADC  TR  54-113  ,  Part  II 


147 


CONFIDENTIAL 


Model  nomenclature  is  explained  by  the  following  example: 

SW-2a-l 

SW  -  denotes  swept-wing  model. 

2  -  denotes  that  the  cross-sectional  dimensions  of  the 

spar  are  the  same  as  for  model  SW-2. 

a  -  denotes  the  first  basic  modification  of  model  SW-2 
(but  no  change  in  spar  dimensions). 

1  -  denotes  the  first  duplicate  model  of  SW-2a. 

D. 2  Design  Data 

Tables  D.l,  D.2  and  D.3  present  design  data  for  all  the 
models  tested.  Therefore,  any  model  can  be  duplicated  by  the  use 
of  Tables  D.l,  D.2  and  D.3  and  the  techniques  discussed  in  Refer¬ 
ences  14,  17,  and  19  for  designing  lead  weights  to  obtain  the 
desired  mass  and  inertial  characteristics  (see  Tables  D. 4,  D.5 
and  D.6).  The  spar  dimensions,  BQ  and  Ho,  were  measured  in  the 
streamwise  direction  for  all  the  planforms.  The  elastic  moduli 
of  the  balsa  were  obtained  experimentally  or  were  estimated 
using  the  methods  of  Reference  43.  Except  where  indicated  in 
Tables  D.4,  D.5  and  D.6,  the  grain  of  the  wood  was  parallel  to 
the  root  chord  for  the  straight  wings,  perpendicular  to  the 
spanwise  centerline  of  the  wing  for  the  swept  planforms,  and 
parallel  to  the  wing  centerline  for  the  delta  models. 

D.3  Mass  and  Stiffness  Data 

Tables  D.4,  D.5  and  D.6  give  mass  and  stiffness  data.  The 


WADC  TR  54-113,  Part  II  148 


CONFIDENTIAL 


CONFIDENTIAL 


mass  per  unit  length  of  the  wing  at  the  root  chord,  m  ,  was 
obtained  by  weighting  the  wing  and  using  the  relation, 

Mw  =  "Vjrf  I  +  M-?Z)  Eq.  (D.4) 

where 

is  the  total  mass  of  the  bare  wing. 

Equation  (D . 4 )  assumes  that  the  mass  varies  as  the  square  of 
the  wing  chord.  The  wing  stiffnesses  at  the  root,  (EI)0  and 
(GJ)q1  and  the  location  of  the  measured  elastic  axis  were  obtain¬ 
ed  by  static  tests,  assuming  that  the  stiffnesses  varied  as  the 
fourth  power  of  the  wing  chord.  The  location  of  the  measured 
elastic  axis  is  that  near  the  tip  of  the  wing,  and  it  varies 
slightly  along  the  span  because  of  sweep  and  root  effects. 
Therefore,  the  locus  of  shear  centers,  which  is  constant  along 
the  span,  was  assumed  to  correspond  to  and  was  calculated  as  the 
locus  of  centroids  of  the  bending  stiffnesses.  This  latter 

’  I 

parameter  is  probably  more  significant  than  the  location  of  the 
measured  elastic  axis  for  the  flutter  engineer.  Calculations 
were  made  to  obtain  the  section  center-of -gravity  location  and 
the  dimensionless  moment  of  inertia  in  pitch,  rK  ,  which  are 
constant  along  the  span.  Spot  tests  showed  that  the  experimen¬ 
tal  values  of  these  parameters  coincided  with  the  calculated 
values  within  the  experimental  accuracy. 

In  Table  D.6,  Models  De-2b,  De-3h  and  De-4a  were  duplicates 
of  previous  models  except  that  lead  was  added  to  lower  the 
second  bending  frequency  without  significantly  changing  the 
first  torsional  frequency  in  an  attempt  to  precipitate  flutter. 
The  location  of  the  lead  depended  on  the  mode  shape  for  the 
second  bending  frequency,  and  lead  was  placed  at  points  of  maxi¬ 
mum  amplitude  for  the  second-bending  mode  along  the  nodal  line 
for  first  torsion.  Therefore,  the  dimensionless  radius  of  gy- 
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ration  in  pitch,  ,  and  the  chordiwse  location  of  the  center 
of  gravity  are  not  quoted  for  these  models. 

D.4  Flutter  Data 


Tables  D.7,  D.8  and  I). 9  present  flutter  data  for  the 
straight,  swept  and  delta  planforms,  respectively.  Enough  data 
are  presented  so  that  other  flow  parameters  can  be  calculated. 
For  example, 


v  =  Ma 


Eq.  (D.5) 


=  /y  RT 


Eq.  (D.6) 


vu*  - 

/  IV /o  £ 
P  =  fO  RT 


Eq.  (D.7) 


Eq.  (D.8) 


RN  v 

IT  q) 

-8  ^/p 

.  _  I  3.059  x  10  °T  ^ 

v  ~  T"  t  m - 


i)  in  ft 


Eq.  (D.9) 


/o  in  slugs/ft3  Eq*  (D-10) 


T  in  C  absolute 


where 

v  is  the  free-stream  velocity 
M  is  the  free-stream  Mach  number 
a  is  the  speed  of  sound 

is  the  ratio  of  specific  heats 
R  is  the  gas  constant 

T  is  the  absolute  static  temperature  of  the  gas 
p  is  pressure 
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t~*-TIME~SEC.  (MEASURED  FROM  START  OF  TUNNEL) 

FIGURE  D .  2  TYPICAL  TIME  HISTORY  OF  FLOW  PARAMETERS  DURING  FLUTTER  TEST 
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-j-  is  the  Refolds  number  per  unit  length 
*0  is  kinematic  viscosity 

The  quantities,  v,  M  and yu.  ,  ar„e  tabulated  in  Tables  D.7,  D.8 
and  D. 9;  the  value  of  the  parameter,  m  ,  can  be  found  in  Tables 
D,4,  D. 5  and  D . 6 ;  for  dry  air, 

Y  -  1.4 

R  =  1718 

the  wing  semi-chord  at  the  root,  b  ,  is  5/12  foot  for  all  the 
models  tested.  The  Reynolds  number  was  on  the  order  of  7.5  x 
106  per  foot  of  reference  length  during  the  flutter  tests  (see 
Reference  13,  p.  51). 

All  legitimate  flutter  points  were  obtained  by  injecting 
the  model  into  a  stable  region  and  approaching  flutter  from  a 
high  Mach  number.  The  Mach  number,  velocity  and  wing  mass  ratio, 
yU.  ,  decreased  during  the  tests.  Figure. D, 2  shows  a  typical 
time  history  of  these  parameters  and  also  illustrates  an  inter¬ 
esting  case.  The  model  fluttered  during  injection  in  a  three- 
degree-of-freedom  flutter  mode  until  a  piece  of  the  trailing 
edge  of  the  aileron,  which  was  of  lead,  was  lost  (see  Figure 
D.5a).  Loss  of  the  trailing  edge  of  the  aileron  changed  the 
aileron  mass  and  inertial  characteristics,  and  the  model  became 

f 

stable  until  it  fluttered  in  a  bending-torsion  mode  at  a  lower 
Mach  number  when  a  piece  of  the  leading  edge  of  the  wing  was 
lost.  The  model  was  again  stable  until  it  fluttered  once  more 
at  a  lower  frequency  and  at  a  lower  Mach  number  in  a  bending- 
torsion  mode  causing  further  damage. 


ft" 

TTo: 


sec 
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(a)  Straight  Wing  with  Free  Aileron,  Cantilever  (Model  ST-le) 


FIGURE  D.3  SOME  HIGH-SPEED  MOVIES  AT  SUPERSONIC  FLUTTER 
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1  wing  with  Tip  Tank,  Free-to-Roll  (First  Stac< 

Model  ST-4c-3)  b 

'I CURE  D.3  (Continued)  SOME  HIGH-SPEED  MOVIES  AT  SUPERSON] 

FLUTTER 
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(c)  Straight  Wing  with  Tip  Tank,  Free-to-Roll  (Second  Stage, 

Model  ST-4c-3) 


FIGURE  D.3  (Continued)  SOME  HIGH-SPEED  MOVIES  AT  SUPERSONIC 


FLUTTER 
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« - -  FLOW  DIRECTION  ;  M>l 


(d)  Straight  Wing  with  Locked  Aileron,  Cantilever  (Model  ST-lf) 

FIGURE  D.3  (Continued)  SOME  HIGH-SPEED  MOVIES  AT  SUPERSONIC 

FLUTTER 
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(e)  Bare  Swept  Wing,  Cantilever  (Model  SW~bs  Injection  Flutter) 

FIGURE  D.3  (Continued)  SOME  HIGH-SPEED  MOVIES  AT  SUPERSONIC 

FLUTTER 

ir 
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FLOW  DIRECTION  ;  M>l 


(f)  Swept  Wing  with  Tip  Tank,  Cantilever  (Model  SW-6) 

FIGURE  D.3  (Continued)  SOME  HIGH-SPEED  MOVIES  AT  SUPERSONIC 

FLUTTER 
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(g)  Bare  Delta  Wing 
FIGURE  D.3  (Continued) 


Cantilever  (Model  De-4c) 

SOME  HIGH-SPEED  MOVIES  AT  SUPERSONIC 
FLUTTER 
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(a)  Straight  Wing  with  Free  Aileron,  Cantilever 


(b)  Straight  Wing  with  Tip  Tank,  Free-to-Roll  (First  Stage) 


(c)  Straight  Wing  with  Tip  Tank,  Free-to-Roll  (Second  Stage) 


FIGURE  D.4  WING-TIP  MOTION  AT  FI, UTTER 


WADC  TR  54-113,  Part  II  160 


CONFIDENTIAL 


CONFIDENTIAL 


INCHES 


RADIANS 


TIME-SECONDS 


(d)  Straight  Wing  with  Locked  Aileron,  Cantilever 
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RADIANS 


TIME-SECONDS 


(e)  Bare  Swept  Wing,  Cantilever  (Injection  Flutter) 


RADIANS 


TIME -SECONDS 


(f)  Swept  Wing  with  Tip  Tank,  Cantilever 
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FIGURE  D.4  (Continued),  WINC-TIP  MOTION  AT  FLUTTER 
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High-speed  movies  were  taken  during  most  of  the  flutter 
tests.  Besides  giving  evidence  of  model  damage,  the  high-speed 
photography  exhibited  clearly  the  exact  nature  of  the  flutter 
mode  shapes.  Figure  D.3  shows  some  high-speed  movies  of  typical 
models  at  flutter.  Timing  lines  could  be  obtained  on  tne  edge 
of  the  film  for  a  check  on  flutter  frequency,  independent  of  the 
oscillograph  records.  On  some  of  the  sequences  there  is  a 
record  of  Mach  number  at  the  bottom  of  each  frame,  which  was 
recorded  optically  by  a  system  of  mirrors  (see  Reference  14). 

Figure  D.4,  which  was  obtained  from  the  high-speed  movies, 
presents  quantitatively  the  time  history  of  the  motion  of  the 
wing  tip  for  some  models  at  flutter.  Measurements  were  taken 
from  the  high-speed  films  with  the  aid  of  a  microreader.  In 
general,  the  straight  wings  fluttered  in  bending-torsion;  the 
swept  wings  sometimes  had  a  small  amount  of  second  bending  in 

m 

the  flutter  mode,  while  it  is  difficult  to  say  what  modes  were 
involved  in  the  flutter  of  the  delta  wings  since  all  the  motion 
was  near  the  tip,  and  very  few  cycles  were  required  to  damage 
the  models. 

Flutter  of  the  bare  straight  wings  was  generally  very  vio¬ 
lent,  rarely  reaching  a  constant  amplitude  before  destruction; 
while  the  bare  swept  wings  usually  fluttered  at  constant  ampli¬ 
tude  for  a  number  of  cycles  before  the  models  were  badly  dam¬ 
aged.  (The  flutter  of  the  bare  swept  wing  which  appears  in 
Figure  D.3e  was  an  injection  flutter,  which  was  generally  very 
violent,  and,  in  this  respect,  is  not  typical  of  a  legitimate 
flutter  for  this  planform) .  The  delta  wings  which  fluttered 
were  usually  so  weak  that  few  cycles  were  required  for  damage. 
Flutter  of  models  with  tip  tanks  was  always  very  violent. 

A  particular  type  of  failure  at  flutter  characterized  each 
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planforra.  Figure  D.5  Shows  some  typical  models  after  flutter. 
For  the  bare  straight  wings  an  outboard  piece  of  the  wing  at  the 
leading  edge  was  lost;  for  the  bare  swept  wings  an  entire  sec¬ 
tion  of  the  model  just  inboard  of  the  wing  tip  gave  way;  while 
for  the  delta  wings  a  large  section  of  the  wing  at  the  trailing 
edge  was  destroyed.  Wings  with  ailerons  often  lost  a  piece  of 
the  aileron  trailing  edge  at  flutter,  and  those  with  tip  tanks 
were  usually  totally  destroyed.  After  datnage  the  models  often 
became  stable,  sometimes  fluttering  again  at  a  lower  Mach  number 
before  retraction. 

D . 5  Tip -Tank  Parameters 

Tables  D.10  and  D.ll  give  pertinent  tip-tank  parameters 
for  the  straight  and  swept  wings,  respectively.  The  total  mass 
of  the  wing  can  be  obtained  by  use  of  equation  (D.4).  Except 
for  the  location  of  the  center  of  gravity  of  the  tip  tank,  the 
models  were  designed  to  have  constant  dimensionless  properties 
including  the  dimensionless  tip-tank  moment  of  inertia  in  pitch 
about  the  wing  elastic  axis.  Some  variation  in  this  latter  par¬ 
ameter  occurred,  however,  because  the  tip.  tanks  were  designed 
for  a  moment  of  inertia  in  pitch  about  an  estimated  elastic- 
axis  position,  for  the  models  were  not  statically  tested  until 
after  the  design  stage.  Later,  after  some  research,  it  was  de¬ 
cided  that  the  locus  of  shear  centers  was  more  significant  than 
the  measured  elastic -axis  location.  The  center  of  gravity  of 
the  tip  tank  was  changed  by  locating  lead  weights  according  to 
the  design  procedure  outlined  in  Reference  (19)  and  it  was 
checked  experimentally  along  with  the  weight  and  the  mass  moment 
of  inertia  in  pitch.  Therefore,  the  geometrical  location  of  the 
tip  tank  with  respect  to  the  wing  was  constant  for  each  model, 
even  though  the  static  unbalance  of  the  tip  tank  varied.  In 
Table  D.10  the  parameter,  "5^,  is  tabulated  for  convenience  in 


WADC  TR  54-113,  Part  II  164 


CONFIDENTIAL 


CONFIDENTIAL 


using  the  theory. 

For  the  free-to-roll  models  a  2  1/2-degree  canted  hinge  was 
used  to  aerodynamically  stabilize  the  wing  (see  Figure  D.6). 

The  models  were  mounted  on  a  simulated  fuselage  which  rolled 
with  .the  model.  Geometric  parameters  for  the  fuselage  were 
chosen  so  as  to  be  typical  of  supersonic  aircraft,  and  the  roll 
axis  corresponded  to  the  centerline  of  the  fuselage.  The  value 
of  the  mass  moment  of  inertia  in  roll  of  the  root  support  was 
chosen  so  that  its  value  corresponded  to  a  fuselage  radius  of 
gyration  in  roll  typical  of  present-day  aircraft.  Its  value 
was  also  checked  experimentally. 

D.6  Aileron  Parameters 


Table  D.12  gives  pertinent  aileron  parameters  for  the  models 
tested.  In  all  cases  the  ailerons  had  sealed  gaps.  Again,  lead 
weights  were  used  to  obtain  the  desired  mass  and  inertial 
characteristics.  The  aileron  was  attached  to  the  wing  by  flex¬ 
ures  whose  size  determined  the  aileron  frequency;  no  mechanism 
was  used  to  adjust  the  aileron  frequency  from  outside  the  tunnel. 
In  the  design  of  the  lead  weights,  allowance  was  made  for  the 
weight  of  the  flexures  so  that  the  dimensionless  mass  and  in¬ 
ertial  parameters  for  the  total  wing  and  aileron  remained  con¬ 
stant  at  all  spanwise  stations.  Less  accuracy  is  quoted  for  the 
elevon  parameters  of  the  delta  wing  (Table  12. b)  since  the  de¬ 
sign  was  complicated  by  the  constant  chord  of  the  elevon.  For 
the  straight  and  swept  wings,  the  chord  and  the  thickness  of  the 
ailerons  tapered  linearly  to  the  tip  so  that  the  lead  weights 
had  to  be  designed  at  only  one  spanwise  station,  for  linear 
taper  to  the  tip  insured  the  mass  variation  of  equation  (D.l) 
and  constant  dimensionless  parameters.  For  the  delta  wings  with 
elevons,  the  advantages  of  linear  taper  are  lost,  and  the  lead- 
weight  design  must  be  made  at  more  than  one  spanwise  station. 
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D . 7  Vibration  Data 


Tables  D.14,  D.15  and  D.16  present  nodal  lines  and  struc¬ 
tural-damping  coefficients  for  the  models  tested.  Only  elastic 
frequencies  are  presented;  for  the  free-to-roll  models ,  the 
first  natural  frequency  is  a  rigid-body  rotation  about  the  roll 
axis.  The  nodal  lines  are  approximate,  and  no  attempt  was  made 
to  obtain  mode  shapes.  The  nodal  lines  were  obtained  optically 
with  salt  while  the  wing  was  vibrating  at  a  natural  frequency. 
Nodal  areas  as  well  as  nodal  lines  appeared  because  of  the  limi¬ 
tation  in  applied  shaking  force.  All  the  frequencies  quoted  are 
coupled,  including  those  for  the  models  with  ailerons,  since  the 
ailerons  were  not  clamped  with  respect  to  the  wings  during  the 
shake  tests.  The  first  two  uncoupled  frequencies  of  all  the 
bare  straight  wings  were  calculated  from  the  coupled  frequencies 
and  the  model  parameters.  Only  the  first  two  natural  frequen¬ 
cies  were  considered.  The  results  show  that  the  first  un¬ 
coupled  frequency  is  about  1$  higher  than  the  first  coupled  fre¬ 
quency  and  that  the  second  uncoupled  frequency  is  about  1  1/2$ 
lower  than  the  second  coupled  frequency. 

D.8  Influence  Coefficients 


Tables  D.17  and  D.18  and  Figures  D.7  and  D.8  present  in¬ 
fluence-coefficient  matrices  with  locations  for  some  of  the 
models  tested.  The  influence  coefficients  were  measured  with 
linear  variable  differential  transformers,  and  their  accuracy 
is  estimated  at  -  0.0005  inch/#  or  less. 

D.9  XF-92A  Airplane 

Table  D.19  gives  some  available  experimental  data  for  the 

XF-92A  airplane  (References  44-46),  since  this  airplane  was  used 
as  a  guide  in  designing  the  delta-wing  models. 
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TABLE  D.l 

DESIGN 

DATA  FOR  STRAIGHT  WINGS 

- 

Bare  Wings 

Model 

£ 

Bo 

H0flnflr  - 

eba 

gba 

spar 

.  %  chord 

spar 

inches 

spar 

inches 

psi 

psi 

ST-1 

41.8 

2.721 

0.171 

0.0045* 

8710* 

22,300* 

ST- Id 

37.0 

2.721 

0.171 

0.0052* 

9590 

29,600 

ST-ld-1 

37.0 

2.721 

0.171 

0.0052* 

9590* 

29,600* 

ST- 2 

42.0 

0.664 

0.251 

0.0045* 

8710* 

22,300* 

ST-3 

42.0 

0.257 

0.300 

0.0045* 

8710* 

22,300* 

ST-4 

42.0 

4.214 

0.160 

0.0043 

6380 

16,800 

ST-4-1 

42.0 

4.214 

0.160 

0.0043* 

6380* 

16,800* 

ST -5 

42.0 

3.141 

0.175 

0.0041 

'  6750 

21,800 

ST-5-1 

42.0 

3.141 

0,175 

0.0041* 

6750* 

21,800* 

ST -6 

42.0 

1.049 

0.229 

0.0048 

9510 

26,900 

ST-7 

42.0 

0.618 

0.252  . 

0.0049 

10,700 

27,400 

ST-7-1 

42.0 

0.618 

0.252 

0.0045* 

8710* 

22,300* 

ST-7-2 

42.0 

0.618 

0.252 

0.0045* 

8710* 

22,300* 

ST-7-3 

42.0 

0.618 

0.252 

0.0045* 

3710* 

22,300* 

ST-8 

42.0 

0.422 

0.266 

0.0039 

5650 

25,200 

ST-12 

41.8 

3.00** 

0.200** 

0.0045 

6720 

23,100  - 

Wings  with  Tip  Tanks 

Model 

£ 

Bo 

Ho 

Sa 

sa 

spar 
% chord 

spar 

inches 

spar 

inches 

psi 

uBBkat 

ST- la 

msm 

0.171 

0.0045* 

8710* 

22,300* 

ST-lc 

1 

0.171 

0.0030 

7800 

17,400 

ST-4a 

ESI  ! 

0.160 

0.0039* 

7880 

19,800 

ST-4b 

0.160 

0.0039* 

7880 

19,800 

ST-4c , -1 

1 

4 .214 

0.160 

0.0045* 

8710* 

22,300* 

ST-4c-2 

4 .214 

0,160 

0.0045* 

8710* 

22,300* 

ST-4c-3 

Hi 

4 .214 

0.160 

0.0045* 

8710* 

22,300* 

Wings  with  Ailerons 

Model 

£ 

Bo 

Ho 

PBA  , 

eba 

spar 
%  chord 

spar 

inches 

spar 

inches 

it-/  in^ 

psi 

ysH 

K  5 

SEP 

wmm 

0.0045* 

8710* 

22,300* 

1 

1 

0.0032 

11 ,800 

20,400 

ESI 

HI 

2.721 

0.0032 

11  ,800 

20,400 

*Assumed  Balsa  Properties 

**Built-up  rectangular  spar;  24ST  Alloy,  skin  thickness  at  root  =  0.023";  Balsa  Core; 
all  dimensions  tapered  linearly  to  tip 
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TAELE  D . 2 


DESIGN  DATA  FOR  SWEPT  WINGS 


spar 
%  chord 


Bare  Wings 


spar 

inches 


Pba 

*/in3 


0.0160** 

0.0160** 

0.0040* 

0.0040* 

0.0040 

0.0038 

0.0052 

0.0055 

0.0053 

0.005? 

0.0048 

0.0060 

0.0060* 


SW-2a 

SW-2a-l 

SW-2b 

SW-2b-l 

SW-2b-2 

SW-3 

SW-3-1 

SW-3c 

SW-3d 

SW-4 

SW-5 

SW-8 

SW-8-1 


spar 
%  chord 


1.355 

1.355 

1.355 

1.355 

1.355 

0.783 

0.783 

0.783 

0.783 

0.661 

0.410 

0.668 

0.668 


spar 

inches 


0.253 

0.253 

0.253 

0.253 

0.253 

0.292 

0.292 

0.292 

0.292 

0.246 

0.306 

0.334 

0.334 


Wings  with  Tip  Tanks 


spar 

inches 


/in~ 


988,500* 

988,500* 

5700* 

5700* 


1 


-3a 
-6 
SW-7 
SW-7a 


0.783 

1.118 

O.896 

O.896 


o'.  292 

0.334 

0.334 

'0.334 


Wings  with  Ailerons 


32,000 
33,600 
20,000 
39,4 


Model 

£ 

spar 
%  chord 

Bo 

spar 

inches 

Ho 

spar 

inches 

PBA  , 

HI  in3 

^A 

psi 

CBA 

psi 

SW-3b 

41.9 

0.783 

0.292 

0.0057 

10,800 

32,600 

Note:  E  for  wood  measured  for  bending  of  the  elastic  axis 

G  for  wood  measured  for  twisting  about  the  elastic  axis 
•Estimated  balsa  properties  **Estimated  pine  properties 
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•TABLE  D.3 


DESIGN  DATA  FOR  DELTA  WINGS 


Bare  Wings 


De-l 

50.0 

- — - - 

5.00 

0.0798 

0.0048* 

730.000* 

De-la 

50.0 

5.00 

0.0798 

0.0048* 

730,000* 

De-la-1 

50.0 

5.00 

0.0798 

0.0048* 

730,000* 

De-la-2 

50.0 

5.00 

O.O798 

0.0048* 

730,000* 

De-2 

40.0 

3.00 

0.0500 

0,0035 

530,000* 

De-2a-l 

40.0 

3.00 

0.0500 

0.0030* 

460,000* 

De-2b 

40.0 

3.00 

0.0500 

0.0030* 

460,000* 

De-2c 

40.0 

3.00 

0 .0500" 

0.0030* 

460,000  * 

De-2d 

40.0 

3.00 

0 .0500 

0.0030* 

460,000* 

De-2d-l 

40.0 

3.00 

0 .0500 

0.0030*- 

460,000* 

De-2d-2 

40.0 

3.00 

0 .0500 

0.0030* 

460,000* 

De-3-1 

40.0 

2.40 

0 .0500 

0.0045* 

690,000* 

De-3a 

40.0 

2.40 

0 .0500 

0.0046 

700,000  * 

De-3b 

40.0 

2.40 

0 .0500 

0.0032 

490,000* 

De-3c 

40.0 

2.40 

0 .0500 

0 .0027 

410,000  * 

De-3d 

26.0 

2.40 

0 .0500 

0.0030 

460,000  * 

De-3e 

40.0 

2.40 

0  .0500 

0.0027* 

410,000  » 

De-3f 

40.0 

2.40 

0 .0500 

0  .0032*' 

490,000  » 

De-3f-l 

40.0 

2.40 

0  .0500 

0 .0032* 

490,000  * 

De-3g 

40.0 

2.40 

0  .0500 

0 .0046* 

700,000  * 

De-3h 

40.0 

2.40 

0  .0500 

0  .0045* 

690,000  * 

De-3i 

40 .0 

2.40 

0  .0500 

0 .0045* 

690,000  * 

De-3j 

40 .0 

2.40 

0  .0500 

0 .0027* 

410,000  * 

De-4 

-- 

-- 

-- 

0  .0032 

490,000  * 

De-4a 

-- 

-- 

-- 

0 .0032* 

490,000  * 

De-4b 

-- 

-- 

-- 

0 .0032* 

490,000  * 

De-4c 

-- 

-- 

0  .0032* 

490,000  * 

De-4c-l 

-- 

-- 

-- 

0  .0032* 

490,000  * 

De-5 

-- 

-- 

-- 

0 .0026 

2300 

De-6 

-- 

-- 

-- 

0  .0026* 

2300* 

Wings  with  Elevons 


25,500* 
25,500* 
25,500* 
25,500* 
17,200* 
14 , 100* 
14 , 100* 
14,100* 
14 , 100* 
14 , 100* 
14 , 100* 
23,500* 
23,600* 
15,500* 


Model 

spar 
%  chord 

B 

0 

spar 

inches 

H 

°spar 

inches 

Pba 

Win3 

^BA. 

psi 

De-2e 

40 .0 

3.00 

0.0500 

0.0044 

670,000* 

♦Estimated  balsa  properties 
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TABLE 

D  .4 

MASS  AND  STIFFNESS  DATA 

FOR  STRAIGHT  WINGS 

Eare  Wings 

Model 

ilu§s/in 

sec Cion 
c.g. 

%  chord 

r  2 
a.42c 

(ED0 

x  lO'l 
a  -in* 

(GJ)0 

x  10'f 
*  -in^ 

calc . 
locus  of 
shear 
centers 
%  chord 

measured 
elastic 
axis 
%  chord 

Remarks 

0.00406 

50.0 

0.250 

0.0119 

0.0249 

42.4 

44.5 

bare  wing 

0.00464 

50.0 

0.250 

0.0122 

0.0226 

37.5 

39-3 

bare  wing, 

Bpar  forward 

ST-ld-l 

0.00464 

50.0 

0.250 

0.0112 

0.0239 

37.5 

42.8 

rebuilt  ST- Id 

ST-2 

0.00378 

50.0 

0.250 

0.0095 

0.0165 

42.6 

48.5 

bare  wing 

ST-3 

0.00447 

50.0 

0.250 

0.0067 

0.0103 

42.9 

48.1 

bare  wing 

ST-4 

0.00387 

50.0 

0.250 

0.0150 

0.0290 

42.4 

46.5 

bare  wing 

ST-4-1 

0.00401 

50.0 

0.250 

0.0159 

0.0307 

42.4 

44.6 

rebuilt  ST-4 

ST- 5 

0.00457 

50.0 

0.250 

0.0147 

0.0297 

42.2 

43.7 

bare  wing 

ST-5-1 

0.00457 

50.0 

0.250 

0.0146 

0.0299 

42.2 

46.2 

rebuilt  ST-5 

ST-6 

0.00433 

50.0 

0.250 

0.0109 

0.0187 

42.3 

51.3 

bare  wing 

ST-7 

0.00427 

50.0 

0.250 

0.0098 

0.0176 

42.4 

48.0 

bare  wing 

ST-7-1 

0.00427 

50.0 

0.250 

0.0090 

0.0171 

42.7 

46.5 

rebuilt  ST-7 

ST-7-2 

0.00427 

50.0 

0.250 

0.0088 

0.0155 

42.7 

43.3 

rebuilt  ST-7-1 

ST-7-3 

0.00428 

50.0 

0.250 

•0.0049 

0.0157 

42.7 

47.8 

rebuilt  ST-7-2 

ST-8 

0.00418 

50.0 

0.250 

0 .0079 

0.0136 

42.3 

47.5 

bare  wing 

ST- 12 

0.00229 

50.0 

0.250 

0.0152 

0.0289 

42.0 

_ 

low-density 
bare  wing 

Wings  with  Tip  Tanks 

Model 

s 

mo 

lugs /in 

section 

c.g. 

%  chord 

r  2 
a.42c 

(“>o 

x  10’f 
«-in“ 

(GJ)0 
x  10*£ 

8  -inc 

calc . 
locus  of 
shear 
centers 
%  chord 

measured 
elastic 
axis 
%  chord 

Remarks 

ST-la 

0.00413 

50.0 

0.250 

0.0116 

0.0232 

42.5 

42.8 

bare  wing  same 
as  ST-1 

ST-lc 

0.00447 

50.0 

0.250 

0.0120 

0.0227 

37.4 

39.1 

bare  wing  same 
as  ST-ld 

ST-4a 

0.00406 

50.0 

0.250 

0.0150 

0.0280 

42.2 

47.5 

bare  wing  same 
as  ST-4 

ST-4b 

0.00406 

50.0 

0.250 

0.0150 

0.0280 

42.2 

47.5 

bare  wing  same 
as  ST-4 

ST-4c 

0.00387 

50.0 

0.250 

0.0147 

0.0303 

42.4 

49.5 

bare  wing  same 
as  ST-4 

ST-4c-1 

0.00387 

50.0 

0.250 

0.0199 

0.0364 

42.4 

42.4 

rebuilt  ST-4c 

ST-4c-2 

0.00387 

50.0 

0.250 

0.-0209 

0.0352 

42.4 

39-9 

rebuilt  ST-4c-l 

ST-4c-3 

0.00387 

50.0 

0.250 

0.0194 

0.0348 

42.4 

41.0 

rebuilt  ST-4c-2 

Wings  with  Ailerons 

Model 

f 

m 

0 

lugs /in' 

section 

c-g- 

%  chord 

r4c 

(EI)0 

x  10 
e  -ind 

<GJ)C 

x  10'f 

8  -  ir. 

calc . 
locus  of 
shear 
centers 
%  chord 

measured 
elastic 
axis 
%  chord 

Remarks 

ST- lb 

0.00407 

50.0 

0.250 

0.0128 

0.0191 

42,4 

47.0 

same  as  ST-1 
but  with  aileron 

ST-  le 

0.00388 

50.0 

0.250 

0.0123 

0.0220 

42.4 

48.0 

same  as  ST-1 
but  with  aileron 

ST- If 

0.00388 

50,0 

0.250 

0.0123 

0.0219 

42.4 

41.7 

same  as  ST-1 
but  with  aileron 
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TABLE  D.5 

MASS 

AND  STIFFNESS  DATA  FOR  SWEPT  WINGS 

Bare  Wings 

Model 

section 

c.g. 

%  chord 

r  2 

a  .42c 

(EI)0 
x  10'6 
#  -  in2 

(GJ)0 

x  10'6 
*  -  in2 

calc . 
locus  of 
shear 
centers 
%  chord 

measured 
elastic 
axis 
%  chord 

Remarks 

SW-2a 

0.00468 

50.0 

.306 

0.0153 

0.0328 

4  5'- 9 

54.2 

made  of  pine,  grain 
perpendicular  to  the 
t  of  wing 

SW-2a-l 

0.00468 

50.0 

.306 

0.015^ 

0.0305 

45.9 

50.0 

repaired  SW-2a 

SW-2b 

0.00452 

50.0 

■  375 

0.0148 

0.0227 

35-4 

41.2 

balsa  grain  perpen¬ 
dicular  to  center- 
line  of  wing 

SW-2b-l  ' 

0.00452 

50.0 

•  375 

0.0211 

35.4 

41.2 

repaired  SW-2b 

SW-2b-2 

0.00455 

50.0 

■  375 

0.0219 

35-4 

47.0 

repaired  SW-2b-l 

SW-3 

0.00401 

50.0 

.375 

0.0135 

0.0163 

42.0 

47.0 

bare  wing 

SW-3-1 

0.00413 

50.0 

■  375 

0.0133 

0.0199 

42.0 

50.8 

rebuilt  SW-3 

SW-3c 

0.00392 

50. 0 

.375 

0.0127 

0.0173 

37.0 

47.5 

same  as  SW-3  but 
elastic  axis  for¬ 
ward 

SW-3d 

0.00398 

50.0 

.375 

0.0140 

0.0191 

m 

55-1 

same  as«8W-3  but 
elastic  axis  aft  , 

SW-4 

0.00406 

50.0 

.3^5 

0.0074 

0.0124 

50.0 

bare  wing 

SW-5 

0.00391 

50. 0‘ 

.375 

0.0078 

0.0095 

1 

52.8 

bare  wing 

SW-8 

0.00412 

50.0 

.375 

0.0173 

0.0199 

46.0 

bare  wing 

SW-8-1 

0.00412 

50.0 

■  375 

0.0173 

0.0199 

46.0 

repaired  SW-8  wing 

Wings  with  Tip  Tanks 

» 

Model 

m 

0 

slugs /in 

section 

c.g. 

% chord 

r  2 
a.42c 

(SUo 

x  10"6 
-in2 

(GJ}0 

x  10'6 

2 

*  -in 

calc . 
locu^  of 
shear 
centers 
%  chord 

measured 
elastic 
axis 
%  chord 

Remarks 

SW-3a 

0.004 06 

50.0 

•  375 

0.0130 

0.0180 

42T6 

52.8 

same  as  SW-3  but 
with  tip  tank 

SW-6 

0.00409 

50.0 

■  375 

0.0274 

0.0343 

42.0 

47.0 

cantilever  tip- 
tank  wing 

SW-7 

0.00391 

50.0 

•  375 

0.0221 

0.0245 

42.0 

48.8 

cantilever  tip- 
tank  wing 

SW-7a 

j 

0.00406 

50.0 

•  375 

0.0213 

0.0288 

42.0 

51.4 

free  to  roll,  tip- 
tank  wing 

Wings  with  Ailerons 

Model 

mo 

slugs/in 

section 

c.g. 

%chord 

r  2 
a.42c 

(EI)0 

x  10"6 
*  -  in2 

(GJ)0 

x  10’5 
0  -  in2 

calc, 
locus  of 
shear 
centers 
%  chord 

measured 
elastic 
axis 
%  chord 

Remarks 

SW-3b 

0.00^26 

50.0 

•  375 

0.0140 

0.0203 

42.0 

51.3 

wing  with  aileron 
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TABLE  D.6 

L 

MASS  AMD  STIFFNESS 

DATA  FOR  DF.LTA  WINGS 

Bare  Wings 

Model 

m 

section 

r  2 

calc . 

Remarks 

T 

slugs/ In 

c.g. 

%  chord 

“,42c 

locus  of 
shear  centers 

%  chord 

De-l 

0.060 

0.00427 

50.0 

0.250 

50.0 

bare  wing 

De-la 

0.040 

0.00398 

50.0 

0.250 

50.0 

same  as  De-l  but 
with  T  reduced 

De- la-1 

0.040 

0.00409 

50.0 

‘0.250 

50.0 

rebuilt  De-la 

De-la-? 

0.040 

0.00409 

50.0 

0.250 

50.0 

rebuilt  De-la-1 

De-2 

0.060 

0.00462 

50.0 

0.250 

49.9 

bare  wing 

De- 2a- 1 

0.060 

0.00459 

50.0 

0.250 

49.9 

bare  wing 

De-2b 

0.060 

0.0050? 

49.9 

same  as  De-2a-l  with 
lead  added 

De-2c 

0.050 

0.00445 

50.0 

0.200 

49  .7 

same  as  De-2a-l  but 
with  T  reduced 

De-2d 

0.040 

0.00430 

50.0 

0.250 

49.4 

same  as  De-2a-l  but 
with  T  reduced 

De-2d- 1 

0.040 

0.00470 

50.0 

0.250 

49.5 

rebuilt  De-2d 

De-@d-2 

0.040 

0.00470 

50.0 

0 . 260 

49.5 

rebuilt  De-2d-l 

De-3-1 

0.060 

0.004  52 

50. 0 

0.250 

49.9 

bare  wing 

De-3a 

0.060 

0.00455 

45.0 

0.250 

49.9 

same  as  De-3-1  but 
with  c.  g.  forward 

De-3b 

0.050 

0.00469 

55.0 

0.250 

49.9 

same  as  De-3-1  but 

tt 

with  c.  g.  aft 

De-3c 

0.060 

0.00276 

50.0 

O'.  250 

49-9 

same  as  De-3-1  but 
less  dense 

De-3d 

0.060 

0.00420 

50.0 

0.250 

49.7 

same  as  De-3-1  but 
with  spar  forward 

De-3e 

0.030 

0.00234 

50. C 

0.250 

49.0 

same  as  De-3c  but 
with  T  reduced 

De-3f 

0.040 

0.00441 

55.0 

0.250 

49.6 

same  as  De-3b  but 
with  T  reduced 

De-3f-l 

0.040 

0.00448 

55.0 

0.240 

49.6 

rebuilt  De-3f 

De-3g 

0.050 

0.00441 

4S.0 

0.240 

49.9 

same  as  De-3a  but 
with  r  reduced 

De-3h 

0.060 

0.00499 

49.9 

same  as  De-3-1  but 
with  lead  added 

De-3i 

0.045 

0.00430 

50.0 

0.2)0 

4q  .8 

same  as  De-3-1  but: 
with  r  reduced 

De-3j 

0.040 

0.00231 

60.0 

0.240 

49  .6 

same  as  De-Jc  but 
with  T  reduced 

De-4 

0.060 

0.00336 

50.0 

0.240 

50.0 

no  spar 

De-4a 

0.060 

0.00381 

50.0 

no  spar,  same  as  De-^ 
with  lead  added 

De-4b 

0.050 

0.00321 

50.0 

0.250 

50.0 

no  spar,  same  as  De-^ 
but  with  T  reduced 

De-^c 

0.040 

0.00307' 

50.0 

0.250 

50.0 

no  spar,  same  as  De-4 
but  with  T  reduced 

De-4c-l 

0.040 

0.00302 

50.0 

0.250 

50.0 

rebuilt  De-4c 

De-5 

0.030 

0.00438 

60.0 

0.300 

50.0 

no  spar,  b.?lsa  grtfin 
parallel  to  root 

De-6 

0.030 

0.00583 

60.0 

0.300 

50.0 

same  as  De-5;  0.00*;  in. 
brass  shim  in  mid- 

plane 

Wings  with  Elevons 

Model 

mo 

section 

? 

calc-. 

Remarks 

T 

slugs/ in 

c.g. 

%  chord 

r° . 4  2c 

locus  of 
shear  centers 

chord 

De-2e 

0.040 

0.00413 

•  >0 . 0 

0.250 

49.7 

aileron  wing 
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EXPERIMENTAL  FLUTTER  DATA  FOR  STRAIGHT  WINCS 


Bate  Wings 


ST- 1 

ST- Id 

ST-ld-1 

ST- 2  • 

ST- 3 

ST-4 

ST-4-1 

ST-5 

ST-5-1 

ST-6 

ST-7 

ST-7-1* 

ST-7-2 » 

ST-7-3 
ST-8 
ST- 12 


6.20 
6.84 
53.2  6.20 
77.0  5.42  5-42 

55.6  5-79  3.25 
270  47.2  4ig3  2.7fi 

370  64.9  6.28  3-45 

36o  54.8  5.30  3.44 

530  76.2  7.46  4.15 

620  81.7  8.05  5-48 

630  83.9  5.76  5.61 

650  72.4  7.73  5.12 

360  27,4  4.30  2.28 


53.0  0.308 
40.5  0.299 
42.0  0.329 
10.0  0.345 
0.382 
0.349 
0.338 
0.332 
0.351 
40.9  0.363 
23.3  0.368 
13.0  0.346 


1.18 

legitimate  flutter 

1.20 

legitimate  flutter 

1.22 

legitimate  flutter 

1.36 

injection  flutter 

1.51 

destroyed  during  Injection 

1.30 

legitimate  flutter 

1.35 

legitimate  flutter 

1.25 

legitimate  flutter 

1.30 

legitimate  flutter 

1.28 

legitimate  flutter 

1.26 

destroyed  by  starting  shock 

1.31 

Injection  flutter 

1.35 

injection  flutter 

1.34 

legitimate  flutter 

1.40 

destroyed  during  injection 

1.19 

legitimate  flutter 

Wings  with  Tip  Tanks 


21.4 

1630 

75.6 

29.1 

4.071 

153-0 

0.3081 

18.1 

1600 

89.5 

33-8 

4.431 

138l.o 

0.3191 

26.3 

1360 

53.2 

19.8 

? 

177  *.3 

0.3231 

29.4 

1590 

79.1 

20.7 

3.431 

177  .u 

0.3231 

- 

1390 

52.6 

- 

3.031 

l75l.o 

0.3431 

33.9 

1290 

59.5 

14.5 

2.821 

l"5^o 

1,2 

0.3*3 

30.0 

1260 

57.5 

I6.0 

2.751 

1,2 
1T5  .0 

1,2 

0.343 

32.7 

1280 

56.8 

15.0 

2.791 

175  l.O 

0.3651 

Model 

r— i 

Hf 

mf 

vf 

cps 

ft/sec 

Wings  with  Ailerons 


ST- lb* 

1.80 

q 

d 

0 

CU 

1620 

84.0 

3.09 

3.19 

194 .0 

ST-  le 

1 .72 

87.9 

1550 

69.9 

6.74 

4.20 

14 1  .0 

ST  If 

1.65 

87.5 

V-T' 

O 

O 

64  .7 

6.55 

_ i 

3.89 

147. 2 

cantilever;  In¬ 
jection  flutter 

cantilever;  In¬ 
jection  flutter 

cantilever;  legiti¬ 
mate  flutter 

cantilever;  in¬ 
jection  flutter 

free-to-roll;  against 
roll  stop  during 
flutter 

free-to-roll; 
legitimate  flutter; 
flutter  mode  shape 
changed  from  pri¬ 
marily  torsion  to 
primarily  bending 
and  roll  at  lower 
Mach  number 


aileron  free;  In¬ 
jection  flutter 

aileron  free;  legi¬ 
timate  flutter;  three 
degrees  of  freedom 
in  .flutter  mode 

aileron  locked; 
legitimate  flutter 


*  Inlection  flutter 


1  Based  on  cantilever  frequencies  of  bare  wing  (see  Table  T.14  for  further 

2  Data  for  ST-4c  used  since  these  are  rebuilt  wings 

3  Coupled  aileron  frequency 


Vibration  Data) 
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TABLE  D.7a 

EXPERIMENTAL 

FLUTTER  DATA  FOR  STRAIGHT  WINGS  (Cont.) 

Bare  Wings 

Model 

v£ 

Vf 

40ns Ho  f/U-i\ 

■MBBM 

a< 

1  ~{ 

9.92 

5-53 

0.268 

7  $5.^' 

10.94 

6.43 

0.258 

ST-ld-1 

9.92 

6. -03 

0.228 

ST-2* 

8.67 

8.67 

0.215 

ST-4 

9.26 

5.20 

0.271 

ST-4-1 

7.89 

4 . 42 

0.251 

ST-5 

10.05 

5.52 

0.261 

ST-5-1 

9.28 

5.  ^0 

0.245 

ST-6 

n.94 

6.64 

0.281 

ST-7-1* 

12.88 

8.77 

0.241 

ST-7-2* 

9.22 

8.98 

0.243 

ST-7-3 

12.37 

8. 19 

0.250 

ST- 12 

6.88 

3-65 

0.258 

Wings  with  Tip  Tanks 

ST-la* 

46.6 

-  - 

ST-lc* 

54.1 

7.09'1^ 

-- 

ST-4a 

31.7 

4.69(1) 

-- 

ST-4b* 

33-1 

5.49^^ 

-- 

ST-4c 

-- 

4.83^ 

ST-4c-l 

23.2 

4.51^^ 

-- 

ST-4c-2 

25.6 

4.4o^ 

-- 

ST-4c-3 

24.0 

4.46^ 

Wings  with  Ailerons 

ST-lb* 

4.94 

5.10 

v  n 

ST-le 

10.8 

6.72 

— 

ST- If 

10.5 

6.22 

0.264 

*  IniectJ.on  Flutter 
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TABLE  D.8 

EXPERIMENTAL  FLUTTER  DATA  FOR  SWEPT  WINGS 


3a.’e  Wings 


Model 

Hf 

“f 

m 

vf 

0 

% 

% 

Remarks 

cps 

■ 

^f 

gg 

al 

cps 

““1 

““1 

SW-2a 

-- 

-- 

97-63 

.. 

225 

0.169 

0.654 

no  flutter 

SW-2a-l 

1.79- 

1.32 

-- 

-- 

74-50 

-- 

217 

0.162 

0.608 

no  flutter 

SW-2b 

1.75- 

1.34 

-■ 

-- 

86-56 

-- 

140 

0.224 

0.743 

no  flutter 

SW-2b-l 

1.73- 

1.74 

-- 

-- 

82-54 

-- 

-- 

137 

0.228 

O.796 

no  flutter 

SW-2b-2 

1.48- 

1.28 

-- 

-- 

76-53 

-- 

— 

140 

0.222 

0.793 

no  flutter 

SW-3 

1.41 

111 

1340 

51.9 

4.61 

3.88 

132 

0.254 

0.902 

legitimate  flutter 

SW-3-1  + 

1.26 

100 

1240 

53.1 

4  .74 

3.24 

146 

O.234 

0.822 

retraction  flutter 

SW-3c 

1.48 

110 

1390 

63.4 

4.83 

3-93 

135 

0.255 

0.844 

"legitimate  flutter 

SW-3d 

1.25 

116 

1230 

50.5 

4  .05 

3.43 

137 

0.258 

0.927 

legitimate  flutter 

SW-4  * 

I.85 

90.4 

1620 

84.3 

6.85 

5.U 

121 

0.225 

0.807 

injection  flutter 

SW-5* 

1.87 

77.8 

1630 

73.2 

8.00 

5.51 

113 

0.247 

0.870 

injection  flutter 

SW-8 

1.62- 

1.26 

“* 

-- 

76-53 

150 

0.252 

0.853 

no  flutter,  first  test 

SW-8* 

1.62 

116 

1470 

58.9 

4.84 

3.74 

150 

0.252 

0.853 

rerun  SW-8  at  lower  p; 
Injection  flutter 

sw-3-1 

1.42 

120 

1330 

49.5 

4.23 

3.48 

146 

0.255 

0.856 

legitimate  flutter 

Wings  with  Tip  Tanks 


Model 

Kf 

'“f 

vf 

n 

vf 

1_ 

vf 

m 

% 

cps 

ft/sec 

■ 

^f 

mg 

G1 

cps 

CD 

Qi 

““1 

SW-3a» 

1.92 

35.5 

m 

87.6 

17.8 

4.351 

1451 

Q.8281 

cantilever ;injection  flutter 

SW-6 

1.30 

31.2 

fga 

4  =  ,o 

15.4 

2.551 

189 1 

0.2771 

0.7991 

cantilever  legitimate 
flutter 

SW-7 

1.44 

30.0 

1350 

48.5 

17.2 

3.201 

1611 

0.3151 

0.8451 

c  antilever ; legitimate 
flutter 

SW-7a 

1.92 

24  .4 

1630 

79  .8 

25-5 

3.641 

1711 

0.2951 

■H 

pH 

t- 

co 

0 

free-to-roll  legitimate 
flutter 

Wings  with  Ailerons 


Model 

Hf 

®f 

cps 

vf 

ft/sec 

B 

vf 

vf 

CD 

Ql 

cps 

% 

““1 

Remarks 

m 

% 

““1 

SW-3b 

1.47 

114 

1370 

53.6 

4.^9 

3.55 

147.6 

0.237 

0.885 

1.002 

aileron  frequency 
taken  same  as  first 
torsion;  legitimate 
flutter;  very  little 
aileron  motion  in 
flutter  mode 

*  Injection  flutter 


+  Retraction  flutter 

1  Based  on  cantilever  frequencies  of  bare  wing  (see  Table  D,15  for  further  vibration  data) 

2  Coupled  aileron  frequency 
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TABLE  D. 8a 


EXPERIMENTAL  FLUTTER  DATA  FOR  SWEPT  WINGS  (Cont.) 


Bare  Wings 


Model 

vf 

vf 

Kis  •Ho  j/Mi  | 

1 

L'lSJIJ 

a<  JIMns 

sw-3 

8.38 

7-05 

0.179 

SW-3-1+ 

8.62 

8.89 

0.193 

sw-3c 

8.78 

7.14 

0.204 

SW-3d 

7-36 

6.24 

0.177 

SW-4* 

12.45 

9.29 

0.227 

sw-5* 

14 . 54 

10.02 

0.205 

SW-8* 

8.80 

6.80 

0.227 

SW-8-1 

7.69 

6.33 

0.196 

Wings  with  Tip  Tanks 

■B 

32.4 

7.91 

■H  i 

28.0 

4.64 

-- 

SW-7 

31.3 

5.82 

-- 

SW-7a 

46.4 

6.62 

-- 

Wings  with  Ailerons 

SW-3b 

8.34 

6.45 

-  0.204 

*  Injection  Flutter 
+  Retraction  Flutter 
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TABLE  D.9 


EXPERIMENTAL  FLUTTER  DATA  FOR  DELTA  WINGS 


Bare  Wings 


De-L  1.72- 
1.3? 


De-2a-l  1.80- 

l.?fi 

De-2b  1 .8',- 

1 .?"' 

De-2c  1.9?- 

1.26 


De-2d-l  1.62- 

1.26 


De-j-1  1.82- 

1.28 

De-3a  1.86- 

1.37 

De-3b  1.87- 

1.33 

De-3c  1.84- 

l.*3i 

De-3d  1.79- 

1.61 


De-jg  1.90 

De-3h  1.84- 

1.26 

De-3i  1.92 


De-’l  1.84- 

1.36 

De-4a  1.90- 

1.26 

De-4b  1.90 

De-4c  1.67 


I 

I 


Model  I  M 


vf 

v£ 

vf 

03 

°1 

cps 

n 

E9 

ft/ sec 

■ 

,Vf£ 

m 

““1 

05 

ai 

- 

86-72 

- 

- 

361 

0.346 

- 

1220 

40.9 

1.53 

1.49 

313 

E«t. 

0.307 

Esc. 

0.621 

Est . 

1400 

72.0 

1.55 

1.71 

313 

0.345 

0.891 

1630 

80.3 

2.39 

2.10 

297 

0.31  > 

0.875 

-- 

90-6? 

-- 

320 

0.4)9 

1.10 

- 

oe-o” 

- 

310 

0.4-n 

1.39 

-- 

98-34 

- 

- 

337 

0.475 

1.25 

98-49 

- 

- 

283 

0.^77 

1.27 

1290 

48.7 

3.57 

2.19 

225 

0.340 

1.22 

-- 

86-59 

-- 

-- 

240 

0.463 

1.23 

1260 

50.4 

2.53 

2.04 

236 

0.445 

1.22 

- 

93-55 

- 

- 

318 

0.712 

1.87 

-- 

93-67 

- 

-- 

320 

0.459 

1.24 

-- 

98-61 

- 

- 

362 

0.486 

1.29 

-- 

4_T 

1 

O 

- 

-- 

357 

0.448 

1.24 

-- 

84-68 

- 

-- 

302 

0.384 

1.06 

1390 

47,2 

3.64 

?.y2 

208 

0.385 

1.15 

1690 

80.4 

4.6o 

2.89 

210 

0.407 

1.21 

1620 

86.5 

3.60 

2.60 

238 

0.399 

1.08 

-- 

-- 

— 

-- 

286 

0.479 

1.22 

— 

.00-56 

- 

- 

333 

0.411 

1.15 

- 

-- 

- 

280 

0.421 

1.19 

1620 

i.eo 

2.45 

253 

0.415 

1.19 

— 

69-44 

-- 

-- 

310 

0.506 

1.25 

-- 

94-44 

- 

- 

314 

0.458 

1.23 

-- 

-- 

- 

- 

275 

0.484 

1.20 

1610 

-;>.i 

2.56 

2.50 

231 

0.437 

1.18 

13=,0 

'12.6 

2.59 

2.21 

233 

0.455 

1.22 

1690 

?6.0 

3.73 

3.23 

189 

0.247 

0.718 

1630 

113.5 

7.94 

129 

0.355 

0.868 

no  flatter 

retraction  flutter 

injection  flutter 
injection  flutter 
no  flutter 

no  flutter 

no  flutter 

no  flutter 

legitimate  flutter 
no  flutter 

same  as  De-2d-l  but  at  lower 
li;  legitimate  flutter 

no  flutter 

no  flutter 

no  flutter 

no  flutter 

no  flutter 

injection  flutter 
legitimate  flutter 
injection  flutter 

model  destroyed  by  start¬ 
ing  shock 

no  flutter 

model  destroyed  by  start¬ 
ing  shock 

injection  flutter 
no  flutter 

no  flutter 

model  destroyed  by  start¬ 
ing  shock 

legitimate  flutter 
legitimate  flutter 
injection  flutter 
injection  flutter  i 


tuf 

cps 

vf 

f  t / s  ec 

1 .92 

27 ‘j 

1620 

Wings  with  Elevons 


vf 


1.*I2  aileron  buzz 

throughout  run 


•Injection  flutter 
+Retraction  flutter 
(  1)  Coupled  aileron  frequency 
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TABLE  D.9a 

EXPERIMENTAL  FLUTTER  DATA  FOR  DELTA  WINGS  (Cont.) 

Bare  Wings 

Model 

vf 

vf 

b0.75 

“u  \{*\ 

b0.75  f 

b0.75  1 

af 

vyWo.75 

De-la+ 

6.12 

5.96 

--  - 

De-la-1* 

6.20 

6.84 

— 

De-la-2* 

9.56 

8.40 

— 

Be -2d 

14.28 

8.76 

0.132 

De-2d-2 

10.12 

8.16 

0.141 

De-3e* 

14,56 

11.68 

0,131 

De-3f 

18.40 

11.56 

0.172 

De-3f-l* 

14 . 40 

10.40 

0.211 

De-3j# 

7.20 

9.80 

0.164 

De-4c 

10.24 

10 . 00 

0.147 

De-4c-l 

10.32 

8.84 

0.133 

* 

j  De-5 

14.92 

12.92 

0.160 

^  * 

De-6 

31.76 

19.32 

0.132 

Wings  with  Elevons 

.....  .  j 

De-2e* 

9. '08 

11.08 

— 

*  Injection  Flutter 

+  Retraction  Flutter 
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TABLE  D.10 


TIP-TANK  PARAMETERS  FOR  STRAIGHT  WINGS 


Model 

m 

0 

slugs/ in 

mt 

( c  •  §  • )  p 

%  wing 
tip  chord 

<<  > 

- .  g .  ip 

dc  „  > 

^  s  o  1  .j 

calc . 
locus  of 
shear 
centers 
%  chord 

V 

ST 

(mb)_£ 

span 

EBUB 

ST-la 

0.00413 

■ 

56.0 

0.155 

42.5 

0.401 

ST-lc 

0.0044? 

-  3 

39.1 

0.149 

'  37.4 

0.0506 

ST-4a 

0.00406 

nu 

25.0 

0.133 

42.2 

-0.511 

ST-4b 

0.00406 

: 

77.5 

0.09T7 

42.2 

1.05 

ST-4c** 

0.00387 

1.40 

25*0 

0.131 

42.4 

-0.517 

ST-4c-l** 

0.0038? 

1 .40 

25.0 

0.131 

42.4 

-0.517 

ST-4c-2** 

0.00387 

1.40 

25.0 

0.131 

42.4 

-0.517 

S"T-4c-3** 

0.00387  • 

1.40 

25.0 

0.131 

42.4 

-0.517 

*about  calc,  locus  of  shear  centers 
**free-to-roll ,  Ig  =  0.00150  slug-ft2 


Geometrical  shape  -  body  of  revolution  obtained 
from  sector  of  circle 

Tip  tank  symmetrically  located  at  wing,  tip  chord 

Fineness  ratio  =0,10 

length  tip  tank  _  ~ 

"  wing  tip  chord  -3 
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TABLE  D. 11 


TIP-TANK  PARAMETERS  FOR  SWEPT  WINGS 


1 

mo 

slugs/ in 

Mx 

>5 

(c .  g .  )  j 

%  wing 
tip  chord 

(r  c.g!T 

'(’■c.g.h 

H,.(bT)S 

calc . locus 
of  shear 
centers 
$  chord 

5t 

ST 

(SB) 

.70  span 

SW-3a 

0.00406 

1.40 

25.0 

0.116 

42.0 

-0.51 

SW-6 

0.00409 

1.39 

25.0 

0.128 

42.0 

-0.51 

SW-7 

0,00391 

1.45 

25.0 

0.125 

42.0 

-0.53 

SW-7  ar* 

0.00406 

1.40 

25.0 

0.119 

42.0 

-0.51 

*free-to-roll,  Ig  =  0.00150  slug-ft2 


Geometrical  shape  -  body  of  revolution  obtained 
from  sector  of  circle 

Tip  tank  symmetrically  located  at  wing  tip  chord 

Fineness  ratio  =0.10 

length  tip  tank  _  _ 
wing  tip  chord”  ~  ^ 
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TABLE  D.12 


AILERON  PARAMETERS  FOR  STRAIGHT  AND'  SWEPT  WINGS 


Model 


ST-Ib 

ST-le 

ST-lf 

SW-3b 


Model 


Half-span  aileron  (outboard) 

Aileron  binge  line 

80.0%  of  wing 
chord  (no  aero¬ 
dynamic  balance) 

Ratio  of  aileron  mass  per  unit 
span  to  total  wing  mass  per  unit 
span 

0.100 

A»ileron  center  of  gravity  in  per 
cent  of  aileron  chord 

30.0% 

Aileron  radius  of  gyration  about 
aileron  binge  line  in  per  cent 
of  aileron  chord 

45.0% 

TABLE  D. 13 

ELEVON  PARAMETERS  FOR  DELTA 

WINGS 

Full-span  elevon.  constant  chord, 

|  no  aerodynamic  balance  | 

Elevon  hinge  line 

perpendicular  to 
root,  14.0%  of 
wing  root  chord 

Ratio  of  average  elevon  mass 
per  unit  span  to  average  total 
wing  mass  per  unit  span 

0.10 

Elevon  center  of  gravity  in  per 
cent  of  elevon  chord 

33% 

Elevon  radius  of  gyration  about 
elevon  hinge  line  in  per  cent  of 
elevon  chord 

45% 

CONFiDENTIAL 


CONFIDENTIAL 


TABLE  D.  14 

EXPERIMENTAL  VIBRATION  DATA  FOR  STRAIGHT  WINGS  j 

MODEL 

1  ST— 1 

ST— Id 

ST-id-H  I 

FIRST 

NATURAL 

FREQUENCY 

y//y  /  /  /  // 

TT 

0J-  47.1 

9 =.006 

01  =  42.0 

g  =  .008 

WoWWJa 

0)  =  46.7 

g  *.031 

SECOND 

NATURAL 

FREQUENCY 

W 

OJ  =  153 

g  *  .007 

W  =  140.5 

9=  .007 

^—j 

W  =  142 

g  =  .015 

THIRD 

NATURAL 

FREQUENCY 

U 

0J  =  181 

g  =  .007 

U 

0J  =  168 

g  =  .007 

u 

01=  173 

g  =  .013 

FOURTH 

NATURAL 

FREQUENCY 

0)=  300 

g  =  .019 

TO  NULL  AREAS  w  =  cps. 
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TABLE  D.I4  (CONT.) 


EXPERIMENTAL  VIBRATION  DATA  FOR  STRAIGHT  WINGS 


MODEL 


ST— 2 


ST-3 


ST- 4 


FIRST 

NATURAL 

FREQUENCY 


SECOND 

NATURAL 

FREQUENCY 


(V-  38.00 
9  =  .015 


0J=  34.38 
9  =  .OI8 


U)*  57.6 

9  «.OII 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 

NATURAL 

FREQUENCE 


Od=  HO 

g*  .012 


UJ  =  150 
9=  .010 


CU  =  235 
g  =  .0i2 


0J«  90.0 

g  ■  .oio 


OJ-  136 
9  =  .019 


NO  NODAL  LINES 
AVAILABLE 


GU=  204 
g  =  .008 


a;»  167 

g  =  .008 


NULL  AREAS 


d)  =  cps. 


CONFIDENTIAL 


TABLE  D.  14  (CONT.) 


EXPERIMENTAL  VIBRATION  DATA  FOR  STRAIGHT  WINGS 


ODEL  ST-4- 1 


ST — 5 


ST—5 —  I 


FIRST 

NATURAL 

FREQUENCY 


SECOND 

NATURAL 

FREQUENCE 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 

NATURAL 

FREQUENC' 


NULL  AREAS 


.013 


0)  =  c  ps. 
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CONFIDENTIAL 


TABLE  D.  14  (CONT) 


EXPERIMENTAL  VIBRATION  DATA  FOR  STRAIGHT  WINGS 


MODEL 


FIRST 

NATURAL 

rREGUENCY 


SECOND 

NATURAL 

FREQUENCY 


THIRD 

NATURAL 

FREQUENCY 


OJ*  140.9 

g  =.oii 


0J  =  128.3 
g  *.oi4 


Gtf*  113 

g  =.017 


FOURTH 

NATURAL 

FREQUENC 


oj  =  lei 
g  =  .ois 


Gl)  =  300 
g=.oi7 


NULL  AREAS 


cu*  162 
g  *  .019 


O)  =  cps. 


GU=  148 
g  * .014 


GU  =  220 
9  =.017 


WADC  Til  Part  II 


CONFIDENTIAL 


TABLE  D.  14  (CONT.) 


EXPERIMENTAL  VIBRATION  DATA  FOR  STRAIGHT  WINGS 


9=  .026 


CU=  241 


9  =  .017 


NULL  AREAS 


( i)  =  cps. 
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TABLE 

:  D.  14  {CONT.) 

EXPERIMENTAL  VIBRATION  DATA  FOR  STRAIGHT  WINGS 

MODEL 

ST-12 

ST- la 
CANTILEVER 

ST— la 
CANTILEVER 

FIRST 

NATURAL 

rREQUENGY 


U-  71.8 
g*  .024 


CJ  *  47.1 

g  =  .006 


(x) «  1 4.3 

g  *  .009 


SECOND 

NATURAL 

FREQUENCY 


CU-  228 
g=  .014 


0J  E  153 

g  =  .007 


Cl)  -  31.0 
9  .017 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 

NATURAL 

FREQUENCE 


OJ  =  272 

g=  .02  2 


CU  =  181 
g  -  .007 


GU  =  125 
3  *  .007 


NO  NODAL  LINES 
AVAILABLE 


Cl) *  2  35 
g  *  .004 


NULL  AREAS 


< j)  =  cp$. 
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CONFIDENTIAL 


CONFIDENTIAL 


CONFIDENTIAL 


TABLE  D.  14  (CONT.) 

EXPERIMENTAL  VIBRATION  DATA  FOR  STRAIGHT  WINGS 

! 

MODEL 

ST—  4c 
CANTILEVER 

ST- 4c 
CANTILEVER 

ST-4c 

FREE  TO  ROLL 

FIRST 

NATURAL 

FREQUENCY 

/////y  ///  y 

/'//////// 

■ - 

1 

i 

j 

i 

1 

0)  =  60.00 

g  =.oo7 

s 

(jU  =  16.46 

g.=  .013 

U)  =  100 

g  =.oo6 

SECOND 

NATURAL 

FREQUENCY 

w 

J 

ffl 

w 

GU=  175 

g  =  .008 

0)  =  39. * 

9  =  .014 

J 

0)  =  176 

g  =  .005 

THIRD 

NATURAL 

FREQUENCY 

U 

U 

01  =  238 

9-  .015 

(i)  =  144 

9  =.026 

CU  =  2  5  4 

9  =  .008 

FOURTH 

NATURAL 

FREQUENCE 

EM  NULL  AREAS  w  =  cps.  J 

EXPERIMENTAL  VIBRATION  DATA  FOR  STRAIGHT  WINGS 


MODE 


FIRST 

NATURAL 

FREQUENCY 


SECOND 

NATURAL 

FREQUENCY 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 

NATURAL 

FREQUENCY 


CU  -  242 

g  =  .005 


NULL  AREAS 


0)  =  cps 
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CONFIDENTIAL 


TABLE  D.l  4  (CONI) 


EXPERIMENTAL  VIBRATION  DATA  FOR  STRAIGHT  WINGS 


MODEL 


ST- 4c -2 
FREE  TO  ROLL 


FIRST 

NATURAL 

:REQUENCY 


SECOND 

NATURAL 

FREQUENCY 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 

NATURAL 

FREQUENCY 


CO  =  99.0 

g=.oio 


G0  =  174 

g=.oos 


(0  =  262 
g=.oo7 


NULL  AREAS 


^ x) *  34.8 

g  *.oie 


(0-  55.3 
g  =  .oi5 


GO*  172 

g  *  .008 


g  =  .oo6 


U)  =  c  ps. 


CO  =  o3.8 

g  =  .026 


C0  =  175 

g  =  .010 


C0  =  235 
g=.oi3 


CO  =  333 

g  =  ,oi4 
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CONFIDENTIAL 


TABLE 

D.  14  (CONT.) 

EXPERIMENTAL  VIBRATION  DATA  FOR  STRAIGHT  WINGS 

MODEL 

ST — lb 

AILERON 

ST—  Id 
AILERON 

ST— If 

AILERON  LOCKED 

FIRST 

NATURAL 

FREQUENCE 


WING  UNCLAMPED  I  WING  UNCLAMPED 


TMKTjTMWJWA 


CD?  49.07 

g  =  ,oio 


AILERON  MODE 


CD  *  52.6 
9  *  .029 


a.1  =  49.3 


6 


SECOND 

NATURAL 

FREQUENCY 


CD  s  1 34 
9  =.009 


CD  =  141 
9  = .007 


CD  =  147.2 
9  =*.013 


THIRD 

NATURAL 

:REQUENCY 


FOURTH 

NATURAL 

FREQUENCY 


CD  =  194 
9  =  .008 


Vi 

CD  =  190 

9  =  .008 

/ 

AILERON  MODE 

NO  NODAL  LINES  AVAILABLE 

Vi 

CD  8  210 

1 

1 

CD  =  192 
9  =  .0I4 


NULL  AREAS 


CD  =  Cps. 
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TABLE  D. i 5  (CONT.) 

EXPERIMENTAL  VIBRATION  DATA  FOR  SWEPT  WINGS 


W  =  CPS. 


CONFIDENTIAL 


table  d.  15  (CONI.) 

EXPERIMENTAL  VIBRATION  OATA  FOR  SWEPT  WIF 


FIRST 

NATURAL 

FREQUENCE 


SECOND 

NATURAL 

FREQUENCY 


SW-2b-l 


W=  31.2 
9  =  .007 


(a)  =  109 
9  «  .006 


SW-2b-2 


U)  =  31.2 
9  =>.015 


cu=  III 

g 


SW-3 


(a)  =  33.5 
9  = .020 


Ct>  «  119 

9  =.013 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 

NATURAL 

FREQUENCY 


a>  =  137 
g  =  .oo5 


O):  214 
9  =.004 


NULL  AREAS 
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U)  =  Cp8. 


(a)  •  132 
S  « .015 


CONFIDENTIAL 


a 


TABLE  D.  15  (CONT.) 


EXPERIMENTAL  VIBRATION  DATA  FOR  SWEPT  WINGS 


FIRST 

NATURAL 

FREQUENCY 


0)*  27.2 
9  *.024 


6 U*  28.4 
g  *.024 


OU  =  37.86 
g  *.oi3 


SECOND 

NATURAL 

FREQUENCY 


U)s  97.6 
9  *.022 


OJ  =  100.0 

g  *.oi8 


OJ  *  128 
9  *.024 


THIRD 

NATURAL 

FREQUENCY 


OJ*  121 
g  *,03i 


tU  *  113 
9  *  .022 


OJ  *  150 
9  *  .Oil 


FOURTH 

NATURAL 

FREQUENCY 


FIFTH 

NATURAL 

FREQUENCY 


CU  =  231 

CU*  237 

g  =  .020 

g  =.03i 

W*  307 
9  *  .024 


NULL  AREAS 


OJ  =  cps. 
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TABLE  D.  15  (CONI) 


EXPERIMENTAL  VIBRATION  DATA  FOR  SWEPT  WINGS 


SW-8-1 

SW-3o 

SW-30 

CANTILEVER 

CANTILEVER 

FIRST 

NATURAL 

FREQUENCY 


CD*  37.18 
9  *  .015 


CD  »  34.78 
Q  *.025 


CD  *  9.93 

g  •  .oio 


SECOND 

NATURAL 

^EQUENCY 


yys.'ysssji 1 


CD*  125 
Q  *.028 


CD*  120 

g  * .012 


CD  =  30.88 

g  *.oio 


i  j 


THIRD 

NATURAL 

REQUENCY 


CD*  146 
g  *.014 


CD  *  145 

g  *.oio 


CD*95.7 
g  *.013 


FOURTH 

NATURAL 

FREQUENCY 


Ci)  *241 
g  =.034 


CD*  237 

g  *.018 


( D  *  200 

g  *  .007 


l'/.VA'.V 


NULL  AREAS 


CD*  295 

g  *.010 


Cl)  1  cps 
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TABLE 

D.  15  (CONT) 

EXPERIMENTAL  VIBRATION  DATA  FOR  SWEPT  WINGS 

MODEL 

SW-6 

SW-6 

SW-7 

CANTILEVER 

CANTILEVER 

CANTILEVER 

FIRST 

NATURAL 

FREQUENCY 


SECOND 

NAfURAL 

FREQUENCY 


CD *  52.4 
9*.027 


U )«  151 

g  =  .050 


CD  =  15.79 

g  =  .021 


(D=  42.1 
g  =  .047 


CU  =  50.7 
g».030 


01  =  136 

g « .030 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 

NATURAL 

FREQUENCY 


FIFTH 

NATURAL 

FREQUENCY 


Ci)  •  189 
g  =.009 


NULL  AREAS 


CD  ■»  115.2 
g  =  .014 


fjJ  =  CpS. 


CD  =  161 

g »  .012 


CD=  273 
9  « .014 
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CONFIDENTIAL 


TABLE  D.I5(C0NT.) 


EXPERIMENTAL  VIBRATtON  DATA  FOR  SWEPT  WINGS 


MODEL 


FIRST 

NATURAL 

FREQUENCY 


SECOND 

NATURAL 

FREQUENCY 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 


SW-7 

CANTILEVER 


SW-7a 

CANTILEVER 


01  =  12.40 
5  .=  .012 


G i)E  50.5 
9  ».027 


0J=  105 
9 =.026 


01*  171 
9=  .Oil 


WADC 


$W--7a 
CANTI  LEVER 


ai  » 12.47 
9  *  .013 


01  =  105 

g  =  .015 


M 

\A 

01 »  2  33 

C 0 =233 

•O 

01=  220 

9  =.016 

9  =  .029 

g  =  .009 

NULL  AREAS 

cu  =  cps. 

L13j  Part  II 
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TABLE 

D.I5  (CONT.) 

EXPERIMENTAL  VIBRATION  DATA  FOR  SWEPT  WINGS 

MODEL 

SW”7g 

FREE  TO  ROLL 

SW-7a 

FREE  TO  ROLL 

SW— 3b 
AILERON 

FIRST 

NATURAL 

FREQUENCY 


CD  *73.70 
9  =.011 


CD  =  29.75 
9  =.011 


fAV»V.’ 


CD  s  35.0 
9  =.021 


SECOND 

NATURAL 

FREQUENCY 


THIRD 

NATURAL 

FREQUENCY 


:U-  175 

g=.oo6 


0D=  231 
9  =.007 


•v.vvavi 


CD  *  53.3 
9  =  .011 


CD  =  125 
9  =.013 


CD  =  1 30.6 
9  =  .014 


AILERON 

FREQUENCY 


a)  =  147.6 
9  =  .015 


FOURTH 

NATURAL 

FREQUENCY 


NULL  AREAS 


CD  =  236 
9  =.009 


CD  =  CPS. 


CD  =  280 
9  =.015 
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TABLE  D.  16  (CONT.) 

EXPERIMENTAL  VIBRATION  DATA  FOR  DELTA  WINGS 


MODEL 


FIRST 
NATURAL 
FREQUENCE 


SECOND 

NATURAL 

FREQUENCY 


la  -2 


THIRD 

NATURAL 

FREQUENCY! 


FOURTH 

NATURAL 

FREQUENCY 


U)=  260 
9  =.03  4 


D®  2 


ie-  2aH 


a )- 146.1 
g  =  .032 


W  297 
9=  .019 


1  A 

\ 

W=3$l 

g  =  -  - 


GU  =  580 

g=  -  - 


Cl)  =  430 
9*  .05 


Ul- 640 
9' .05 


NULL  AREAS 


(t)  -  cps. 
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TABLE  D.iS  (CONT.) 


EXPERIMENTAL  VIBRATION  DATA  FOR  DELTA  WINGS 


MODEL  De-2b 


FIRST 

NATURAL 

FREQUENC 


CONFIDENTIAL 
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TABLE  D.  16  (CONT.) 


EXPERIMENTAL  VIBRATION  DATA  FOR  DELTA  WINGS 


MODEL 


De  —  3i 


FIRST 

NATURAL 

FREQUENCY 


SECOND 

NATURAL 

FREQUENCY 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 

NATURAL 

FREQUENCY 


NULL  AREAS 


(d  =  cps 


CONFIDENTIAL 


CONFIDENTIAL 


|  TABLE 

D.  IS  (CONT.) 

EXPERIMENTAL  VIBRATION  DATA  FOR  DELTA  WINGS 

MODEL 

59 

<<$• 

1 

» 

O 

De  —  4b 

De-4c 

FIRST 

NATURAL 

FREQUENCT 


•  SECOND 
NATURAL 
FREQUENCY 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 
NATURAL | 

frequency! 


OJ=  144 
g=.036 


OJ=  385 
9=.026 


CU*  133 
5=  .036 


W*  330 
9=.034 


0)  =  101 
g  =  — 


W  =  272 
9 


NULL  AREAS 


U)  =  cps. 


WADC  TR  54-113  Part  II 


CONFIDENTIAL 


TABLE  D.I6  (CONT.) 


EXPERIMENTAL  VIBRATION  DATA  FOR  DELTA  WINGS 


MODEL  De-4c- 


De-  5 


De— 6 


FIRST 

NATURAL 

FREQUENCY 


SECOND 

NATURAL 

FREQUENCY 


THIRD 

NATURAL 

FREQUENCY 


FOURTH 

NATURAL 

FREQUENCY 


NULL  AREAS 


0)  =  CDS. 
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FIGURE  D.7  LOCATION  OF  INFLUENCE-COEFFICIENT  POINTS  FOR  STRAIGHT 

WING  MODELS 
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Units  are  inches/# 
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WING  MODEL  WITHOUT  ELEVON 


EXPERIMENTALLY  DETERMINED  FORCE-DEFLECTION 
INFLUENCE- COEFFICIENT  MATRIX  FOR  MODEL  De-la- 
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Units  are  Inches/ 


EXPERIMENTALLY  DETERMINED  FORCE -DEFLECTION 
INFLUENCE- COEFFICIENT  MATRIX  FOR  MODEL  De-3f- 
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